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1.0  INTRODUCTION 


Responding  to  a  well  recognized  need  by  many  in  both 
the  modeling  and  operational  communities  for  an  im¬ 
proved  capability  in  the  collection  and  assessment  of 
whole  sky  cloud  characteristics,  a  new  generation  of 
video  based  imaging  systems  was  developed  and  fielded 
by  the  Marine  Physical  Laboratory  (MPL)  of  Scripps 
Instituticffi  of  Oceanography  (Johnson,  et  al,  1989).  This 
Final  Report  summarizes  the  development  and  use  of  a 
resultant  data  base,  generated  by  these  previously  fielded 
systems,  containing  approximately  900  Gigabytes  of 
digitized  raw  imagery.  These  images  represent  approxi¬ 
mately  4600  data  days  geographically  distributed  among 
seven  widely  separated  CONUS  sites. 

There  were  several  inter-related  tasks  associated  with 
the  program  pursued  during  the  forty-six  months  covered 
in  this  Final  Report.  In  their  simplest  form  they  were  to, 
a)  provide  technical  coordination  and  operational  main¬ 
tenance  for  the  several  on-site  video  systems,  b)  establish 
an  image  oriented,  data  management  and  archival  facility 
at  the  Marine  Physical  Laboratory,  c)  conduct  a  cloud 
image  analysis  program  directed  toward  the  determina¬ 
tion  of  Cloud-Free  Line  of  Sight/Cloud  Free  Arc  (CFLOS/ 
CFARC)  probability  distributions  and  the  refinement  of 
analytical  models  for  estimating  these  distributions,  and 
d)  initiate  preliminary  studies  into  the  development  of 
experimental  procedures  to  enable  temporally  short  term 
predictions  of  those  cloud  field  dynamics  associated  with 
CFLOS/CFARC  applications  to  the  GBL  mission.  Each 
of  these  tasks  has  been  addressed  in  accordance  with  the 
programs  on-going  adjustment  of  goals  and  priorities, 
and  is  discussed  separately  in  the  following  sections  of 
this  Final  Report. 

2.0  AUTOMATED  SYSTEMS  FOR  CLOUD 
ASSESSMENTS 

The  Whole  Sky  Imager  (WSI),  as  an  emerging  opera¬ 
tional  system  has  been  described  previously  in  GL-TR- 
89-0061,  Johnson,  et  al,  1989.  In  its  final  field  ready 
configuration,  reviewed  briefly  in  this  introductory  sec¬ 
tion,  the  WSI  system  provided  the  data  acquisition  mecha¬ 
nism  for  generating  the  image  oriented  data  base  which 
is  discussed  in  the  body  of  this  final  report. 

2.1  The  Whole  Sky  Imager 

The  Whole  Sky  Imager  (WSI)  is  a  ground-based 
electronic  imaging  system,  wliich  monitors  the  upper 
hemisphere.  It  is  a  passive,  i.e.  non-emissivc  system, 
which  acquires  calibrated  multi-spcctral  images  of  the 
sky  dome.  The  system,  shown  in  Fig.  2-1,  views  the  sky 
through  a  scries  of  spectral  and  neutral  density  filters. 


Fig.  2-1.  Whole  Sky  Imager  Camera  Assembly 


using  a  fisheye  lens  to  acquire  most  of  the  sky  dome.  A 
fixed  gain  CID  (charge  injection  device)  solid  state 
camera  is  utilized,  and  a  full  set  of  radiometric  and 
geometric  calibrations  are  acquired  prior  to  fielding  the 
system.  Data  are  acquired  in  512  x  480  format,  which 
yields  1/3  degree  spatial  resolution.  This  corresponds  to 
a  17  meter  footprint,  for  a  cloud  layer  at  3  km  height. 

In  the  field,  cloud  images  are  acquired  under  control 
of  an  IBM  AT-class  microcomputer.  This  is  a  stand¬ 
alone  unit,  requiring  essentially  no  user  intervention; 
control  of  all  peripheral  functions  is  fully  automatic. 
Four  digital  images  are  acquired  every  minute,  and 
archived  on  8  mm  tape,  for  post-processing.  Approxi¬ 
mately  1 . 2  gigabytes  are  acquired  and  archived  per  week 
at  each  site. 

A  sample  radiance  image  is  shown  in  Fig.  2-2.  In  this 
image,  the  center  is  at  the  zenith  overhead,  and  the  edges 
are  near  the  horizon.  The  black  square  is  the  sun  occultor, 
which  shades  the  lens  in  order  to  minimize  stray  light.  A 
tower  may  also  be  seen  in  the  field  of  view  near  the 
occultor.  In  the  processing,  a  determination  is  made  at 
each  point  in  the  image,  yielding  a  cloud  decision  image 
at  full  spatial  re.solution.  The  cloud  decision  image  is 
illustrated  in  Fig.  2-3.  In  this  illustration  the  areas 


ri}!.  2-2.  Siunplc  WSI  Image,  acquired  at  6.50  nm 


I-'ig.  2..<.  Sanqile  CI'and/No  cloud  Decision  Image. 
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d  :::  1  'I'i'jMqiie  climd  are  light  gre>  through  w  hile  respec- 
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Fig.  2-4.  Sample  One-Minute  CIoud/No  cloud 

Decision  Image.  This  shows  the  subset  saved 
1  minute  intervals. 


In  summary,  the  Whole  Sky  Imager  is  a  system  for 
cloud  a.s.se.ssmeni.  specifically  for  the  archival  of  cloud 
Held  spatial  and  temporal  dynamics.  The  WSI  acquires 
calibrated  multi -spectral  sky  radiances  every  minute. 
From  thc.se,  automatic  cloud  determinations  arc  pro¬ 
duced  with  l/.l  degree  spatial  resolution.  The  system  was 
originally  deployed  at  7  sites  throughout  the  country, 
illustrated  in  Fig.  2-5  and  as  of  1  Jan  9 1 .  had  acquired  an 
average  of  26  months  of  data  per  site.  Tlte  data  will  be 
furtlier  discussed  in  Section  3. 


•  WSMRCiNMi  •  KAF0I^Ml  •  WSMRI-t  ,SMl  •C..SWC.CA 

•  UATBiMTl  •  COLUMBlAtMO;  •WALABARtFl' 

Fig.  2-5.  Whole  Sky  Imager  data  sites 

2.2  The  Portable  WSI 

7hc  portable  \V  SI  unit  is  show'n  in  Fig.  2-6.  This  unit, 
designed  for  short  term  deployments  and/or  special  ap¬ 
plications,  is  simi  lar  to  the  WSI ,  but  has  a  more  transptrrt 
able  sup[xm  and  environmental  control  housing.  The 
portable  system  was  initially  inst.allcd  at  the  I’niversiiv 


of  Wisconsin,  for  a  cooperative  study  with  Eloranta  and 
Gmnd,  using  their  High  Spectral  Resolution  lidar  (HSRL) 
lidar  system.  The  intent  was  to  acess  techniques  to  relate 
the  WSI  cloud  determinations  to  the  cloud  optical  depth 
determined  by  the  lidar  system. 


Fig.  2-6.  Portable  Whole  Sky  Imager 


2.3  Supporting  Instrumentation 

The  WSI  cloud  data  base  is  stored  in  digital  format  on 
8mm  video  tape  cassettes  (Sony  P6-120MP,  NTSC). 
These  small  cassettes,  when  used  with  the  EXABYTE 
8200  streaming  tape  sub-system,  can  hold  up  to  2.2 
Gigabytes  of  digital  information  each.  All  WSI  data 
manipulations  conducted  at  MPL  are  built  around  this 
small  compact  tape  system  under  the  control  of  an  IBM/ 
AT  class  microcomputer.  All  software  is  MS-DOS 
compatible,  and  is  generally  written  in  either  FORTRAN 
or  C  programming  language. 

The  basic  hardware  configuration  required  for  WSI 
image  processing  is  essentially  the  same  as  that  required 
for  the  initial  acquisition.  That  is,  a  system  able  to  read 
the  8mm  data  tapes,  an  image  processing  system  able  to 
handle  the  512  X  480  pixel  image  format,  and  a  host 
computer  suitable  for  overall  control.  At  MPL,  the  data 
acquisition  system  shown  in  Fig.  2-7  is  also  appropriate 
for  alt  subsequent  image  processing  tasks.  For  bulk 
processing,  however,  the  basic  computer  unit  is  supple¬ 
mented  with  an  Eighteen  Eight  Labs  PL  1 250  floating 
point  array  processor  to  accelerate  the  system  through¬ 
put.  and  dual  EXABYTE  drives  are  attached  to  each 
control  computer.  The  complete  WSI  system  is  illus¬ 
trated  in  Fig.  2-6,  which  combines  the  photographic 
description  with  the  basic  system  block  diagram. 

3.0  THE  WSI  IMAGE  DATA  BASE 

The  WSI  cloud  data  base  was  acquired  primarily  to 
provide  a  reliably  calibrated,  computer  friendly  image 


library  that  could  be  used  in  a  variety  of  model  validation 
and  model  development  activities.  Of  specific  interest 
was  its  applicability  to  the  generation  of  joint  multi-site 
cloud  free  line  of  sight  (CFLOS)  and  cloud  free  arc 
(CFARC)  probabilities.  Tliis  interest  was  in  turn  driven 
by  concerns  related  to  CFLOS/CFARC  impacts  upon 
ground  based  laser  (GBL)  applications.  Within  this 
context,  it  was  important  that  the  imagery  within  the  data 
base  provide  computer  access  to  a  "micro-climatology" 
of  cloud  field  characteristics,  well  defined  in  space  and 
time,  and  at  relatively  high  spatial  and  temporal  resolu¬ 
tions.  The  fully  automatic  WSI  system  provided  this 
acquisition  capability. 

3.1  Background 

As  can  be  deduced  from  the  entries  in  Table  3.1,  not 
all  field  systems  produced  equal  shares  of  the  resultant 
data  base.  On  average,  field  system  downtimes  were  at 
about  a  20%  rate,  ranging  from  a  low  of  about  14%  at  the 
Columbia  NWS  site  to  an  unfortunate  high  of  about  28% 
at  the  China  Lake  NWC  site.  For  the  most  part,  the  no¬ 
data  intervals  at  China  Lake  were  due  to  power  outages 
that  exceeded  the  capacity  of  the  WSI  back-up  battery 
pack.  Early  in  the  data  collection  interval,  continuing 
problems  related  to  the  8mm  tape  subsystem  and  the 
video  frame  grabber  were  experienced.  Both  problems 
were  eventually  resolved  by  rcho-fitting  additional  cool¬ 
ing  and  air  filtering  capacity  to  the  interior  computer 
housing.  Air  conditioner  failures  in  the  external  housing 
were  a  continuing  irritation  at  the  high  humidity  sites, 
leading  to  their  eventual  replacement  with  a  chilled 
liquid  coolant  system.  The  efficacy  of  these  basic  system 
upgrades  was  reflected  in  the  ai'cragc  months  per  main¬ 
tenance  visit  which  improved  from  an  early  two  to  three 
month  average  interval  to  a  four  to  five  month  interval 
during  the  later  portions  of  the  program. 

3.2  Current  Extent  of  the  Data  Base 

The  Whole  Sky  Imager  (WSI)  data  network,  while 
fully  implemented,  consisted  of  seven  independently 
operating,  fully  automatic  electronic  camera  systems. 
Each  system  collected  whole  sky  imagery  for  twelve 
hours  a  day,  seven  days  a  week.  Each  resultant  weekly 
data  tape  currently  in  archive  contains  up  to  1 200  mega¬ 
bytes  of  image  data,  plus  around  200  megabytes  of 
housekeeping  space.  Thus,  at  the  end  of  the  data  collec¬ 
tion  interval  31  Dec  90,  the  WSI  cloud  data  base  con¬ 
tained  approximately  900  Gigabytes  of  raw  image  data. 
About  675  Gigabytes  of  these  data  have  been  processed 
through  varying  portions  of  the  sequence  discussed  in 
Section  4.  The  data  quantities  are  substantial,  and  have 
required  carefully  conceived  and  efficiently  executed 
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Table  3.1 


WSI  DATA  TAPE  STATUS  AS  OF  31  JANUARY  1991 

C-STA 

30  MAR  88 

HELSTF 

29  MAR  88 

KIRTLANO 

17  MAY  88 

CHINA  LAKE 

23  JUN88 

MALMSTROM 

29  AUG  88 

MALABAR 

18  NOV  88 

COLUMBIA 

9  FEB  89 

A 

TOTAL  DAYS  IN  FIELD 

995 

1008 

953 

897 

478 

735 

691 

B 

TOTAL  DAYS  OF  DATA 

763 

819 

749 

641 

406 

602 

592 

C 

TOTAL  FULL  DATA  DAYS 
(12  HR) 

652 

623 

607 

453 

332 

515 

561 

D 

TOTAL  PARTIAL  DATA  DAYS 

111 

196 

142 

188 

74 

87 

31 

E 

TOTAL  DAYS  DOWN 
(NO  DATA) 

232 

189 

204 

256 

78 

133 

99 

DATE  WITHDRAWN 

FROM  SERVICE 

19  Dec  90 

31  Dec  90 

25  Dec  90 

6  Dec  90 

20  Dec  89 

19  Nov  90 

31  Dec  90 

A  =  +  E 

B  =  C  +  D 

processing  procedures  to  insure  optimum  downstream 
uiilization.  Elements  of  these  procedures  are  summa¬ 
rized  in  Section  4.0. 

The  WSI  cloud  data  base  presently  resides  in  the 
custody  of  the  Marine  Physical  Laboratory  ( MPL),  Scripps 
Institution  of  Oceanography,  San  Diego.  C  A,  92093.  As 
of  1  Jan  91 ,  imagery  representing  4614  days,  geographi¬ 
cally  distributed  as  shown  in  Table  3.1  and  Fig.  3-1  have 
been  collected  and  returned  to  MPL  for  processing  and 
analysis.  A  reasonable  estimate  is  that  lliese  raw  data 
reside  on  approximately  900  data  tapes. 

The  preparation  of  these  900  raw  data  tapes  plus  their 
appropriate  output  products  with  the  procedural  docu¬ 
mentation  deemed  necessary  for  subsequent  usages  is 
essentially  the  task  outlined  in  Section  4.0. 

It  should  be  noted  that  these  data  quantities  do  not 
include  several  relatively  small  case  study  sets  collected 
by  the  portable  system.  The  data  listed  in  Table  3.1  and 
discussed  in  this  report  arc  only  those  assoc  iated  with  the 
multi-station  network  sponsored  under  Air  Force  con¬ 
tract  F19628-88-C-0005.  An  inventory  ol  the  complete 
WSI  Image  Database,  segregated  by  site,  date  and  archi¬ 
val  tape  number  is  included  as  Appendix  A.  The  data  in 
Appendix  A  have  been  reviewed  in  accordance  with  our 
established  QC  procedures,  and  have  been  purged  of  all 
sets  deemed  inappropriate  for  further  application.  The 
material  included  in  Apfiendix  A  has  been  previously 
released  as  Optical  Systems  Group,  Technical  Note  No. 
228,  M.  L.  Ciandro.  1991. 
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4.0  DATA  PROCESSING  PROCEDURES 

The  general  concept  guiding  the  post-acquisition  pro¬ 
cessing  of  the  WSI  data  base  is  illustrated  in  Fig.  4-1, 
"WSI  BASIC  IMAGE  PROCESSING  FLOW  CHART  . 
This  conceptual  game  plan  is  implemented  through  a 
series  of  proprietary  software  routines  designed  to  im¬ 
pose  strict  quality  and  calibration  controls  upxtn  the  raw 
data,  while  simultaneously  reducing  the  number  of  im¬ 
ages  required  for  insertion  into  the  optimized  cloud/no 
cloud  decision  processes.  Carried  to  its  nominal  conclu¬ 
sion,  this  processing  plan  reduces  the  basic  imagery  (that 
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Fig.  4-1.  WSI  Basic  Image  Processing  Flow  Chart 


created  by  each  camera  in  the  system)  to  about  one  fourth 
of  its  original  volume  at  the  composite  ratio  stage  and 
then,  if  desired,  by  an  additional  increment  depending 
upon  the  output  format  specified  by  the  user  of  the 
derived  products.  The  cloud/no  cloud  images  arc  created 
at  the  same  eight  bit  resolution  as  arc  the  raw  radiance  and 
composite  ratio  images. 

4.1  Data  Proces.sing  Procedural  Sequence 

With  the  basic  data  tapes  in  hand,  the  proces.sing 
hardware  and  software  on-line,  and  the  data  base  man¬ 
agement  system  implemented,  one  can  address  the 
stepwise  manipulation  of  all  component  elements  from 
raw  tape  to  delivered  end  product.  In  the  following 
paragraphs  this  stepwise  process  will  be  described  on  the 
basis  of  a  single  data  tapes  .sequential  manipulation,  and 
from  this  unit  tape  estimate,  the  overall  program  time 
estimates  will  be  produced. 

In  assessing  the  following  estimates  of  the  proces.sing 
times  associated  with  each  procedural  sequence,  one 
should  be  aware  that  throughput  times  arc  often  heavily 
influenced  by  different  machinc/softwarcinterfaccscvcn 
when  determined  on  systems  containing  nominally  "iden¬ 
tical "  hardware  and  code  versions  With  this  caveat  in 
mind,  the  estimates  included  herein  and  indicated  on  Fig. 
4-2,  were  derived  from  actual  full  length  processing 
episodes  rumting  in  the  hardware/softwarc  configuration 
listed  in  Table  4.1. 

4.2  Data  Processing  Run  Time  Estimates 

Processing  the  complete  stcpwi.se  procedure  required 
to  convert  a  field  tape  confining  seven  days  of  raw 
unedited  cloud  imagery  into  a  fully  processed  tape  con¬ 


taining  the  same  seven  days  of  data  in  the  form  of  cloud- 
no  cloud  imagery  is  illustrated  in  Fig.  4-2,  and  the  run¬ 
time  summary  for  both  machine  time  and  supporting 
analy.st  time  is  shown  in  Table  4.2. 
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Table  4.1 

WSI  Standard  Data  Processing  Configurations 


HARDWARE 

SOFIAV'ARE 

TM!  Model  20(11  A  Computer 

Proprietary  Software  Programs 

I.  S()386  CPf 

1.  TAPE(}C,  Vcrsion.xxx 

2,  803H7  Coprocessor 

2.  EXCOPY,  Version  xxx 

x  PL  1250  Array  Processor 

.1  TAPRAT,  Version  xxx 

1.  ASC-K8  SCSI  interface 

4.  CLDDEC.  Version  xxx 

5,  EXABYTE  8200. 2  ea 

\^/4S25  PROM 

t).  Seagate  60  .MB  Hard  Disc 

Table  4.2 

Raw  to  Decision  Conversions 
Run  Time  Summary 


■ 

Sciii.crce  No, 

.Machine  Hours 

A-oalyst  Hours 

iFromFig,  4-2) 

- 1 

Primary 

File  Copy 

Total 

Housekeeping  AOveisight 

1 

7 

2 

9 

2 

2 

5 

.5 

s  5 

0.5 

-3 

D 

,5 

25 

25 

4* 

iNT 

,5*5 

In 

05+0,5 

2‘2 

2+.2 

44 

05+05 

6  1 

1 

!  2 

05 

1 

5-6 

5+,5 

1.' 

0.5*05 

»•  1 

!  +  l 

.2+.2 

24 

0,5+05 

0 

1 

'•S 

1 2 

0,5 

f'llals  T!0 

5.8 

52 

(a: 

10 

1  or  10 

V 

L_il_ 

4!  S 

8 

'  S  '" ;ra;c  r',irs  rcc'u'rcii  !('r  cath  on:  minuir  or  ici'  irr.ouic  uragc  sa 
r  :•  arc  hasa!  upon  "clear.  '  r..nc  'Aiihoui  appreciable  sofpearcor 

■I'lrJctarc  p'::c.Vc _ 


In  cslimaiing  the  overall  processing  task,  as  outlined 
‘  e  4  2  &  Table  4.2,  it  is  important  to  note  that  the 
i.ii.  I.r  shoun  in  Table  4  2  are  really  lor  best  elfort  data 
jiioccssmg”  only.  They  do  not  include  the  very  substan- 
;i  1  i.'ll'mi  'Ahich  is  heavily  related  to  "housekeeping  and 
('.  rsighl".  but  which  is  more  specilically  onented  to- 
.'.ard  on-line  data  "interpretation  and  analysis"  (I&.A).. 
\('rdo  they  include  the  more  styli/cd  but  highly  impor¬ 
tant  implementation  and  maintenance  of  the  appropnate 
Ikita  Hase  Management  (unctions. 

rile  I'.iblc  4  2  comment  related  K'  "clean  "  mns  is  a 
tommeiu  igni'ietl  at  extreme  pcnl.  One  must  keep  in 
mind  dial  .idol  'A'SI  llcid  svsiems  'acre monitorcif  b\  on¬ 


site  host  peiMinnel  with  wadely  van  ing  degrees  of  tech¬ 
nical  background  and  program  related  interests.  Not  all 
received  the  imposition  ol  the  WSI  system  into  their 
daily  routines  with  equal  e.xuberancc.  The  vagane.s  of 
local  weather  extremes,  power  lluctuations.  harda  ire 
faults  and  conflicting  on-site  pnonties  all  contnbuted  to 
sometimes  substantial  variations  in  the  normally  full 
automatic  mode  that  the  WSI  attempted  to  maintain. 
With  this  "real-world"  context  in  mind,  it  is  easy  to 
visualize  the  nature  of  the  on-line  data  interpretation  and 
analysis  that  must  accompany  the  more  fonmalized  data 
"processing"  outlined  in  Fig.  4-2. 

The  primary  tools  employed  by  the  analyst  to  diag¬ 
nose  and  sanitize  those  raw  data  tapes  which  fail  either 
thcirinitialTAPEQCruns.orlheirsubsequentTAPRAT 
runs,  arc  the  regularly  produced  TAPEQC  output  files, 
the  Local  Apparent  Noon  (LAN)  image  files,  and  the  raw 
data  tapes  themselves,  all  in  concert  with  the  analysts 
highly  honed  understanding  of  the  WSI  hardware  and 
software  operating  regimes.  These  diagnostic  skills  are 
not  readily  transferable  and  they  are  absolutely  essential 
to  the  orderly  and  efficient  data  flow  shown  in  Fig.  4-2, 

Some  estimate  of  the  additional  machine  time  re¬ 
quired  to  support  the  "interpretation  and  analysis"  (l/.-\) 
function  plus  the  aitcndai.i  requirement  for  data  re-runs 
should  be  provided  as  pan  of  the  overall  ta.sk  definition. 
As  a  first  approximation,  one  finds  about  7%  of  the  data 
designated  as  special  handling.  During  the  first  quaner 
of  calendar  1990,  an  estimated  10%  of  the  accumulated 
data  runs  were  repicaied  due  to  spcci  fic  I/A  decisions  and 
another  10%  were  rc-runduc  to  mi.scellaneous  hardware 
peculiarities,  primarily  EXABYTE  incompatibilities. 
Thus  one  might  assume  a  worst  case  requirement  for  up 
to  27%  of  miscellaneous  re- runs,  however  as  "production 
line'  processing  settles  dow  n,  a  lOto  15%  re-runestimatc 
seems  more  reasonable. 

Descriptive  summanes  of  the  software  programs 
shown  in  Fig.  4-2,  and  tlicir  primary  sub-routines  are 
contained  in  .Appendix  B.  Piogram  listings  and  addi¬ 
tional  illustrative  matcnaN  as  are  required  by  the  ('on- 
tract  Data  Requirements  List  arc  provided  under  separate 
cover 

4.3  Data  Processing  .Status 

As  noted  in  Section  2. 1 ,  there  are  approximately  9(K) 
raw  data  tapes  currently  in  archive  status,  representing 
over  4fi()0  separate  data  days  Some  ot  these  data  tapes 
have  been  processed  ihnnigh  all  ot  the  sequences  illus¬ 
trated  in  Fig.  4-2,  however  the  majority  have  been  through 
steps  one  and  two  only.  .An  illustration  of  the  data  base 


status  as  of  30  Sep  91  is  shown  in  Fig.  4-3.  A  slightly 
more  detailed  representation  broken  down  by  major 
software  program  is  shown  in  Table  4.3. 
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Table  4.3 

Processed  Data  Tape  Summary 


Data 

1  talion 

Raw 

Data 

Tapes 

Major  SoRware  Programs 

TAPEQC 

CHECKLAN 

TAPI 

1 

lAT 

10 

DEC 

1 

TMF 

DEC 

10 

C  Sla 

190 

190 

70 

70 

70 

70 

70 

llc'^f 

166 

166 

0 

Kjiland 

142 

142 

69 

69 

69 

69 

69 

China  Ijke 

135 

135 

0 

.W  Imslrorn 

71 

71 

32 

32 

32 

32 

32 

.MTabar 

104 

IM 

27 

27 

27 

27 

27 

Ccjmbia 

95 

95 

59 

59 

60 

59 

60 

Teals 

903 

903 

257 

257 

258 

257 

258 

The  processing  sequence  illustrated  in  Fig.  4-1  is 
descriptive  of  the  preliminary  run  sequence  used  during 
the  mid- 1990  processing  of  the  "one-minute"  data  sub¬ 
set.  The  output  product  for  this  preliminary  sequence  is 
a  data  base  containing  derived  cloud/no  cloud  images 
ba.sed  upon  a  fixed  threshold  decision  algorithm  (Apr. 
■90). 

A  major  activity  during  late  1990  and  early  1991  has 
been  the  development  of  an  improved  cloud  decision 
algorithm.  This  effort  has  produced  a  new  composite 
decision  program  called  CMPDEC.  The  algorithm  im¬ 
proves  the  thin  cloud  discrimination  by  correcting  for 
aerosol  scattering  processc.i,  as  described  in  Section  5. 
Thcm'''.cd  is  particularly  effective  at  reducing  the  over¬ 
determination  of  thin  cloud  found  in  regions  near  the  sun 
and  close  to  the  horizon  that  was  made  by  the  earlier  fixed 
threshold  algorithm. 

CMPDEC  was  used  in  place  of  CLDDEC  in  the 
proces.sing  of  the  "ten-minute"  data  sub-sets.  Fig.  4-2 
and  Tables  4. 1  and  4.3  are  still  valid  if  this  substitution  is 
recognized. 

4.4  An  Alternative  Product:  Calibrated  Radiance 

An  alternative  procedure  involves  the  modification  of 
the  Fig.  4-2  .sequence  to  generate  calibrated  radiance 
imagery,  either  in  addition  to  or  in  lieu  of  the  cloud/no 
cloud  imagery.  This  sequence  modification  is  not  trivial 
since  the  calibrated  radiance  imagery  is  not  normally 
held  intact  in  image  memory,  but  for  maximum  compu¬ 
tational  throughput  is  used  piecemeal  and  not  retained 
once  the  ratio  computations  arc  complete.  An  estimated 
time  estimate  for  implementing  this  alternative  can  also 
be  deduced  from  the  estimates  in  Fig.  4-2. 


Modifying  the  processing  sequence  shown  in  Fig.4-2 
to  produce  output  tapes  of  calibrated  radiance  imagery 
has  not  yet  been  implemented,  but  does  not  represent  a 
major  programming  effort  provided  one  makes  a  clean 
interrupt  at  the  point  of  insertion  within  the  software 
program  TAPRAT.  There  will  be  however,  changes  to 
the  logical  flow  within  the  program  which  may  result  in 
substantial  increases  in  processing  times.  In  any  case,  Utc 
creation  of  a  calibrated  radiance  library  would 
procedurally  involve  the  re-running  of  modified  steps  4, 
5,  &  6  shown  in  Fig.  4-2.  These  modifications  would  be 
of  the  form  shown  in  Fig.  4-4,  with  programs  ABSRAD 
and  SELECT  representing  the  necessary  new  code. 

Since  the  absolute  radiance  values  which  arc  desired 
as  output  reside  temporarily  within  the  existing  sequence, 
and  their  accumulation  should  be  reasonably  simple,  die 
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machine  lime  necessary  lo  execute  llie  ABSRAD  task 
should  represent  a  relatively  short  interval.  A  reasonable 
first  approximation  for  the  ABSRAD  execution  lime  is 
its  being  equivalent  to  or  less  than  the  CLDDEC  runs.  To 
the  degree  with  which  this  approximation  is  true,  then  the 
incremental  elTori  lo  execute  the  sequence  shown  in  Fig. 
4-4  can  be  estimated. 


Assuming  a  clean  insenion,  it  is  rca.sonable  lo  assume 
approximately  3  man-mo  to  code,  test  and  debug  the 
ABSR.\D  &  SELECT  software.  Housekeeping  and 
Oversight  (11/0)  analyst  hours  required  for  processing 
should  not  exceed  the  TAPRAT  rate,  nor  should  the  1/A 
time.  Thus  an  e.slimaie  of  around  4  analyst-hours  per 
radiance  tape,  about  half  the  TAPRAT  rale,  seems  rea¬ 
sonable. 

Rerunning  the  nominal  100  raw  data  tapes  through  the 
Sfl  ECT/ABSRAD  sequence  rcpre.sents  an  incremental 
increase  of  600  machine  hours  and  400  analyst-hours 
(approximately  l.Osystem-mo. and 2. 3 analyst-mo.).  An 
:i. tutted  incremental  co.st  for  implementing  this  alter- 
,  ..JVC  proecssing  sequence  is  readils  determinable. 

-s.O  THE  AinOMATED  CLOUD  DECISION 
A1  CORITHMS 

this  section  discusses  the  methods  used  lo  convert  the 
W.SI  images  lo  quantitative  decisions  regarding  cloud 
cover.  Whereas  it  is  relatively  straightforward  to  acquire 
WSI  images  which  are  pleasing  and  readily  intcrprciablc 
(()  the  human  user,  the  task  of  acquiring  imagery  which 
may  be  automatically  assessed  through  image  processing 
techniques  is  more  demanding.  Requiring  that  llicse 
image  processing  techniques  be  both  fast  and  autono¬ 


mous  is  also  an  important  consideration.  This  section 
discus.ses  first  the  acquisition  and  processing  of  the  data 
to  yield  the  ratio  images,  and  then  the  cloud  algorithms 
used  to  yield  the  cloud  decision  images. 

5.1  Raw  Data  Characteristics 

The  CID  (Charge  Injection  Device)  sensor  has  a 
number  of  charactensiics  which  arc  advantageous  in 
acquiring  full  sky  imagery'.  Most  imporiruitly  is  the  fixed 
gain,  allowing  the  measured  signal  to  be  directly  related 
toUicsky  radiance  fora  given  filler  and  aperture  setting. 
The  sensors  also  have  very  low  blooming,  so  that  even  if 
a  portion  of  the  scene  is  significantly  offscalc  bright, 
other  portions  of  the  scene  are  not  corrupted  radiomeiri- 
cally.  The  noise  is  also  reasonably  low,  typically  2 
counts  on  a  scale  of  0  to  255,  in  a  reasonably  new  camera. 

The  sensors  arc  some  what  limited  in  sensitivity  range. 
We  did  not  wish  lo  use  signals  below  20  counts,  where 
the  least  count  change  of  1  becomes  a  significant  fraction 
of  the  signal.  With  this  limitation,  the  overall  sensitivity 
range  is  somewhat  over  1  log  (a  factor  of  10).  Even  with 
the  sun  filtered  by  a  4  log  neuU'al  density  filter  in  the  sun 
occulior,  the  radiance  range  of  the  upper  hemisphere  is 
largcrthan  1  log.  Forihisrca.son, two pairsoffilters were 
used:  one  set  to  acquire  the  darkest  parts  of  the  sky  on 
.scale,  and  the  other  buffered  by  a  .5  log  neutral  density 
filler  that  puts  the  brighter  portion  on  scale.  (An  offset  of 
.5  log  was  chosen,  rather  Uian  1  log,  to  give  a  region  of 
overlap  lo  be  used  for  accuracy  as.scssmcni.)  Thus  the 
spectral  filter  wheel  held  a  blue  filler  centered  at450nm. 
a  red  at  650  nm,  a  blue  with  .5  log  offset,  and  a  red  with 
.5  log  offset. 

The  sensor  package  al.so  included  a  neutral  density 
filiei  wheel,  lo  allow  the  sensor  to  adjust  lo  the  changing 
light  levels  during  the  day.  Each  minute,  the  sy.stem 
acquired  images  in  the  four  neutral  density  positions,  and 
made  adecision  regarding  which  filicrto  use  forthcnexl 
data  acquisition.  The  llux  control  algorithm  was  de¬ 
signed  to  keep  the  minimum  flux  level  above  40  counts. 
If  there  was  insufneieni  radiant  energy  in  the  open  hole 
(no  attenuation)  position,  the  afX'riurc  was  automatically 
of)cncd,  m  increments  of  I  /2  stop.  The  aperture  w  as  kept 
at  1/16  when  possible,  however,  to  optimize  the  image 
quality  by  maximizing  the  depth  of  focus.  (The  depth  of 
field  with  this  system  is  essentially  infinite,  with  both 
clouds  and  the  bolts  holding  on  the  protective  dome  in 
focus.  However,  the  analogous  depth  of  focus  at  the 
image  plane  was  quite  low ,  resulting  in  defocusing  near 
the  edge  of  the  image  it  the  aperture  w  as  fully  open.) 

Following  data  acquisition,  die  data  tapes  were  pro¬ 
cessed  through  a  quality  control  program ,  which  veri  tied 
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that  data  were  onscalc,  and  that  the  signals  in  the  comers 
of  the  imager’s  active  area  which  arc  outside  the  optical 
image  were  dark,  noise  free,  and  balanced  (i.e.  no  signifi- 
ca:  t  bias  from  left  to  right  orlop  to  bottom).  The  program 
also  verified  that  all  four  spectral  images  were  acquired, 
the  llux  control  sequence  was  completed  successfully, 
and  no  times  were  skipped.  Any  such  problems  were 
noted,  so  that  systems  could  be  repaired  as  necessary  and 
spe  cial  processing  requirements  could  be  addressed. 

5.2  Processing  to  the  Ratio  Stage 

Prior  to  the  automated  ratio  processing,  the  data  were 
checked  in  several  additional  respects.  Images  were 
ev;  luated  to  determine  whether  the  clock  time  was  cor¬ 
rect,  based  on  expected  vs  actual  sun  position,  and  to 
det  .'imine  whether  the  occultor  arm  was  properly  aligned 
with  respect  to  its  readout.  Special  problems  such  as 
debris  on  the  protective  dome  or  use  of  an  incorrect 
occultor  arm  arc  noted  Any  of  these  special  consider¬ 
ations  are  noted  digitally  in  an  input  file.  Abnormalities 
such  as  time  errors  were  corrected;  others  such  as  debris 
on  the  dome  could  not  be  corrected,  but  they  were  noted 
in  the  header  of  the  ratio  tape  and  thence  the  cloud  tape, 
to  make  users  aware  of  the  anticipated  data  quality. 

The  first  step  in  the  ratio  processing  involves  the 
linearity  calibration.  Whereas  the  camera  signal  is  ap¬ 
proximately  linear  with  respect  to  the  input  flux,  this 
calibration  measures  the  response  more  precisely.  The 
linearity  correction  is  quite  fast  using  a  look  up  table 
(firmware  on  the  image  processing  board),  which  re¬ 
places  each  measured  signal  with  the  corrected  signal 
which  is  linearly  related  to  the  flux.  The  calibration  is 
normally  acquired  prior  to  fielding  the  system,  and  the 
look  up  table  is  computed  by  a  program  designed  to 
process  the  linearity  calibration  data. 

Next,  i  f  the  signal  shaveavaluclcss  than  20  counts  out 
of  255,  the  ratio  is  assigned  a  value  of  0.  This  0  value 
represents  regions  external  to  the  optical  image,  and  any 
other  regions  which  may  be  offscalc  dark.  For  tho.se 
picture  elements  (pixels)  which  arc  oascale,  the  ratio  is 
computed  from  the  equation 

Ratio  =  128  *  SPR  *  Radiance  (2)/Radiancc(l)  (5.1) 

The  constants  in  this  equation  will  be  discus.scd  below. 
For  those  pixels  which  arc  offscalc  bright  in  filter  1  or  2 
(blue  or  red),  the  ratio  mu.st  be  computed  using  filters  3 
and  4  (blue  and  red  with  offset),  with  the  following 
equation. 

Ratio  =  128  *  SPR  *  NDR  Radiance  (4)/Radiancc(3) 

(5.2) 


This  creates  a  composite  image,  using  the  2/1  pair  at 
darker  pixel  locations,  and  the  4/3  pair  at  brighter  pixel 
locations.  Those  pixels  which  are  offscalc  bright  in  filter 
3  or  4  are  so  indicated  with  a  specific  byte  value  in  the 
ratio  image  and  later  in  tlic  cloud  dcci.sion  image. 

In  the  above  equations,  the  SPR  is  the  spectral  ratio, 
which  corrects  for  the  small  differences  in  the  passband 
of  the  red  and  blue  filter  (both  nominally  70  nm),  as  well 
as  the  differences  in  the  sensitivity  of  the  CID  at  450  and 
650.  The  NDR  is  the  neutral  density  ratio,  which  corrects 
for  the  small  error  induced  by  the  fact  that  the  .5  log  offset 
filter  is  not  quite  neutral,  and  also  corrects  fordifferenccs 
between  the  3,4  passbands  in  comparison  with  the  1,2 
passbands.  Both  the  SPR  and  the  NDR  ratios  arc  nor¬ 
mally  determined  from  ab.solute  calibrations,  and  veri¬ 
fied  or  adjusted  using  the  field  data.  The  neutral  density 
filters  used  for  flux  control  are  also  slightly  spectrally 
selective;  the  correction  for  this  is  included  in  the  SPR 
and  NDR  terms.  For  this  reason  the  SPR  and  NDR 
factors  are  dependent  on  the  neutral  density  filter  se¬ 
lected  during  flux  control,  and  may  therefore  vary  from 
image  to  image. 

The  above  steps  also  include  a  correction  for  image 
size.  Rather  than  ratioing  the  same  pixels  on  the  active 
array  in  image  space,  the  program  takes  into  account  the 
size  and  placement  of  the  optical  image  on  the  sensor 
array  so  that  it  ratios  the  pixels  which  arc  at  the  same  look 
angle  in  object  space.  For  example,  if  one  imager  is  set 
up  so  that  the  optical  image  is  offset  1 0  pixels  to  the  right 
of  the  centerof  the  active  array,  the  radiance  pixels  10  to 
the  right  of  center  will  be  placed  at  the  center  in  the 
computed  ratio  image.  Thus  all  imagers  arc  corrected  to 
a  standard  size  and  location.  The  exact  placement  of  the 
radiance  images  vary  .somewhat  from  system  to  sy.stcm, 
and  also  from  filter  to  filter  due  to  a  2%  size  difference 
between  the  red  and  blue  images.  Thus  the  image  size 
correction  is  both  system  and  filter  dependent. 

In  the  above  equations,  the  choice  of  the  128  multi¬ 
plier,  and  the  choice  to  use  rcd/bluc  rather  than  bluc/rcd, 
were  driven  by  several  considerations.  The  equation 
should  map  the  theoretically  possible  range  ol  ratios  into 
the  0  -  255  byte  value  range;  a  change  of  1  in  output  byte 
value  .should  correspond  to  roughly  the  same  change  in 
cloud  density  independently  of  cloud  brightness;  mea¬ 
surement  noise  should  not  be  overly  amplified  in  the  ratio 
image;  yet  radiometric  resolution  in  the  radiance  image 
should  not  be  lost  in  converting  to  the  ratio  image. 

The  next  step  in  the  ratio  computation  is  the  overlay  of 
the  mask.  Those  pixels  which  should  be  obscured  by  the 
occultor,  as  well  as  any  which  arc  obscured  by  objects  in 
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ihe  field  of  view  (primarily  at  the  Columbia  site)  ane  set 
10  a  value  of  0.  Thus  they  are  treated  as  “no  data”,  just  as 
the  offscalc  dark  values  are.  Following  this  step,  the  data 
arc  smoothed  with  a  3  x  3  kemal,  to  moderate  the  effects 
of  system  noise. 


The  inverse  equations,  for  obtaining  the  pixel  position 
for  a  given  angular  position,  arc 

X  =  255  +  0.8R  sin4)  (5.6) 

y  =  240  +  R  cos(!)  (5.7) 


Finally,  the  data  are  written  to  tape,  along  with  a 
header  which  contains  identifying  information  regarding 
the  data,  its  quality,  and  its  processing.  The  filter  4  red 
image  is  also  written  to  tape,  to  provide  a  means  to 
evaluate  the  efficacy  of  the  ratio  and  the  cloud  decision 
process. 

In  the  resulting  ratio  image,  the  clear  sky  has  a  lower 
value  than  the  clouds.  Most  of  the  variance  in  the  clear 
sky  radiance  has  been  removed  through  ratioing,  al¬ 
though  the  ratio  is  still  higher  in  the  aureole  and  near  the 
horizon.  Essentially,  the  brightness  variation  caused  by 
changes  in  the  scattering  angle  and  changes  in  the  optical 
depth  of  the  path  (such  as  near  the  horizon)  have  been 
removed.  But  color  variation  caused  by  differences 
between  the  Mic  scattering  properties  at  450  and  at  650 
nm  have  not  been  removed.  The  opaque  clouds  have 
significantly  higher  ratios  than  the  clear  sky  ratios  in 
general.  Grey  storm  clouds,  even  those  darker  than  the 
adjoining  sky,  are  also  characterized  by  high  ratios. 

5.3  Geometric  Calibration  of  the  Images 

Both  the  ratio  images  and  the  cloud  decision  images 
generated  in  the  next  stage  of  processing  have  been 
corrected  to  the  standard  image  size  and  placement.  The 
identification  of  the  look  angle  associated  with  a  given 
pixel  location  is  then  the  same  for  all  sites.  To  a  first 
approximation,  the  fractional  distance  of  a  pixel  to  the 
edge  of  the  image  is  linearly  related  to  the  zenith  angle  of 
the  path  of  si  ght  associated  with  the  pixel. Amoreprecise 
calibration  has  yielded  a  parabolic  fit  of  the  relation 
between  this  fractional  distance  on  the  image  plane  and 
llic  zenith  angle.  The  result  is  the  following  equations, 
which  give  the  zenith  and  azimuth  angles  as  a  function  of 
pixel  position. 


e=585  (l-Vl-.0O1128T) 


(l)  =  tan 


1.25(x-255) 

(y-240) 


where 


(5.3) 


(5.4) 


R=  y[l.25(x-255)f+(y-240)^ 


(5.5) 


where  r  is  now  defined  in  terms  of  0,  by 

R  =  0  (3.032  -  0.00259  0)  (5.8) 

If  the  angular  relation  were  truly  linear,  such  that  r 
were  proportional  to  0,  where  r  is  the  fractional  distance 
to  the  radius,  then  the  solid  angle  seen  by  the  pixel  would 
be  related  to  the  area  of  the  pixel  by  the  following 
proportionality. 

dQ  sin0 

dA  0  (5.9) 

With  the  parabolic  relation,  this  proportionality 
becomes 

—  a  —  ((3.032  -  0.00259  0)(3. 032  -  0.00518  0)]'’ 

^  ®  (5.10) 

An  estimate  of  cloud  cover  may  be  obtained  by 
counting  the  number  of  cloud  vs  no  cloud  pixels.  It  is 
more  accurate  to  weight  each  pixel  by  dQ/dA,  how¬ 
ever  it  can  be  shown  that  the  maximum  error  in  the 
total  sky  cover  which  occurs  from  ignoring  this 
weighting  function  is  approximately  2. 1  %.  For  most 
purposes,  this  weighting  correction  is  not  necessary. 

5.4  Fixed  Threshold  Algorithms 

As  previously  noted,  two  automated  cloud  discrimi¬ 
nation  algorithms  were  developed.  The  fixed  threshold 
algorithm  (CLDDEC)  was  u.sed  in  the  1-min  cloud 
decision  processing,  while  the  more  complicated  com¬ 
posite  fixed-variable  threshold  algorithm  (CMPDEC) 
performed  the  lO-min  decision.  Descriptions  of  thc.sc 
two  techniques  were  presented  by  Koehler,  cl  al  ( 1 99 1; 
at  the  ClDOS-91  conference. 

The  simplest  ratio  cloud  decision  algorithm  assigns  a 
certain  sky  state  to  fixed  ranges  of  ratio  value.  In  our 
application,  we  want  to  discriminate  between  clear  sky, 
thin  cloud  and  opaque  cloud.  Two  values,  the  thin  and 
opaque  thresholds,  need  to  be  defined.  All  pixels  with  a 
red/bluc  ratio  exceeding  the  opaque  threshold  arc  then 
assumed  opaque  cloud.  Those  pixels  with  ratios  between 
the  thresholds  are  assumed  thin  cloud  cover.  Finally, 
those  pixels  with  ratios  below  the  thin  threshold  arc 
assumed  clear. 
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The  fixed  threshold  algorithm  does  a  reasonable  job  in 
defining  opaque  cloud,  but  in  many  circumstances  fails 
to  properly  distinguish  the  true  thin-clear  boundary.  This 
problem  is  illustrated  in  the  example  shown  in  Fig.  5-1. 
Two  contrails  are  evident  in  this  image:  a  denser  one 
oriented  roughly  from  top  to  bottom,  and  a  narrower  one 
angling  from  the  upper  left  to  the  lower  right.  The  region 
just  below  the  solar  occultor  is  brighter  than  other  parts 
of  the  image  due  to  its  proximity  to  the  sun’s  aureole,  but 
should  probably  be  classified  as  clear.  Figure  5-2  shows 
two  different  thin  threshold  slices  through  the  same  ratio 
image.  The  cloud  decision  with  a  thin  threshold  of  1 20 
fails  to  identify  significant  parts  of  the  contrails.  When 
the  thin  threshold  is  brought  down  to  70  to  bring  out  the 
remainder  of  the  contrails,  the  clear  regions  near  the 


Fig,  5-1.  WSI 650  nm  radiance  image  from 
Kirtland  AFB,  8  March  1989 


occultor  and  the  horizon  are  misidentified  as  thin  cloud. 
This  occurs  because  the  clear  sky  ratio  background  is  not 
constant  The  sky  appears  lighter  blue  near  the  horizon 
andnearthesim  than  it  appears  in  the  downsun  direction. 
Thus,  a  thin  cloud  in  the  downsun  part  of  the  sky  may 
actually  have  a  smaller  ratio  than  clear  sky  near  the  sun 
or  horizon.  An  improved  algorithm  would  have  to  take 
this  sky  background  variation  into  account. 

5.5  Variable  Threshold  Algorithm 

Several  similar  contrail  cases  were  studied  in  an  effort 
to  determine  the  relationship  between  thin  clouds  and 
their  background  sky  conditions.  We  found  that  the  ratio 
of  thin  cloud  red/blue  radiance  ratio  to  the  clear  sky 
background  red/blue  radiance  ratio  remains  fairly  con¬ 
stant  along  sections  of  contrails  exhibiting  uniform  opti¬ 
cal  properties.  If  the  clear  sky  background  ratio  were 
known,  a  thin  cloud  discrimination  could  then  be  based 
on  whether  the  observed  ratio  exceeds  the  background 
ratio,  and  the  degree  of  “thinness”  could  be  estimated  by 
the  observed/background  ratio.  The  thin  cloud  decision 
“problem”  then  becomes  one  of  determining  a  reason¬ 
able  estimate  of  the  clear  sky  background  ratio. 

The  rcd/blue  clear  sky  ratio  distribution  is  influenced 
by  many  factors  the  most  important  of  which  are  the  solar 
zenith  angle,  and  the  haze  feamres  of  the  atmospheric 
boundary  layer.  Figure  5-3  shows  the  variation  of  the 
ratio  distribution  as  a  function  of  solar  zenith  angle.  Note 
how  the  clear  sky  background  “whitens”  (the  ratio  in¬ 
creases)  as  the  sun  approaches  the  horizon.  Figure  5-4 
illustrates  the  boundary  layer  haze  influence,  with  the 
pristine  February  arctic  air  mass  yielding  “bluer”  ratios 
than  the  haze  laden  air  mass  from  July. 


Fig.  5-2.  Fixed  threshold  cloud  decisions  with  two  different  thin  cloud  thresholds. 
Lighter  shades  indicate  cloud,  darker  indicate  clear  sky. 
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Our  new  cloud  decision  algorithm  uses  clear  sky 
background  fields  derived  from  clear  day  observations. 
The  haze  effect  is  normalized  by  dividing  the  image  by 
the  reference  red/blue  ratio  from  the  intersections  be¬ 
tween  the  45°  look  zenith  circle,  and  the  locus  of  points 
separated  by  a  45°  scattering  angle  from  the  sun.  The 
effect  of  the  normalization  is  shown  in  Figure  5-5.  In 


Fig.  5-4.  Clear  Sky  Haze  Layer  Dependence: 
Columbia,  MO,  1989 


Fig,  5-5.  Clear  Sky  Normalized  to  Minimize  Haze 


Layer  Dependence:  Columbia,  MO,  1989 


practice,  a  set  of  normalized  clear  sky  ratio  tables  are 
extracted  from  many  clear  days  at  a  specified  site  for  each 
5°  solar  zenith  angle.  The  table  consists  of  normalized 
ratios  at  5°  look  zenith  intervals  from  0  to  75°,  and 
azimuths  relative  to  the  solar  azimuth  at  15°  intervals. 
Tables  for  each  solar  zenith  angle  set  are  then  averaged 
together  to  provide  the  best  normalized  clear  sky  esti¬ 
mate. 

In  the  cloud  decision  algorithm,  the  clear  sky  back¬ 
ground  for  a  particular  image  is  detemiined  as  follows. 
First,  the  solar  zenith  angle  is  determined  for  the  date, 
time  and  site  location  of  the  observation.  A  linear 
interpolation  is  then  made  between  the  two  closest  solar 
zenith  angle  tables  to  estimate  the  normalized  back¬ 
ground  table  for  the  image.  Given  the  zenith  and  relative 
azimuth  of  a  particular  pixel  in  the  image,  its  normalized 
clear  sky  ratio  estimate  is  then  computed  by  bilinear 
interpolation  of  the  table  values.  This  normalized  value 
is  then  multiplied  by  a  reference  ratio  value  determined 
from  an  interactive  procedure  that  examines  the  refer¬ 
ence  value  variation  for  that  day  and  time. 

The  cloud  decision  is  then  made  by  comparing  the 
observed  ratio  to  the  fixed  opaque  threshold,  and  the  clear 
sky  ratio  estimate.  The  pixel  is  then  given  a  numerical 
value  that  defines  whether  it  is  categorized  as  clear,  thin 
cloud,  opaque  cloud,  obscured  from  view,  or  indetermi¬ 
nate.  This  last  category  occurs  when  the  clear  sky  ratio 
estimate  exceeds  the  opaque  cloud  threshold.  The  cloud 
decision  image  from  the  new  algorithm  for  the  contrail 
case  is  shown  in  Figure  5-6.  Note  how  the  contrails  are 
identified  without  misidentifying  the  points  near  the  sun 
and  horizon. 


Fig.  5-6.  Variable  threshold  cloud  decision  image 


In  the  resulting  cloud  decision  image,  those  points 
identified  in  the  ratio  process  as  “no  data”,  offscale 


Fig.  5-3.  Clear  Sky  Solar  Zenith  Angle  Dependence: 
1  l/Feb/89  -  Columbia,  MO 
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bright,  or  indetenninatc,  arc  identified  by  specific  byte 
value.  The  indeterminate  category  represents  cases  where 
the  clear  ratio  estimate  exceeds  the  opaque  threshold  due 
to  Jie  presence  of  thick  haze.  This  is  analogous  to  the 
case  when  the  observer  cannot  visually  distinguish  the 
white  cloud  in  the  haze. 

The  clear,  thin,  and  opaque  cloud  identifications  arc 
each  given  a  specific  range  of  values  in  the  cloud  decision 
imjge,  i.c.  clear  =  1  to  99,  thin  =  100  to  1 39,  and  opaque 
=  140  to  200.  In  the  clear  and  opaciuc  categories,  the 
variance  in  the  final  representation  is  proponional  to  the 
raiiO.  But  in  the  thin  category,  the  variance  is  related  to 
the  perturbation  of  the  ratio  with  respect  to  the  sky 
background  ratio,  and  thus  related  to  the  optical  depth  of 
the  thin  cloud. 

This  approach  was  used  to  satisfy  two  needs.  If  it  is 
desired  to  simply  count  ilic  number  of  pixels  in  the 
various  categories  .such  as  clear  or  opaque,  this  may  be 
done  unambiguously.  Yet  the  variance  within  each 
category  is  retained,  so  that  tlic  image  can  be  false  color 
displayed  with  its  texture  retained,  and  so  that  the  cloud 
image  may  be  reprocessed  to  investigate  the  effect  of 
different  threshold  selections.  The  specific  rules  used  in 
the  clouddecision  image  format  are  given  in  Appendix  B. 

It  should  also  be  noted  that  the  cloud  decision  images 
include  headers,  much  like  tlio.se  used  in  the  ratio  images, 
which  give  information  on  the  data,  its  quality,  and  its 
processing. 


tcnthsoftoial  .sky  covercategory  differences  between  the 
WSl  and  the  human  ob.server  at  Columbia  from  the  14- 
month  sample.  (For  example,  if  the  WSl  detects  50% 
cloud  cover  and  the  ob.server  60%,  this  is  a  -1  category 
difference. lvalues  from  both  the  old,  fixed  threshold, 
and  new,  variable  threshold  algorithms  are  shown.  The 
new  algorithm  shows  better  agreement  with  the  human 
observer.  This  comparison  only  indicates  cases  where 
the  improvement  was  due  to  changing  the  WSl  total  sky 
cover  amount.  In  many  ca.ses,  the  total  sky  cover  did  not 
change  but  the  spatial  distribution  of  that  percentage 
improved. 

6.0  DATA  BASE  APPLICATIONS 

It  is  a  well  accepted  truism  that  as  more  and  more  vital 
elements  of  our  national  defense  become  increasingly 
sophisticated,  they  also  become  increasingly  sensitive  to 
the  influences  imposed  upon  them  by  their  surrounding 
operational  environment.  This  fact  is  nowhere  more 
clearly  illustrated  than  in  the  arena  of  electro-optical 
system  performances  in  turbid,  cloudy  atmospheres  where 
cloud  attenuation  of  the  sensor/target  path  of  sight  can 
result  in  the  complete  and  near  instantaneous  negation  of 
system  performance.  It  is  therefore  only  prudent  that 
intense  and  innovative  research  into  the  development  of 
techniques  for  the  as.sessment  and  prediction  of  appropri¬ 
ate  cloud  population  characteristics  and  cloud  field  dy¬ 
namics,  particularly  at  small  spatial  and  temporal  scales, 
be  actively  pursued. 


5.6  Data  Base  Processing  with  the  New  Algorithm 

The  new  variable  threshold,  directionally-dcpendent 
algorithm  has  been  used  to  process  the  10-minute,  cloud 
decision.  The  results  arc  encouraging.  The  more  reliable 
definition  of  thin  cloud  illustrated  in  the  contrail  example 
is  the  rule,  not  the  exception.  A  measure  of  this  improve¬ 
ment  is  shown  in  Figure  5-7.  This  figure  illustrates  the 
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CLOUD  COVER  CATEGORY  DIFFERENCE  (WSl  OBSERVED) 


Fig.  5-7.  Total  Cloud  Cover: 

Columbia.  MO  9  Feb  89  -  28  Mar  90 


The  WSl  image  oriented  data  base,  by  virtue  of  its 
relatively  high  spatial  and  temporal  resolution,  is  particu¬ 
larly  appropriate  for  applications  associated  with  the 
acquisition  and  tracking  of  earth  orbiting  devices.  The 
statistical  characterization  of  local  cloud  effects  upon 
earth  to  space  optical  communication  links  is,  of  course, 
the  reason  the  automatic  multi-station  data  collection 
network  was  originally  established.  The  generalized 
output  product  from  tliis  network,  the  temporally  varying 
pixel  by  pixel  definitit)n  of  cloud  thickness  throughout 
the  observable  sky  dome  is  directly  applicable  to  all  cloud 
free  line  of  sight  (CFLOS)  and  cloud  free  arc  (CFARC) 
determinations  either  as  real  time  operational  support 
products,  or  as  statistically  predictive  model  inputs. 
Ground  ba.sed  laser  operations  which  are  particularly 
impacted  by  cloud  field  dynamics  can  employ  the  appro¬ 
priate  variations  of  CFLOS/CFARC  persistence  and  re¬ 
currence  statistics  directly  into  mis.sion  planning  .sce¬ 
narios  with  increasingly  short  respon.se  times. 

Several  examples  illu.strating  the  varied  use  of  the 
WSl  image  data  are  presented  in  tlie  following  summa¬ 
ries.  They  are  of  course  user  specific  applications  con- 


ducted  undcrlhe  general  scope  of  this  particular  contrac¬ 
tual  effort,  but  readily  typify  the  types  of  tasks  for  which 
the  data  base  is  deemed  appropriate. 

6.1  Concepts  for  Anti-Satellite  (AS AT)  Engage¬ 

ment  Scenarios 

With  regard  to  determining  the  impact  o f  cloud  e f feels 
upon  definition  of  a  specific  ASAT  systems  engagement 
space,  it  is  clear  that  there  are  several  contributing  ele¬ 
ments.  The  characteristics  of  the  satellites  orbit,  the 
characteristics  of  the  intervening  atmosphere  specifi¬ 
cally  including  the  nature  of  the  prevailing  cloud  field, 
and  the  nature  of  the  ground  based  laser  itself  must  all  be 
considered.  Whereas  it  is  appropriate  to  be  fully  cogni¬ 
zant  of  the  methodology  intended  for  the  specification  of 
these  various  characteristics,  we  feel  that  the  optimum 
analytic  scheme  would  employ  as  complete  a  decoupling 
between  the  major  sub-elements  as  is  practical.  This  is 
particularly  true  with  regard  to  the  atmospherics.  To  the 
degree  that  this  decouoling  of  the  cloud  field  statistics 
and  related  characterizations  from  the  other  perfomiance 
characterizations  within  the  general  ASAT  system  per¬ 
formance  model  can  be  achieved,  the  more  readily  a  truly 
definitive  set  of  CFARC  statistics  can  be  generated.  In 
the  long  term  view,  we  feel  that  atmospheric  decoupling 
is  an  extremely  important  issue  and  should  be  addressed 
early  and  creatively  in  defining  the  analytic  approach. 

For  the  application  at  hand,  i.c.,  meteorological  sup¬ 
port  products  for  application  within  a  mlatively  generic 
ASAT  Scenario,  we  will  assume  meteorological 
decoupling,  and  the  production  of  statistical  products 
!fom  the  WSl  data  ba.se  for  White  Sands  Missile  Range. 
The  statistics  have  been  extracted  for  a  sample  two  week 
interval  for  procedural  illustration  only.  i.e.  for  mo.st 
realistic  applications  substantially  longer  sample  inter¬ 
vals,  sea.sonally  weighted  would  be  more  appropriate, 
riicsc  concepts  have  been  previously  discu.ssed  in  Opti¬ 
cal  Systems  Group,  Technical  Note  No.  219,  T.  L 
Koehler  and  J.  E.  Shields  (1990). 

For  this  demonstration,  an  arbitrary  arc  trajectory, 
apresented  by  pixel  row  324  on  the  image,  was  chosen 
10  [provide  a  simulated  satellite  track,  and  a  variety  of 
sample  statistics  related  to  this  arc  were  generated.  The 
gens' ral  t|iiestion  to  be  an.swcred  was:  What  is  the 
probability  of  finding  acloud  free  arc  length  sufficient  to 
allow  a  successful  deposition  ol  energy  uixm  the  salcl- 
lite In  this  example,  ".successful "  was  to  be  defined  by 
the  users  system  performance  model  in  a  manner 
decoupled  from  the  meteorological  detemiinations. 

I  he  first  lour  figures  illustrate  the  statistical  data 
extracted  Irom  the  two  week  sample  of  WSl  imagery. 


The  last  three  illustrate  an  assumed  procedure  for  apply¬ 
ing  these  statistical  data  to  a  generic  ASAT  engagement. 

6.1.1  Graphics  &  Illustrations 

Unconditional  Prohabihty  Plots 

Figure  6.1-1  shows  the  frequency  distribution  for  start 
and  end  point  of  the  longest  cloud  free  arc  occurring 
along  the  selected  arc  f  row  324)  within  each  image.  Only 
those  points  which  are  neither  completely  clear  nor 
completely  cloudy  are  included  on  the  plot.  (The  fre¬ 
quencies  of  dear  and  overcast  cases  are  printed  on  the 
plot.)  The  shape  of  this  curv'c  appears  reasonable;  even 
for  this  small  sample  used  in  its  derivation. 

The  .second  plot  on  this  page.  Fig.  6.1-2,  gives  the 
frequency  distribution  for  arc  length.  (Again,  cleai  and 
overcast  are  not  included.)  This  is  followed  by  Fig. 
6.1-3,  a  plot  of  the  cumulative  frequency  for  arc  length, 
bolhforthe  total  sample,  and  for  A.M.  (the  3  hours  before 
local  apparent  noon)  and  P.M..  (the  three  hours  after), 
illustrating  a  simple  temporal  sorting  option. 


Unconditional  Prohuhilin  Interpretation 

From  the  combined  cumulative  frequency  plot,  Fig. 
6.1-3,  one  can  see  that  the  20')  value  corresponds,  for 
this  sample,  to  an  arc  length  of  4.3  degrees.  Although 
these  numbers  should  not  tx'  trusted  at  this  point  as 
definitive  statistics,  they  can  be  u.scd  for  illustrative 
purposes.  In  this  sample,  il  one  wishes  an  80'^  success 
rate,  one  must  be  able  to  input  sufficient  energy  in  45 
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Table  6.1.1 
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Fig.  6.1-3 

degrees  of  transit  time.  Note  that  the  afternoon  was 
considerably  more  cloudy  than  the  morning.  In  the 
morning  one  has  over  60  degrees  available,  but  in  the 
afternoon  only  30  degrees  arc  available,  assuming  a 
desired  80%  success  rate. 

Conditional  Probability  Plot 

The  start  point  along  a  selected  arc  can  also  be  quite 
critical  in  determinations  of  whether  it  is  possible  to  input 
satisfactory  energy  levels  onto  the  target  vehicle.  There¬ 
fore  we  also  computed  some  conditional  probabilities. 
One  can  first  compute  the  probability  that  the  start  point 
is  a  given  position.  Then  for  a  given  start  point,  one  can 
compute  the  conditional  probability  that  the  arc  length  is 
a  given  amount.  The  application  of  this  approach  is 
illustrated  in  the  remaining  plots  and  tables. 

The  first  plot.  Fig.  6. 1  -4  and  Table  6.1.1,  illustrate  the 
arc  length  conditional  probability,  based  on  this  limited 
data  sample,  given  the  start  pixel  position  shown.  For 
illustration  only ,  the  arc  has  been  broken  into  9  segments. 


CONDITIONAL  CLOUD  FREE  ARC  PROBABILITY 


Fig.  6.1-4 


(These  computations  would  normally  be  done  at  full 
pixel-by-pixel  resolution.)  Thus,  forexamplc,  .7%  of  the 
cases  had  a  start  point  within  category  1  and  arc  length 
within  category  1  (roughly  0-12  degrees). 

Application  of  the  Conditional  Probability 

Once  similar  conditional  probabilities  are  computed 
for  sufficiently  long  periods,  these  data  could  be  used  in 
the  following  way.  The  ASAT  model  could  be  used  to 
generate  the  amount  of  energy  which  can  be  input  for  a 
given  start  point  and  arc  length.  This  would  be  computed 
for  all  combinations  of  start  point  and  arc  length.  This  has 
been  done  for  a  simulated  energy  profile,  in  the  figure 
labeled  "Simulated  Energy  Matrix",  Fig.  6. 1  -5  and  Table 
6.1-2. 


Fig.  6.U 


Ifonc  knows  the  energy  required  formission  success, 
the  energy  matrix  is  thrcsholded  at  this  known  required 
energy.  The  resulting  thrcsholded  matrix  is  0  if  a  given 
start  point  and  arc  length  do  not  yield  sufficient  energy. 
A  sample  thresholded  energy  matrix  has  been  plotted  as 
Fig.  6.1-6  and  listed  in  Table  6.1.3  for  a  successful 
engagement  energy  of  7000  relative  units.  This  binary 
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matrix  is  then  multiplied  by  the  conditional  cloud  free  arc 
length  matrix  shown  earlier,  to  determine  how  often 
success  is  achieved.  The  sum  of  this  final  matrix  yields 
percent  success  rate. 

Table  6.1.2 


Simulaied  Eoergy  Mitnx 
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Table  6.1.3 
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The  final  plot,  Fig.  6.1-7,  compares  the  resulting 
success  rate  for  the  unconditional  and  conditional  prob¬ 
abilities.  The  unconditional  probability  gives  a  some¬ 
what  more  conservative  answer.  'Fhat  is,  for  a  given 
required  energy,  the  unconditional  probability  will  yield 
a  slightly  lower  estimate  of  the  success  rate.  For  a 
required  energy  of  7000,  the  conditional  probabihty 
predicts  a  2%  higher  success  rate. 
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For  operational  purposes,  ihc  important  question  is: 
how  much  will  we  overestimate  the  energy  output  re- 
quirementsof  file  GBL  i  f  we  use  unconditional  probabili¬ 
ties.  Aaswering  this  directly  requires  multiple  runs  with 
different  input  energy  profiles.  However  we  can  ap¬ 
proach  the  problem  indirectly  by  asking  a  slightly  differ¬ 
ent  question;  what  is  the  highest  cumulative  incident 
energy  that  could  be  achieved  for  a  given  fixed  input 
energy.  At  the  80%  point,  the  unconditional  probability 
implies  an  energy  of  only  47(X)  units  would  be  available 
on  target;  the  conditional  probability  implies  that  7200 
units  of  energy  would  be  available,  for  the  same  input 
profile.  That  is,  the  unconditional  probability  estimate 
implies  at  least  7200  energy  units  would  be  available 
H()%  of  the  time.  Thus  in  temis  of  energy,  the  conditional 
probability  gives  a  much  more  favorable  answer. 

It  should  tv  noted  that  the  stair-step  nature  of  the 
unconditional  percent  success  vs.  required  energy  plot  is 
an  artifact  of  the  cc'arse  arc  length  resolution  used  for  the 
Illustration  of  concept  in  tliis  sample.  This  artifact  would 
disappear  when  the  conputations  are  made  at  full  resolu¬ 
tion 

Summary 

The  cumulative  frequency  of  arc  length  based  on  the 
unconditional  probabililv  of  arc  length,  yields  a  very 
convenient  result.  One  can  determine  immediately  the 
arc  length  in  w  hich  sufficient  energy  must  be  input.  This 
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is  a  very  convenient  format,  as  it  uncouples  the  cloud 
statistics  hx)m  the  ASAT  model;  these  distributions  may 
be  run  for  as  many  cases  as  desired,  and  input  as 
distributions  to  the  ASAT  model. 

The  conditional  probabilities  also  decoupled  the  cloud 
impact  computations  from  the  ASAT  computations. 
That  is,  the  cloud  statistics  may  be  determined  indepen¬ 
dently  from  many  cases,  and  input  as  a  distribution  to  the 
ASAT  model.  They  also,  in  our  opinion,  yield  a  some¬ 
what  more  accurate  result  than  the  unconditional  prob¬ 
abilities.  They  are  not  as  convenient  to  interpret  how¬ 
ever,  because  determination  of  the  final  percent  success 
rate  requires  use  of  an  energy  matrix  from  the  ASAT 
model. 

6J!  Preliminary  Cloud  Free  Arc  (CFARC) 
Statistics  from  WSI  Imagery 

Introduction 

The  study  described  in  Section  6. 1  was  designed  to 
illustrate  the  procedure  for  combining  CFARC  statistics 
with  laser  energy  profiles  to  compute  the  probability  of 
success  in  an  ASAT  scenario.  The  WSI  data  for  that 
study  was  extracted  from  an  arbitrary  easi-west  row  for 
two  one-week  samples  that  exhibited  reasonable  cloud 
cover  variations.  In  contrast,  the  following  section 
describes  a  brief  study  using  CFARC  infonmation  from 
WSI  data  collected  at  the  C-Slation  site  during  the  two 
extreme  months  of  the  year  in  terms  of  the  climatological 
cloud  cover  at  WSMR,  along  a  north-south  arc  that  better 
simulates  a  realistic  satellite  track.  The  preliminary  one- 
minute  cloud  decision  imagery  generated  with  the  fixed 
threshold  algorithm  available  at  the  time  of  this  study 
were  employed. 

6.2.1  White  Sands  Cloud  Cover  Climatology 

Considering  our  desire  to  limit  the  resources  allocated 
to  this  task,  and  the  preliminary  nature  ol  the  cloud/no 
cloud  data  available  at  the  time  of  this  study,  we  decided 
to  focus  on  two  monthly  samples  that  illustrate  the 
extremes  in  the  cloud  cover  distributions  at  WSMR. 
Since  the  degree  of  cloud  cover  along  an  arc  in  the  sky  is 
strongly  dependent  on  the  total  cloud  cover,  the  selection 
of  which  monthly  periods  to  examine  was  based  on  the 
climatological  sky  cover  data  from  Stallion  Site  at  WSMR 
illustrated  in  Fig.  6.2- 1 .  The  strong summerdiumal  cycle 
with  clearer  mornings  and  cloudier  afternoons  is  clearly 
evident.  The  extreme  periods  (black  circles)  occur  in 
consecutive  months,  with  June  AM  being  the  clearest  on 
average,  and  June  PM  being  cloudiest  on  average.  The 
(+)  symbols  show  the  corresponding  C  Station  WSI 
average  cloud  covers  from  1989.  The  WSI  one  year 
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Fig.  6.2-1 


average  matches  the  Stallion  1 3-year  site  climatology  in 
the  July  PM  cloudy  period.  In  contrast,  the  June  AM 
averages  differ  significantly. 

Fig.6.2-2  illusU’ates  the  cloud  cover  distributions  from 
WSI  and  climatology  for  the  periods  of  interest.  The 
completely  clear  and  completely  cloudy  categories  are  in 
close  agreement  for  the  June  AM  results.  However,  the 
sky  cover  amounts  in  the  3-9  tenths  categories  are  smaller 
in  the  WSI  results  compared  to  climatology,  appearing 
instead  in  the  1  tenth  category,  leading  to  the  lower  WSI 
average  cloud  cover.  Similarly,  the  WSI  results  from 
July  PM  ate  smaller  in  the  same  3-9  tenths  categories  than 
those  from  climatology,  but  both  the  1  tenth  and  10  tenths 
cases  have  increased,  resulting  in  no  difference  in  aver¬ 
age  total  cloud  cover.  The  differences  noted  can  arise 
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from  several  sources,  including  natural  year-to-year  vari¬ 
ability,  different  sites  (C-Station  is  approximately  100 
miles  south  of  Stallion  Site)  and  surrounding  terrain,  and 
known  biases  in  the  human  observations  vs.  WSI  (see, 
for  example,  McGuffie  and  Hendcrson-Scllcrs,1989). 

6.2.2  CFARC  Data  Reduction 

The  basic  data  employed  in  the  study  is  a  subset  of  the 
preliminary  one-minute  cloud/no  cloud  data  release  made 
earlier  to  Ihe  Analytic  Sciences  Corp.  (TASC).  Soft¬ 
ware  was  written  to  extract  a  337  pixel  cloud/no  cloud 
sample  from  image  column  255  starting  at  row  72  at  5 
minute  intervals.  This  column  runs  from  south  to  north 
through  zenith,  and  the  sample  is  bounded  by  the  60° 
zenith  circle  (Fig.  6.2-3).For  each  subset  extracted,  the 
maximum  arc  free  of  either  opaque  (white)  or  thin 
(shaded)  cloud  pixels  was  determined,  and  the  corre¬ 
sponding  starting  pixel  and  arc  length  values  saved.  The 
maximum  CFARC  in  Fig.  6.2-3  is  highlighted  in  dashed 
white.  All  available  data  from  the  C-Station  site,  27  days 
for  June  1989  and  28  days  from  July  1989,  were  pro¬ 
cessed  in  this  manner. 


Fig.  6.2-3 


The  results  that  follow  are  based  only  on  the  morning 
data  from  June  and  the  afternoon  data  from  July.  Any 
period  with  an  arc  containing  occulter  mask  pixels  was 
ignored. 

6.2.3  Resultant  CFARC  Distributions 

General  CFARC  Distributions 

The  first  step  in  the  data  analysis  was  to  classify  the 
maximum  CFARC  starling  and  ending  pixel  pocitions 
and  the  resulting  arc  lengths  into  24  categories,  each 
category  representing  an  arc  of  5°  zenith  angle.  The  June 
AM  results  are  shown  in  Fig.  6.2-4.  All  three  di.stribu- 
tions  do  not  include  the  completely  cloudy  arc  (4.6%)  or 


completely  clear  arc  (75.4%)  cases.  As  expected,  the 
start  pixel  distribution  has  a  maximum  at  the  south  edge, 
and  the  end  pixel  distribution  has  a  maximum  at  the  north 
edge.  The  start  pixel  distribution  also  has  an  unusual 
secondary  maximum  near  the  center  of  the  arc,  due 
primarily  to  spurious  cloud  identifications  made  near  the 
sun,  but  outside  the  occulter  mask  in  the  20  to  30  minute 
periods  just  before  and  after  the  mask  intersects  the  arc  of 
interest  A  similar  peak  occurs  in  the  center  of  the  arc 
length  distribution,  implying  that  these  are  probably 
completely  clear  arc  cases  that  have  been  misidentified 
as  partial  arc  cases.  Thus,  the  completely  clear  case 
frequency  would  be  increased  by  almost  2%.  Note  also 
the  the  end  pixel  maximum  exceeds  the  start  pixel  maxi¬ 
mum  by  over  6%,  indicating  an  azimuthal  preference  for 
clouds  to  the  south  compared  to  the  north  during  the  June 
mornings. 


As  might  be  expected,  the  July  PM  results  (Fig.  6.2- 
5)  have  greater  frequency  values  than  those  from  June 
AM,  due  to  fewer  completely  clear  cases.  Maximum 
CFARCs  starting  at  the  southern  edge  now  exceed  those 
ending  on  the  northern  edge,  suggesting  a  reversal  in  the 
cloud  preference  compared  to  June  AM.  The  strong  peak 
in  the  first  arc  length  category  (0°  -  5°)  is  probably  caused 
by  anom  alous  clear  pixels  appearing  in  completely  cloudy 
CFARC  situations. 

Cumulative  CFARC  Length  Distributions 

Users  of  CFARC  information  often  ask  for  the  prob¬ 
ability  of  having  an  arc  length  of  at  least  some  specified 
value,  or  alternatively,  the  arc  length  that  is  met  or 
exceeded  some  specified  percentage  of  the  time.  Fig. 
6.2-6  provides  just  such  information,  with  the  vertical 
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axis  depicting  the  probability  of  exceeding  the  maximum 
CFARC  specified  along  the  horizontal  axis.  For  ex¬ 
ample,  the  maximum  CFARC  along  the  specified  track  is 
grcatcrthan40°90%  ofthetimcinthe  June  AM  sample, but 
only  48%  in  the  July  PM  sample.  Note  also  that  the 
completely  clear  arc  frequency  appears  at  the  right  side 
of  the  diagram,  and  the  completely  cloudy  frequency  can 
be  obtained  by  subtracting  the  value  on  the  extreme  left 
from  100%.  Another  interesting  feature  of  these  results 
is  that  fewer  than  80%  of  the  July  PM  cases  have  clear 
arcs  of  any  length. 

The  following  four  figures  provide  additional  insight 
needed  for  interpreting  the  results  in  Fig.  6.2-6.  Fig.  6.2- 
7  combines  the  WSI  cloud  determination  image  from 
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Fig.  6.2-6 


Fig.  6.2-7 


Fig.  6.2-3  with  the  June  AM  curve  from  Fig.  6.2-6.  The 
maximum  CFARC  for  the  image  shown  in  Fig.  6.2-7  was 
74°,  as  indicated  by  the  dotted  vertical  line  in  the  inset 
graph.  This  length  was  exceeded  roughly  80%  of  the 
lime  in  the  June  AM  sample.  In  other  words,  80%  of  the 
cases  had  maximum  CFARCs  clearer  than  this  case. 

In  the  worst  case  scenario  (July  PM),  the  80%  level 
must  include  completely  cloudy  arcs,  as  shown  illus¬ 
trated  in  Figs  6.2-8  and  6.2-9.  While  it  is  easy  to  see  that 
a  selected  arc  will  be  completely  cloudy  in  an  overcast 
situation  (Fig.  6.2-8),  one  must  also  realize  that  com¬ 
pletely  cloudy  arcs  also  arise  in  cases  with  partial  cloud 
cover  (Fig.  6.2-9).  In  further  contrast  with  the  best  case 
example  in  Fig.  6.2-7,  Fig.  6.2- 1 0  illustrates  a  case  with 
a  35°  maximum  CFARC.  Only  50%  of  the  July  PM  cases 
exceeded  this  value. 


Fig.  6.2-8 
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Fig.  6.2-9 


Conditional  Probability  Distributions 

In  many  laser  scenarios,  the  position  of  the  maximum 
CFARC  can  be  almost  as  important  as  its  length.  Tables 
6.2. 1  and  6.2.2  give  the  percentage  occurrences  of  maxi- 


Tabu  6.2.1  •  Conditlooal  nuilmum  ('FARC  probability  (%]  given  arc  sUrt 
poaiUoo  for  the  June  A.M  sample  touling  IJ65  cases. 
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Table  6.2  J  ■  Conditional  maiimum  ('FARC  probabQily  (^1  given  arc  start 
povilion  for  the  July  FM  sample  totaling  1572  cacses. 
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mum  CFARC  length,  given  that  the  CFARC  begins  at 
certain  positions  along  the  arc.  Tlie  data  have  been 
organized  in  twelve  10°  categories.  For  example,  0.4% 
of  the  1572  July  PM  maximum  CFARCs  started  in  the 
southern  30°  -  40°  zenith  angle  interval  and  extended  into 
the  30°  -  40°  arc  length  interval,  as  did  the  case  shown  in 
Fig.  6.2-10. 


Fig.  6.2-10 


Fig.  6.2-1 1  gives  a  graphical  illustration  of  the  condi¬ 
tional  probability  tables.  The  maxima  in  stan  pixel 
categories  6  and  7  (the  region  1 0°  cn  cither  side  of  zenith) 
in  the  June  AM  panel  correspond  to  the  secondary 
maximum  noted  previously  in  the  start  and  end  pixel 
di.stributions.  ITic  July  PM  maximum  in  arc  length 
category  1  (0°  -  5°)  suggests  that  the  peak  in  this  category 
noted  carherin  the  bottom  panel  of  Fig.  6.2-5  is  localized 
in  one  part  of  the  arc,  and  may  be  produced  by 
nonuni formitics  in  the  imaging  chip.  Except  for  the 
problems  just  noted,  the  conditional  probability  values 
arc  quite  reasonable,  and  arc  consistent  in  structure  to 
those  from  other  periods  (not  shown). 

6.2.4  Concluding  Remarks 

Potential  u.sers  of  the  WSI  data  base  should  llnd  the 
results  from  this  study  encouraging.  Previous  compari¬ 
sons  reported  in  Shields,  ct  al,  ( 1 990)  showed  reasonable 
agreement  between  overall  sky  coveresiimates  made  by 
the  WSI  and  the  human  ob.server.  Tlie  CFARC  results 
presented  here  imply  tliat  useful  information  regarding 
finer  scale  features,  such  as  single  arcs,  can  also  be 
extracted  from  the  WSI  data. 

Tliese  data  and  illustrations  have  been  previously 
discu.sscd  in  Optical  Systems  Group  Technical  Note  No. 
222.  T.  L.  Koehler,  1990, 
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Conditional  Cloud  Free  Arc  Probability 
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Fig,  6.2-11 


6.3  Comparisons  of  WSI  Cloud  Determinations 
with  Trained  Observers 


63.1  Total  Cloud  Cover  Comparison 

Much  of  the  WSI  data  from  the  Columbia  site  was 
proces.sed  using  the  preliminary  (fixed  threshold)  algo¬ 
rithm,  in  order  to  give  a  preliminary  assessment  of  the 
accuracy  of  these  techniques.  Fig.  6.3-1  shows  the  cloud 
cover  distribution  from  7  months  of  WSI  data,  compared 
with  the  values  of  total  sky  cover  reported  on  the  National 
Weather  Service  Form  lO's.  The  WSI  values  arc  from  the 
one  minute  image  at  the  reported  time  of  the  weather 
observation.  This  plot  is  for  the  six  hours  surrounding 
local  apparent  noon.  The  comparison  between  WSI  and 
observer  is  in  general  quite  good.  The  cloud  algorithm 
identifies  some  clear  cases  as  1/10  cover,  but  in  all  other 
cloud  categories  the  match  in  frequency  of  observance  is 
excellent. 
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Fig.  6.3-1 


Another  indication  of  data  quality  is  a  direct  casc-by 
case  compari.son  between  WSI  and  observer.  For  Fig. 
6.3-2,  the  difference  between  WSI  and  observer  has  been 
computed  for  each  case.  Thatis.a  WSI  value  of  7/10  and 
observer  value  of  5/10  would  be  a  difference  of  2 


Whencvcrautomatic  machines  arc  developed  to  mimic 
.some  spcci  lie  element  of  human  pcrfomiancc,  there  is  an 
immediate  and  necessary  demand  to  demonstrate  cred¬ 
ibility.  The  human's  ability  to  perceive  and  interpret 
visible  spectrum  imagery  in  real  time  is,  and  will  surely 
remain  for  some  while,  unchallenged.  However,  there 
arc  a  myriad  of  human-like  image  acquisition  and  inter¬ 
pretation  tasks  that  profit  immensely  from  the  ability  to 
store  and  later  retrieve  images  in  a  reliable  manner.  One 
of  ihe.se  tasks  is  obviously  WSI  cloud  determination  and 
documentation.  A  preliminary'  as.scssmcnt  of  the  WSI 
systems  ability  to  mimic  the  human  in  this  particular  task 
is  summarized  below.  This  information  has  previously 
appeared  in  Shields,  ct  a1.  1990. 
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categories.  The  distribution  of  category  differences  is 
shown  in  Fig.  6.3-2.  The  majority  of  the  cases  show  a 
category  difference  of  0.  The  average  difference  is  less 
than  half  a  category,  i.e.  much  less  than  1/10  cloud  cover. 

6.3.2  Specific  Temporal  Dynamics  Case  Study 

Although  the  observer  to  WSI  comparison  is  gener¬ 
ally  good,  the  WSI  data  can  show  much  more  variance 
than  the  observer  data,  simply  due  to  the  limited  temporal 
frequency  of  the  observer  values.  One  particularly 
dynamic  day,  14  April  at  Columbia,  Mo.,  is  illustrated  in 
Figs.  6.3-3  and  6.3-4.  In  Fig.  6.3-3,  which  shows  the 
cloud  cover  determinations  for  both  the  WSI  and  the 
observer,  the  observer  values  are  consistent  with  the 
WSI,  but  they  certainly  do  not  show  the  true  variability, 
due  to  their  limited  temporal  frequency.  For  example,  at 
1350  and  1450,  both  WSI  and  observer  show  approxi¬ 
mately  20%  cloud  cover,  but  during  the  intervening  hour 
the  cloud  cover  increased  to  nearly  80%.  The  cloud 
images  at  the  middle  and  end  of  this  hour  are  shown  in 
Figs.  6.3-5  and  6.3-6.  The  fastest  rate  of  change  during 
this  hour  occurred  at  between  1400  and  1413,  when  the 
cloud  cover  changed  from  20  to  74%. 


Fig.  6.3-3 


Fig.  6.3-4 


Fig.  6.3-6 


One  can  also  see,  in  Figs.  6.3-5  and  6.3-6,  a  horizontal 
line  part  way  down  the  image  in  the  Northern  sky.  This 
line  represents  an  arc  traveling  from  horizon  to  horizon, 
rising  to  a  45  degree  zenith  angle.  The  maximum  cloud 
free  arc  length  (CFARC)  along  this  arc  has  been  plotted 
in  Fig.  6.3-4.  Comparing  Figs.  6.3-3  and  6.3-4,  one  can 
see  that  in  general  there  is  a  tendency  for  short  CFARCs 
tooccur  with  high  cloud  cover.  As  expected  however,  the 
maximum  CFARC  is  not  always  well  related  to  the  cloud 
cover  at  a  given  point  in  time.  For  example,  at  1450,  the 
cloud  free  arc  length  was  quite  short,  in  spite  of  the  low 
cloud  cover.  Examination  of  the  cloud  image  for  this 
time  (Fig.  6.3-6)  shows  that  this  occurred  because  the 
cloud  band  happened  to  lie  over  the  selected  arc. 

6.4  Applications  to  Model  Development  and 
Evaluation 

The  emergence  of  the  WSI  cloud  data  base  has 
prompted  a  variety  of  cloud  model  evaluations  and 
reviews,  many  of  which  have  been  discussed  within  the 
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user  community  via  the  annual  Proceedings  of  the  Goud 
Impacts  on  DoD  Operations  and  Systems  (CIDOS)  con¬ 
ferences.  Two  short  studies  conducted  at  MPL  related  to 
th(  modeling  of  CFLOS/CFARC  statistical  properties 
have  been  presented  previously  as  part  of  this  currently 
reported  research  effort.  Portions  of  these  earlier  reports 
have  been  abstracted,  and  arc  included  in  the  following 
paragraphs. 

6.^  .1  Evaluation  of  Stochastic  Models  for  Estimat¬ 
ing  the  Persistence  Probability  of  Cloud  Free  Lines- 
of-Sight 

The  report  titled  above,  as  is  Sectior.  6.4. 1 ,  was  issued 
as  GL-TR-89-0275,  Scientific  Report  No.  1,  30  Sep  89, 
authored  by  Wayne  S.  Hcring.  Its  abstract  and  introduc¬ 
tory  paragraphs  are  reproduced  below. 

Abstract 

Stochastic  models,  based  upon  the  Omstein-Uhlenbeck 
(O-  U)  class  of  the  simple  Markov  process  have  been  used 
effectively  to  estimate  the  conditional  probability  of  a 
variety  of  weatherevcnts(Gringoncn,  1972).  This  study 
is  directed  toward  extension  of  the  analytical  form  of  the 
0-U  Markov  model  to  yield  estimates  of  the  joint  occur¬ 
rence  probability  and  duration  of  sky  cover  and  cloud 
free  lines  of  sight  (CFLOS).  As  a  first  step,  model 
estimates  are  made  of  the  recurrence  probability  of 
CFLOS  for  specific  categories  of  sky  cover  in  tenths  as 
a  function  of  time.  Then  model  calculations  are  made  for 
the  persistence  probability  of  CFLOS  and  sky  cover  as  a 
function  of  sky-covcrcategory  and  lime  using  an  analyti¬ 
cal  representation  of  the  mathematic  solution  of  persis¬ 
tence  probability  for  the  0-U  Markov  process  given  by 
Kielson  and  Ross  (1975).  The  model  calculations  of 
CLFOS  recurrence  and  persistence  are  compared  with 
the  Columbia,  MO,  data  base  (Lund,  1973).  Again  using 
an  analytical  form  of  the  0-U  Markov  model,  calcula¬ 
tions  arc  made  of  the  joint  occurrence  probability  of  sky- 
cover  at  two  or  more  sites  and  compared  with  the  obser¬ 
vations  from  the  central  United  Slates.  Finally,  analytic 
techniques  arc  presented  for  calculating  probability  esti¬ 
mates  of  the  duration  of  cloud-free  or  cloudy  lines  of 
sight  from  one  or  multiple  ground  sites  to  points  in  space. 
Trial  calculations  based  upon  climatic  summaries  of  sky- 
cover  are  made  for  a  selected  group  of  sites  in  south 
western  United  States. 

Introduction 

Minute  by  minute  determinations  of  sky-cover  and 
cloi'd-frce-linc-of-sight  (CFLOS)  arc  being  made  in  an 
experimental  field  program  initiated  in  1988  by  the 
University  of  California,  San  Diego.  Solid-state,  whole- 


sky  imagery  (WSl)  systems  have  been  installed  in  a 
network  configuration  consisting  of  6  field  sites  in  the 
western  United  States.  The  primary  goal  is  to  obtain 
detailed  observations  of  the  joint  occurrence  frcquenc) 
of  CFTOS  in  time  and  .space  that  arc  required  to  evaluate 
and  extend  sky-cover  models.  The  stochastic  sky-cover 
models  are  to  be  used  for  estimating  the  impact  of  clouds 
on  ground-based  systems  that  depend  upon  unobscured 
paths  of  sight  to  satellites  in  space. 

A  comprehensive  statistical  model  has  been  devel¬ 
oped  by  Boehm,  ct  al,  (1986)  that  spxtcifically  provides 
estimates  of  the  duration  of  cloud  free  lines  of  sight  from 
multiple  ground  sites  to  orbiting  and  geostationary  satel¬ 
lites.  The  innovative  method  establishes  the  climatic 
probabilities  through  repetitive  simulations  of  sky-cover 
distributions  with  the  multidimensional  Boehm  Sawtooth 
Wave  Model.  The  simulation  model  and  its  many  com¬ 
ponent  approximations  arc  the  prime  candidates  for  test 
and  evaluation  with  the  new  WSl  data  base. 

Prominent  among  alternate  approaches  to  the  problem 
arc  models  based  on  the  Omstein-Uhlenbeck  (0-U)  class 
of  the  simple  Markov  process.  This  approach  has  been 
applied  successfully  both  in  purely  analytical  form  and 
by  Monte  Carlo  simulation  of  probability  distributioas  to 
estimate  the  joint  occurrence  and  duration  of  a  variety  of 
weather  events  (Gringorlen  1966,  1967,  1968,  1972). 
This  study  is  directed  toward  extension  of  the  analytical 
form  of  the  0-U  Markov  model  to  yield  estimates  of  the 
joint  occurrence  and  persistence  probability  of  cloud  free 
lines  of  sight  in  time  and  space. 

The  ultimate  objective  of  the  modeling  process  is  to 
determine  the  joint  climatological  probability  of  the 
duration  of  cloud  free  lines  of  sight  to  a  point  in  space 
from  one  or  more  preselected  ground  sites,  given  the 
climatic  summaries  of  sky-cover  at  each  site.  The  solu¬ 
tion  to  a  number  of  intermediate  modeling  problems  is 
required  to  achieve  the  desired  objective.  In  Section  2, 
the  application  of  existing  models  to  estimate  cloud  free 
line  of  sight  (CFLOS)  probability  as  a  function  of  sky- 
cover  and  the  zenith  angle  of  the  path  of  sight  is  de¬ 
scribed.  The  basic  0-U  Markov  modeling  concepts  are 
reviewed  in  Section  3.  Estimates  of  the  single  station 
duration  of  CFLOS  as  a  function  of  sky-cover  using  the 
Keilson-Ross  procedure  are  discussed  in  Section  4.  Es¬ 
timates  of  the  joint  conditional  probability  of  sky-cover 
at  multiple  sites  is  discu.ssed  in  Section  5.  Procedures  for 
estimating  the  frequency  of  downtime  intervals  caused 
by  cloud  obscured  lines  of  sight  persisting  concurrently 
at  designated  network  sites  arc  described  in  Section  6. 
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6.4.2  Probability  Estimates  of  Cloud-Obscured 
Line-of-Sight 

The  report  titled  above,  as  is  Section  6.4.2,  was  issued 
as  Scientific  Report  No.  2,  30  May  90,  authored  by 
Wayne  S.  Hering.  Its  abstract  and  introductory  para¬ 
graphs  arc  reproduced  below. 

Abstract 

The  purpose  of  this  study  is  to  explore  the  basic 
relationships  between  the  climatic  frequency  distribu¬ 
tion  of  sky  cover  and  the  frequency  and  duration  of 
Cloud y-Lincs-of-Sight  (CLOS)  from  one  ormore  ground 
sites  to  points  in  space.  Through  a  series  of  model 
calculations,  estimates  are  made  of  the  incremental 
changes  in  CLOS  persistence  and  recurrence  that  corre¬ 
spond  to  changes  in  the  climatic  frequency  of  overcast, 
broken,  scattered  and  clear  sky  conditions. 

The  determination  of  CLOS  recurrence  and  persis¬ 
tence  probability  were  made  using  the  bivariate  normal 
distribution  (see  Gringorten,  1972).  The  analytic  sto¬ 
chastic  model  is  based  upon  the  Omstein-Uhlenbeck 
class  of  the  simple  Markov  process.  A  discussion  of  the 
procedures  for  modeling  CLOS  behavior  is  given  in  a 
previous  report  (Hering,  1989). 

Model  results  are  presented  for  CLOS  as  a  function  of 
sky  cover  and  path  of  sight,  CLOS  recurrence  and  pcrsi.s- 
icnce  probability  relative  to  fixed  categories  of  sky  cover, 
and  both  single  and  multiple  site  CLOS  recurrence  and 
persistence  probabilities. 

Introduction 

The  purpose  of  this  study  is  to  explore  the  basic 
relationships  between  the  climatic  frequency  distribu¬ 
tion  of  sky  cover  and  the  frequency  and  duration  of 
Qoudy-Lines-of-Sight  (CLOS)  from  one  ormore  ground 
sites  to  points  in  space.  The  determination  of  CLOS 
probability  is  of  particular  importance  for  the  assessment 
of  expected  downtime  of  ground-ba.sed  laser,  optical  data 
link  and  orbital  surveillance  systems.  Through  a  scries  of 
rruHlel  calcuidtions,  •’stimaics  are  made  of  the  incremen¬ 
tal  changes  in  CLOS  persistence  and  iccurrcncc  that 
correspond  to  changes  in  the  climatic  frequency  of  over¬ 
cast,  broken,  scattered  and  clear  sky  conditions, 

The  model  calculations  of  CLOS  behavior  arc  of  tlie 
following  type: 

a.  Given  the  climatic  frequency  of  sky  coverat  a  selected 
ground  site,  what  is  the  likelihood  that  a  path  of  sight  with 
/criith  angle  0 

1 )  will  be  obscured  by  clouds? 


2)  will  be  obscured  by  clouds  at  a  particular  time 
and  then  be  obscured  by  clouds  again  after  time 
interval  t? 

3)  will  be  obscured  by  clouds  continuously  over 
time  interval  t? 

b.  Given  the  climatic  frequency  distribution  of  sky  cover 
at  a  group  of  selected  ground  sites,  what  is  the  joint 
occurrence  frequency  and  joint  persistence  probability  of 
CLOS  events  such  as  those  listed  in  ( 1 ),  (2)  and  (3)  above 
as  a  function  of  site  separation  distance  d? 

Modeling  Procedures 

The  determination  of  CLOS  recurrence  and  persis¬ 
tence  probability  in  this  study  were  made  using  the 
bivariate  normal  distribution  (see  Gringorten,  1972). 
The  analytic  stochastic  model  is  based  upon  the  Omstein- 
Uhlenbeck  class  of  the  simple  Markov  process.  A 
discussion  of  the  procedures  for  modeling  CLOS  behav¬ 
ior  is  given  in  a  previous  report  (Hering,  1989). 

Markov  model  determinations  of  event  probabilities 
are  made  through  initial  conversion  of  climatic  fre¬ 
quency  distributions  of  sky  cover  into  corresponding 
values  of  Equivalent  Normal  Deviate  (END).  The  trans¬ 
formation  for  any  cumulative  frequency  distribution  can 
be  made  using  standard  statistical  tables  or  by 
computationally  fast  algorithms. 

The  fundamental  modeling  assumption  is  that  the 
autocorrelation  coefficients  for  sky  cover  decay  expo¬ 
nentially  with  time  or  distance.  Also,  the  correlation 
coefficient  for  CLOS  under  fixed  sky  cover  conditions  is 
assumed  to  decay  exponentially  with  time.  The  expres- 
•sion  for  the  temporal  autocorrelation  coefficient  is  sim¬ 
ply 

p,  =  exp  (-t/x) ,  (1) 

where  t  is  elapsed  time  and  i  is  the  preselected  relaxation 
time.  The  corresponding  expression  for  the  spatial 
autocorrelation  coefficient  is 

pjj  =  exp  f-d/D) ,  (2) 

where  d  is  the  site  separation  distance  and  D  is  the  input 
value  for  ihc  relaxation  distance.  These  correlation 
functions  arc  used  in  detenniiiing  mcurrcncc  probabili¬ 
ties. 

In  this  study,  nominal  input  values  were  used  for  sky 
cover  relaxation  time  (16  hours),  sky  cover  relaxation 
distance  (5(X)  mi)  and  CLOS  relaxation  time  (30  min). 
The  cffccLs  of  uncertainties  in  these  values  will  be  as¬ 
sessed  in  association  with  some  of  the  results  as  tabulated 
below. 
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Tlic  modeling  procedures  (llcring,  1989)  arc  simple 
and  easy  lo  apply.  The  translormaiion  lo  ihc  equivalent 
normal  deviate  values  arc  exact  for  any  climatic  cumu¬ 
lative  frequency  distribution.  Once  made,  the  model 
estimates  of  joint  event  (recurrence)  probability  are 
given  by  a  simple  regression  equation  for  the  bivariate 
normal  distribution  as  follows; 

Yt  =  Pt  Yo  +  ni  (3) 

where  rit  is  the  END  corresponding  to  the  conditional 
probability  Pr  (y  ^  yt  1  yo).  and  yo  and  yt  are  END 
equivalents  of  the  cumulative  frequency  of  the  initial  and 
final  events,  respectively.  Thus,  the  expression  yields 
the  conditional  probability  of  a  weathcrevent  following 
a  prescribed  initial  condition. 

The  value  of  y^  in  Eq.  3  is  defined  uniquely  for 
continuous  variables  by  the  cumulative  frequency  distri¬ 
bution  at  the  initial  time.  For  variables  expressed  only  in 
categories,  such  as  sky  cover  amount,  it  is  important  to 
subdivide  the  category  probability  range  into  subsets 
with  smaller  but  equal  ranges.  The  calculations  of 
conditional  probability  should  be  carried  out  in  turn 
using  each  of  the  subset  midpoints  as  y^  and  the  results 
averaged  to  yield  the  composite  result  for  the  sky-cover 
category.  Experience  shows  that  division  into  6  subsets 
is  sufficient  for  good  results. 


the  bivariate  normal  distribution  is  given  by  Keilson  and 
Ross  (1975).  Since  the  computer  routine  for  the  formal 
solution  is  lengthy,  an  alternative  analytical  solution  was 
developed  for  the  CLOS  determinations  which  provides 
reliable  approximation  of  the  formal  solution  over  the 
desired  range  of  output.  The  Keilson-Ross  solution  for 
the  case  where  the  climatic  cumulative  frequency  of  the 
weather  event  is  50  percent  is  simply 

Fo  =  (l/7i)  sin  ‘  (exp -t/x)  (4) 

where  Fq  is  the  unconditional  probability  that  y  ^  yo 
throughout  the  time  interval  t. 

The  .solutions  for  y^  0  in  this  study  were  approxi¬ 
mated  by 

fr(yo)  =  f(Fo)  +  yo(i+oi3e-'’-^‘^") 


-2  <  V  <2 
0<e-«‘>-"^  <3 


where  fr(yo)  is  the  END  corresponding  to  the  probability 
that  y  ^  yo  throughout  the  time  interval,  f  ( Fo)is  the  END 
of  Fo  in  Eq.  4,  and  yo  is  the  unconditional  probability  of 
the  weather  event.  In  contrast  with  Eq.  2  for  recurrence 
probability,  the  persistence  expression  (Eq.  5)  assumes 
that  the  climatic  frequency  of  the  event  remains  the  same 
throughout  the  interval. 

The  sequence  of  four  steps  for  the  determination  of 
CLOS  recurrence  and  persistence  for  both  single  and 
multiple  sites  is  as  follows: 

a.  calculation  of  CLOS  probability  as  a  function  of 
zenith  angle  q  and  particular  category  of  sky  cover 
amount, 

b.  calculation  of  CLOS  recurrence  and  persistence 
probability  as  a  function  of  time  interval  for  specific 
individual  categories  of  sky  cover, 

c.  calculation  of  sky-cover  amount  recurrence  and 
persistence  probability  as  a  function  of  time  interval,  site 
separation  distance,  and  the  climatic  frequency  distribu¬ 
tion  of  sky  cover. 

d.  given  a,  b  and  c,  calculation  of  the  recurrence  and 
persistence  probability  of  CLOS  at  a  single  site  and  for 
simultaneous  recurrence  and  persistence  of  CLOS  at 
multiple  sites  as  a  functionof  sky-cover  frequency  distri¬ 
bution,  time  interval  and  site  separation. 

7.0  RECOMMENDATIONS  FOR  FURTHER 
RESEARCH 

The  WSI  cloud  data  base,  summarized  in  Tables  3.1 
and  4.3,  provides  an  extensive  selection  of  digitized 
cloud  imagery  suitable  for  a  broad  variety  of  modeling 
and  ca.se  study  applications.  It  has  been  compared  with 
data  derived  from  classical  observations  by  meteorologi¬ 
cally  trained  personnel,  and  found  to  produce  equivalent 
statistical  characterizations.  Thus,  it  is  well  linked  to 
other  long  term  historical  archives,  for  which  it  provides 
a  useful  high  resolution  supplement.  Experience  with  the 
data  to  date  clearly  illustrates  the  technical  viability  of 
fully  automatic  digital  .systems  as  reliable  vehicles  for  the 
acquisition  and  analysis  of  whole  sky  cloud  cover  char¬ 
acterizations. 

Model  validation  studies,  particularly  those  related  to 
relatively  short  term  phenomena,  such  as  CFLOS  or 
CLOS  (Cloudy  Line  of  Sight)  recurrence  and/or  persis¬ 
tence  statistics,  arc  an  appropriate  arena  for  WSI  class 
data.  As  abstracted  from  Hcring  (1990); 

"The  data  base  generated  by  the  Whole  Sky  Imagery 
(WSI)  network  will  enable  validation  of  the  many  factors 


involved  in  ihe  CLOS  modeling  process.  The  trial 
calculations  made  in  this  study  provide  a  first  approxi¬ 
mation  of  CLOS  l>chavior  relative  to  the  climatic  fre¬ 
quency  distribution  of  sky  cover.  However,  certain 
aspects  of  the  modeling  process  remained  to  be  analyzed 
and  refined  on  the  basis  of  more  extensive  CLOS  data. 

The  relationships  between  CLOS  frequency  and  the 
zenith  angle  of  the  path  of  sight  will  be  examined  in  detail 
using  the  new  WSI  data  base.  Important  refinements 
include  the  definition  of  CLOS  dependence  on  zenith 
angle  as  a  function  of  cloud  type.  In  addition,  it  is 
important  to  take  a  close  look  at  the  relative  importance 
of  the  azimuthal  dependence  of  the  climatic  frequency  of 
CLOS  as  revealed  by  the  data  from  individual  sites  in  the 
WSI  network. 

The  temporal  and  spatial  variations  in  sky  cover  have 
been  analyzed  rather  extensively  using  the  conventional 
climatic  data  base.  Thus,  tlie  unccnaintics  in  the  model 
input  values  of  sky  cover  relaxation  time  and  distance 
and  the  resultant  impact  on  CLOS  probability  estimates 
can  be  estimated.  However,  the  more  detailed  WSI 
network  data  arc  necessary  for  comprehensive  analysis 
of  cloud-clement  relaxation  time  as  defined  in  paragraph 
4  0,  tuid  its  variations  with  season  and  site  location.  Yet 
another  factor  requiring  study  is  the  potential  variation  in 
cloud-elcmcnt  relaxation  time  witli  the  zenith  angle  of 
the  path  of  sight.  As  shown  by  time  lapse  sequences  of 
ground-based  cloud  imagery,  apparent  cloud  motion 
overhead  is  more  rapid  than  near  the  horizon  " 

Whereas  the  existing  WSI  imagery'  has  proven  its 
utility  in  a  variety  of  daytime  applications,  its  usefulness 
could  be  substantially  broadened  if  it  were  to  contain 
nighttime  cloud  data  in  addition  to  its  current  daytime 
inventory.  A  new  day/night  fisheye  system  has  been 
devclofxjd  at  the  Marine  Physical  Laboratory  which  will 
enable  the  addition  of  these  low  flux  level  images.  The 
data  obtained  by  these  24  hour/day  imagers  is  directly 
applicable  to  extending  the  characicri/alion  of  cloud 
ficid  statistics  throughout  the  24  hour  diurnal  cycle. 
Sl\idies  of  cloud  field  dynamics  dunng  ihc  dawn  and 
twilight  transitions  would  be  particularly  informative, 
and  readily  achievable  via  an  expanded  24  hr/day  image 
archive.  Two  prototype  systems,  developed  under  coop¬ 
erative  tri-service  support  arc  currently  operational  at 
MPL  under  interactive  computer  control,  with  a  techni¬ 
cal  as.se.ssment  of  full  automatic  operation  under  review. 
Applying  systems  of  this  class  to  the  extension  of  the 
existing  cloud  data  base  provides  exciting  research  op¬ 
portunities  iliat  should  promptly  be  pursued. 


Continued  research  into  the  characrerization  of  thin 
cirrus  clouds,  both  daytime  and  nighttime,  is  an  appropri¬ 
ate  application  for  the  WSI  data  base.  With  the  imple¬ 
mentation  of  the  new  cloud  discrimination  algorithm 
desenbed  in  Section  5  of  this  report,  the  separation  of 
clear  sky  background  and  iliin  cloud  definition  is  sub¬ 
stantially  enhanced.  Thus,  the  opponunity  for  thin  cloud 
optical  depth  approximations  is  likewise  improved,  par¬ 
ticularly  w  hen  conducted  as  a  joint  exercise  w'ith  active 
probing  systems.  Nighttime  optical  depth  estimates 
using  background  starfields  as  the  baseline  reference 
level  is  an  attractive  and  achievable  experimental  en¬ 
deavor. 
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Summary 


Responding  to  a  well  recognized  need  by  many  in  both  the  modeling  and 
operational  communities  for  an  improved  capability  in  the  collection  and 
assessment  of  whole  sky  cloud  characteristics,  a  new  generation  of  video 
based  imaging  systems  has  been  developed  and  fielded  by  the  Marine 
Physical  Laboratory.  One  of  these  systems,  the  Whole  Sky  Imager,  has 
been  deployed  at  several  widely  separated  portions  of  the  United  States,  and 
has  gathered  several  million  images  appropriate  for  determining  cloud  cover 
at  very  high  spatial  and  temporal  resolution.  Cloud  cover  estimates  derived 
from  a  7-month  sample  of  these  cloud  images  shows  very  good  agreement 
with  observed  sky  cover  amounts.  The  capabilities  of  the  Whole  Sky 
Imager  and  an  overview  of  the  current  status  and  quality  of  the  WSI  data 
base  are  discussed  in  Technical  Note  No.  226. 

An  inventory  of  the  complete  WSI  Image  Data  base,  segregated  by  site, 
date  and  archival  tape  number,  is  included  herein.  The  data  in  this  catalog 
has  been  reviewed  in  accordance  with  our  established  QC  procedures,  and 
has  been  purged  of  all  sets  deemed  inappropriate  for  further  application. 
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2.7  National  Weather  Facility,  Columbia,  Missouri  (COL) .  62 
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F'g-  No.  Title  Page  No. 

1  Whole  Sky  Imager  Sites  .  34 

2  WSI  Data  Status  &  Processing  Summary .  35 
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APPENDFX  A 


1.0  INTRODUCTION 

This  catalog  is  a  complete  listing  of  all  available  data  deemed  of  adequate  quality  for 
follow-on  processing  and  analysis  from  each  of  the  seven  WSI  sites.  Figure  1  shows  the 
location  of  each  data  site.  Processing  quality  extends  from  excellent  (new  camera)  to 
degraded  but  useable  (old  camera).  As  a  camera  ages  the  imagery  becomes  grainy  in 
appearance,  but  this  does  not  normally  interfere  with  Ratio  or  Cloud/No  Cloud  processing 
or  determination.  Figure  2  shows  the  status  of  each  site  and  the  level  of  processing 
completed  for  that  site.  A  summary  of  data  quality  for  any  site  and  time  are  available  upon 
request. 

Throughout  the  catalog  there  are  gaps  between  tape  sequence  numbers.  The  data 
contained  on  these  unlisted  tapes  is  not  deemed  suitable  for  general  distribution.  In  .some 
cases  the  data  can  be  retrieved  but  would  require  additional  special  processing.  Further 
evaluation  of  this  data  is  available  upon  request. 

Each  site  was  given  a  3  character  identifier  as  noted  on  Figure  1.  This  identifier  is 
written  to  the  header  of  every  image  on  the  tape. 

When  a  station  tape  was  received  from  the  field  it  was  labelled  and  assigned  the  next 
available  sequence  number  starting  with  001  for  each  station.  The  tape  label  consists  of  a  3 
character  station  identifier  and  a  3  digit  sequence  number,  e.g.  "WSCl  12"  identifies  a  tape 
from  the  White  Sands  C-Station  site  (WSC),  which  was  the  1 12th  tape  received. 

U.sers  should  be  aware  that  sequential  data  days  need  not  necessarily  occur  on  tapes  of 
consecutive  sequence  numbers  due  to  occasional  disruptions  at  the  host  facilities. 
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WHOLE  SKY  IMAGER  DATA  SITES 


•  WSC  •  KAA  •  CL4  •  BAR 

•  WSH  •  MAG  •  COL 


Figure  1 
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APPENDIX  A 


2.1  White  Sands  Missile  Range/C-Station,  New  Mexico 


TAPE  • 

DATE 

START 

STOP 

TAPE  1 

DATE 

START 

STOP 

YEAR  = 

JIM 

W3C052 

5-11/AUG 

13:20 

01:19 

WSC034 

12-18/APR 

05:50 

17:48 

WSC053 

12/AUG 

14:40 

15:20 

WSC035 

19-25/APR 

05:50 

17:48 

HSC054 

21/AUG 

13:10 

09:11 

22/AUG 

13:10 

15:09 

VSC036 

26-30/APR 

05:50 

17:48 

1-2/APR 

05:50 

17:48 

WSC055 

23/AUG 

14:50 

01:19 

24-25/AUG 

13:10 

01:09 

WSC037 

3-5/MAY 

05:50 

17:48 

26/AUG 

13:10 

07:11 

6/MAY 

05:50 

07:48 

WSC056 

26/AUG 

16:00 

01:11 

irsc038 

7-13/MAY 

05:50 

17:48 

WSC057 

27/AUG 

17:40 

01:09 

WSC039 

14-17/MAY 

05:50 

17:48 

28/AUG 

13:10 

14:32 

WSC040 

18-23/MAY 

05:50 

17:48 

WSC071 

3/NOV 

17:40 

00:49 

24/MAY 

05:50 

09:00 

4-9/NOV 

12:50 

00:49 

*J^LLJ?AYS  FROH  THIS  PQIMT.PM 

ARE  ON  ZULU -TIME  YEAR  = 


WSC043 

10/JUN 

14:30 

01:09 

11-15/JUN 

13:10 

01:08 

WSC084 

27/JAN 

14:40 

01:19 

16/JUN 

13:10 

15:26 

20-31/JAN 

13:20 

01:19 

1-2/FEB 

13:20 

01:19 

WSC044 

16/JUN 

18:10 

01:08 

17-20/JUN 

13:10 

01:00 

WSC085 

3/FEB 

15:40 

01:19 

4-6/FEB 

13:20 

01:19 

WSC045 

29/JUN 

16:50 

01:11 

7/FEB 

16:50 

17:49 

30/JUN 

13:10 

01:08 

1/JUL 

23:09 

00:58 

WSC086 

8 /FEB 

23:30 

01:19 

9-10/FEB 

13:20 

01:19 

WSC046 

8/JUL 

14:11 

01:09 

9/JUL 

13:10 

17:46 

WSC087 

11/FEB 

14:20 

01:19 

12-17/FEB 

13:20 

01:19 

W3C047 

11/JUL 

19:20 

19:56 

WSCO80 

10/FEB 

13:20 

01:19 

WSC048 

14-18/JUL 

13:20 

01:18 

19/JUL 

13:20 

17:56 

WSC089 

19/FEB 

23:30 

01:19 

20/FEB 

13:20 

01:19 

WSC049 

20/JUL 

20:10 

01:19 

21/FEB 

13:20 

23:45 

21-25/JUL 

13:20 

01:19 

26/JUL 

13:20 

15:20 

WSC090 

22-25/FEB 

13:20 

01:19 

26/FEB 

13:20 

16:47 

W3C050 

27/JUL 

14:20 

01:19 

28/JUL 

13:20 

15:19 

WSC091 

26/FEB 

23:20 

01:19 

27/FEB 

13:20 

20:54 

W3C051 

2/AUG 

16:30 

01:19 

3/AUG 

13:20 

01:19 

WSC092 

2/MAR 

22:20 

01:19 

4/AUG 

13:20 

20:48 

3-5/MAR 

13:20 

01:19 
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Ai  l’F'NDIX  A 


WliJli'  Sand':  Missil*'  Kan>’,o /('-S  t  a  I  i  on  ,  N('w  Moxico 


TAPE  1 

date _ 

START 

STOP 

lAEfiJL 

date _ 

START 

STOP 

WSC093 

6-9/MAR 

13:20 

01:19 

WSC112 

27/JUN 

16:20 

01:09 

10/MAR 

13:20 

23:11 

28-30/JUN 

13:10 

01:09 

1-3/JUL 

13:10 

01:09 

W3C094 

10/MAR 

23:40 

01:19 

11-16/MAR 

13:20 

01:19 

WSC113 

6/JUL 

22:20 

01:09 

7-12/JUL 

13:10 

01:09 

WSC095 

17/MAR 

13:20 

01:19 

18/MAR 

13:20 

23:17 

WSC114 

13/JUL 

14:40 

23:08 

1fSC09S 

21 /MAR 

13:50 

01:19 

HSC115 

16-21/JUL 

13:10 

01:09 

22-27/MAR 

13:00 

01:19 

WSC117 

23-29/JUL 

13:10 

01:09 

K3C098 

5-11/APR 

13:10 

01:09 

WSC118 

30-31/JUL 

13:10 

01:09 

WSC099 

19/APR 

14:10 

01:09 

1-3/AUG 

13:10 

01:09 

20-25/APR 

13:10 

01:09 

4 /AUG 

13:10 

19:51 

WSClOO 

26-30/APR 

13:10 

01:09 

WSC119 

5-9/AUG 

13:10 

01:09 

1-2/MAY 

13:10 

01:09 

10/AUG 

13:10 

19:55 

WSClOl 

3-9/MAY 

13:10 

01:09 

WSC120 

10/AUG 

21:40 

01:09 

11-16/AUG 

13:10 

01:09 

WSC102 

10-13/MAY 

13:10 

01:09 

l-l/MAY 

14:00 

01:09 

WSC121 

17/AUG 

14:20 

01:09 

15-16/MAY 

13:10 

01:09 

18/AUG 

13:10 

01:09 

WSC103 

17-23/MAY 

13:10 

01:09 

WSC122 

19-22/AUG 

13:10 

01:09 

23/AUG 

13:10 

23:03 

WSC104 

24/MAY 

22:30 

01:09 

25-30/MAY 

13:10 

01:09 

WSC124 

24/AUG 

23:02 

01:09 

25-30/AUG 

13:10 

01:09 

WSC105 

31/MAY 

13:10 

01:09 

1-3/JUN 

13:10 

01:09 

WSC125 

31 /AUG 

13:10 

01:09 

4/JU(4 

13:10 

15:10 

1-4/SEP 

13:10 

01:09 

5-6/SEP 

13:00 

00:59 

WSC106 

4/JUN 

17:40 

01:09 

5-10/JUN 

13:10 

01:09 

WSC126 

7-13/SEP 

13:00 

00:59 

W3C107 

11-12/JUN 

13:10 

01:09 

WSC127 

14-20/SEP 

13:00 

00:59 

13/JUN 

13:10 

20:39 

WSC128 

21-27/SEP 

13:00 

00:59 

W3C108 

14-19/JUN 

13:10 

01:09 

WSC129 

28-30/SEP 

13:00 

00:59 

WSC109 

20/JUN 

18:00 

01:09 

1/OCT 

13:00 

00:59 

21/JUN 

13:10 

15:36 

2-4/OCT 

12:50 

00:49 

WSCllO 

21/JUN 

18:30 

01:09 

WSC130 

5-11/OCT 

12:50 

00:49 

22-23/JUN 

13:10 

01:09 

WSC131 

12-15/OCT 

12:50 

00:49 

USClll 

4-5/JUL 

13:10 

01:09 

6/JUL 

13:10 

21:47 

WSC132 

16-22/OCT 

12:50 

00:49 
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APPENDIX  A 


White  Sands  Missile  Range/C-Statlon,  New  Mexico 


TAPE  1 

DATE 

START 

STOP 

TAPE  1 

DATE 

START 

STOP 

W3C133 

23-29/OCT 

12:50 

00:49 

WSC151 

6-12/FEB 

13:20 

01:19 

irsci34 

30-31/OCT 

12:50 

00:49 

WSC152 

13/FEB 

15:10 

01:19 

1-5/NOV 

12:50 

00:49 

14-19/FEB 

13:20 

01:19 

WSC135 

6-12/NOV 

12:50 

00:49 

WSC153 

20/FEB 

15:20 

01:19 

21-26/FEB 

13:20 

01:19 

1fSCl36 

13-19/NOV 

12:50 

00:49 

WSC154 

27/FEB 

13:20 

01:19 

WSC137 

20/NOV 

12:50 

00:49 

28/FEB 

13:20 

20:06 

21/NOV 

12:50 

19:09 

HSC155 

28/FEB 

20:50 

01:19 

WSC138 

21/NOV 

19:40 

00:49 

1/MAR 

13:20 

01:19 

22-27/NOV 

12:50 

00:49 

WSC156 

2-6/HAR 

13:20 

01:19 

WSC139 

28/NOV 

14:40 

00:49 

7/MAR 

13:20 

18:18 

29/NOV 

12:50 

00:49 

30/NOV 

13:00 

00:59 

KSC157 

7/MAR 

19:00 

01:19 

1-4/DEC 

13:00 

00:59 

8-13/MAR 

13:20 

01:19 

WSC140 

5-10/DEC 

13:00 

00:59 

WSC158 

14/MAR 

14:40 

01:19 

15/MAR 

13:10 

01:09 

1fSC141 

12-18/DEC 

13:00 

00:59 

16/MAR 

13:10 

20:37 

WSC142 

19/DEC 

13:00 

18:20 

WSC159 

16/MAR 

22:00 

01:09 

17-22/MAR 

13:10 

01:09 

lfSCl43 

19/DEC 

18:50 

00:59 

20-22/DEC 

13:00 

00:59 

WSC160 

23-29/MAR 

13:10 

01:09 

23-25/DEC 

13:10 

01:09 

WSC161 

30/MAR 

15:30 

01:09 

WSC144 

26-31/DEC 

13:10 

01:09 

31/MAR 

13:10 

01:09 

1-5/APR 

13:10 

01:09 

■Ilia 

WSC162 

6-9/APR 

13:10 

01:09 

10/APR 

13:10 

17:58 

WSC144 

1/JAN 

13:10 

01:09 

WSC163 

10/APR 

18:20 

01:09 

11-16/APR 

13:10 

01:09 

WSC145 

2-8/JAH 

13:10 

01:09 

WSC1G4 

17/APR 

13:10 

01:09 

trSC146 

9-15/JAN 

13:10 

01:09 

18-23/APR 

13:00 

00:59 

WSC147 

16-22/JAN 

13:20 

01:19 

WSC165 

24-25/APR 

13:00 

00:59 

26/APR 

13:00 

00:00 

WSC148 

23/JAN 

13:20 

17:41 

WSC167 

1/MAY 

18:30 

00:59 

IfSCl49 

23/JAN 

19:00 

01:19 

2-7/MAY 

13:00 

00:59 

24-29/JAN 

13:20 

01:19 

WSC168 

8-14/MAY 

13:00 

00:59 

WSC150 

30/JAN 

20:30 

01:19 

31/JAN 

13:20 

01:19 

WSC169 

15/MAY 

13:20 

00:59 

1-5/FEB 

13:20 

01:19 

16-21/MAY 

13:00 

00:59 
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Ai’i’f'Nnrx  A 


WtilLe 

TAPE  i 

Sands  Miss 

DATE 

i  1  f  RnnRf“/c;-SL.i  t  i  ('ll , 

5TMT  STOP 

New  Mexico 

TAPE  1 

DATE 

START 

STOP 

WSC170 

22-24/HAy 

13:00 

00:59 

WSC188 

29-30/NOV 

13:00 

00:59 

25/HAY 

13:00 

23:45 

1-5/DEC 

13:00 

00:59 

wscni 

18-23/JUL 

13:10 

01:09 

WSC189 

6-12/DEC 

13:00 

00:59 

WSC172 

24-30/JUL 

13:10 

01:09 

WSC190 

13-19/DEC 

13:00 

00:59 

WSC173 

21/ AUG 

21:10 

01:09 

22/AUG 

13:10 

01:09 

23/AUG 

13:10 

14:24 

WSC174 

23/AUG 

16:20 

01:09 

24-29/AUG 

13:10 

01:09 

WSC175 

30-31/AUG 

13:10 

01:09 

1-3/SEP 

13:10 

01:09 

4-5/SEP 

13:00 

00:59 

WSC176 

6-8/SEP 

13:00 

00:59 

9/SEP 

13:00 

13:09 

W3C177 

10/SEP 

14:10 

00:59 

11-16/SEP 

13:00 

00:59 

lf3Cl78 

17/SEP 

14:20 

00:59 

18-23/SEP 

13:00 

00:59 

WSC179 

30/SEP 

13:00 

00:59 

1-6/OCT 

12:50 

00:49 

W3C180 

7-9/OCT 

12:50 

00:49 

lO/OCT 

12:50 

18:06 

WSC181 

lO/OCT 

18:20 

00:49 

11-16/OCT 

12:50 

00:49 

WSC182 

17-23/OCT 

12:50 

00:49 

WSC183 

24-30/OCT 

12:50 

00:49 

WSC184 

31/OCT 

12:50 

00:49 

1-6/NOV 

12:50 

00:49 

WSC185 

7-13/NOV 

12:50 

00:49 

WSC186 

14/NOV 

15:30 

00:49 

15-20/NOV 

12:50 

00:49 

W3C187 

22/NOV 

15:50 

00:49 

23-28/NOV 

12:50 

00:49 
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2.2  White  Sands  Missile  Range/HELSTF ,  New  Mexico 


TAPE  t  DATE 

■START  STOP 

TAPE  •  DATE 

START  STOP 

YEAR  =1988 

WSH021  18/JUL 

16:00  01:19 

19/JUL 

13:20  01:18 

WSHOOl  29/MAR 

16:00  18:00 

20/JUL 

13:20  15:30 

W3H002  30-31/MAR 

06:10  18:09 

WSH022  20/JUL 

17:10  01:19 

1-4/APR 

06:10  18:09 

21-26/JUL 

.  13:20  01:19 

WSH004  11/APR 

16:50  17:51 

WSH023  27/JUL 

14:20  23:20 

12/APR 

06:10  08:42 

WSH024  28/JUL 

19:00  01:19 

WSH005  12/APR 

11:30  18:09 

29/JUL 

13:20  21:00 

13-18/APR 

06:10  18:09 

WSH025  29/JUL 

23:20  01:19 

WSH008  2/MAY 

09:10  18:09 

30/JUL 

13:20  17:36 

3-8/MAY 

06:10  18:09 

WSH026  1/AUG 

15:00  01:19 

WSH009  9/MAY 

08:20  18:09 

2/AUG 

13:20  00:12 

10-14/HAY 

06:10  18:09 

15/MAY 

06:20  12:12 

WSH027  3/AUG 

15:40  00:40 

4 /AUG 

13:20  20:23 

IfSHOlO  18/MAY 

08:20  18:09 

19-24/MAY 

06:10  18:09 

WSH028  5/AUG 

19:20  01:00 

6-7/AUG 

13:20  01:19 

IfSHOll  25/MAY 

07:10  18:09 

8-9/AUG 

13:20  00:50 

26-31/MAY 

06:10  18:09 

irSH029  12/AUG 

20:10  01:09 

WSH012  1/JUN 

07:10  18:09 

13/AUG 

13:10  13:09 

2-7/JUN 

06:10  18:09 

14-15/AUG 

13:10  01:09 

MSHOn  8/JUN 

08:30  12:28 

WSH030  16/AUG 

14:10  00:40 

17/AUG 

13:10  01:09 

;_ALL  DAYS  FROM  THIS  POINT  ON 

18/AUG 

13:10  13:10 

ARE  ON  ZULU-TIME 

WSH031  18/AUG 

15:00  01:09 

WSH014  9/JUN 

21:30  01:09 

19/AUG 

14:10  01:09 

10-15/JUN 

13:10  01:08 

20/AUG 

13:10  01:09 

21/AUG 

13:10  00:30 

WSH015  16/JUN 

14:10  01:09 

22-24/AUG 

13:10  01:09 

17/JUN 

13:10  13:38 

RSH032  26/AUG 

18:20  01:09 

»rSH017  29/JUN 

15:20  01:09 

27/AUG 

13:10  23:32 

30/JUN 

13:10  01:08 

WSH033  29/AUG 

16:10  01:09 

WSH018  5/JUL 

16:30  20:46 

30-31/AUG 

13:10  01:09 

1/SEP 

13:10  01:09 

WSH019  6/JUL 

14:00  01:09 

2/SEP 

13:10  15:10 

7-8/JUL 

13:10  01:08 

WSH034  5/OCT 

21:10  00:59 

WSH020  11/JUL 

18:50  01:09 

6-11/OCT 

13:00  00:59 

12-17/JUL 

13:20  01:18 

WSH035  12/OCT  20:00  00:59 
13-18/OCT  13:00  00:59 


WliiLe  S;imis  Mi'^silp  Kat)f>,o/IIF.I„STI' ,  New  Moxico 


APPKNDIX  A 


lAEUL 

DATE 

START 

STOP 

TAPE  1  DATE 

START 

STOP 

H3H036 

19/OCT 

16:40 

00:49 

WSH051  16/FEB 

16:50 

01:19 

20/OCT 

12:50 

00:49 

17-21/FEB 

13:20 

01:19 

22/FEB 

13:20 

10:19 

WSH037 

26/OCT 

18:30 

00:49 

27-31/OCT 

12:50 

00:49 

WSH052  27/FEB 

19:40 

01:19 

1/NOV 

12:50 

00:49 

28/FEB 

13:20 

01:19 

1/MAR 

13:30 

01:08 

irSH038 

3/NOV 

16:10 

00:49 

WSH053  2/MAR 

17:00 

01:19 

WSH039 

23/NOV 

20:50 

00:59 

3-8/MAR 

13:20 

01:19 

24-29/NOV 

13:00 

00:59 

WSH054  10/MAR 

15:50 

01:19 

WSH040 

30/NOV 

17:10 

00:59 

11-16/MAR 

13:20 

01:19 

1-6/DEC 

13:00 

00:59 

WSH055  17/MAR 

16:20 

01:19 

NSH041 

7/DEC 

16:20 

00:59 

18-20/MAR 

13:20 

01:19 

8-13/DEC 

13:00 

00:59 

WSH056  21/MAR 

17:30 

01:19 

MSH042 

14/DEC 

16:20 

00:59 

22/MAR 

13:20 

01:19 

15-16/DEC 

13:00 

00:59 

23/MAR 

13:20 

23:10 

17-20/DEC 

13:10 

01:09 

WSH057  30/MAR 

21:40 

01:09 

VSH043 

21/DEC 

15:30 

01:09 

31/MAR 

13:10 

01:09 

22/DEC 

13:10 

20:44 

1-5/APR 

13:10 

01:09 

WSH058  6/APR 

16:30 

01:09 

Jlftl 

7-11/APR 

13:10 

01:09 

12/APR 

13:10 

19:09 

W3H044 

3/JAN 

20:50 

01:09 

WSH059  13/APR 

15:40 

01:09 

4-5/JAN 

13:10 

01:09 

14-15/APR 

13:10 

01:11 

6/JAN 

13:10 

15:09 

16/APR 

14:30 

01:09 

17/APR 

13:10 

01:09 

WSH045 

6/JAN 

17:10 

01:09 

18/APR 

13:50 

01:09 

7-8/JAN 

13:10 

01:09 

19/APR 

13:10 

01:09 

9-10/JAN 

13:20 

01:19 

11/JAN 

13:20 

18:10 

WSH060  20/APR 

15:30 

01:09 

21/APR 

13:10 

01:20 

WSH046 

12/JAN 

19:30 

01:19 

22-26/APR 

13:10 

01:09 

13-18/JAN 

13:20 

01:19 

WSH061  28/APR 

15:40 

16:00 

WSH047 

19/JAN 

17:00 

01:19 

20-25/JAN 

13:20 

01:19 

WSH062  2/MAY 

16:00 

01:09 

3/MAY 

13:10 

21:19 

VSH048 

26/JAN 

18:10 

01:19 

27-31/JAN 

13:20 

01:19 

WSH063  5/MAY 

15:00 

01:09 

1/FEB 

13:20 

01:19 

6/MAY 

13:10 

01:09 

7/MAY 

13:10 

23:09 

WSH049 

2/FEB 

16:20 

01:19 

3-0/FEB 

13:20 

01:19 

WSH064  8/MAY 

15:20 

01:09 

9/MAY 

13:10 

01:09 

♦fSHOSO 

9-15/FEB 

13:20 

01:19 

A  I 


APPENDIX  A 


Wlilte  Sands  Missile  Range/HELSTF,  New  Mexico 


IA££JL 

■  III  ■■ 

NSH065 

10/MAY 

15:20 

00:42 

HSH066 

18/MAY 

18:20 

01:09 

19/MAY 

13:10 

01:09 

20/MAY 

13:10 

15:09 

WSH067 

22/MAY 

15:40 

01:09 

23/MAY 

13:30 

01:09 

24/MAY 

15:30 

01:09 

WSH068 

25/MAY 

16:20 

01:09 

26-29/HAY 

13:10 

01:09 

30/HAY 

13:10 

19:00 

WSH069 

1/JUN 

15:10 

01:09 

2-7/JUN 

13:10 

01:09 

HSH070 

12/JUN 

15:10 

01:09 

13/JUN 

13:10 

01:09 

14/JUN 

13:10 

19:54 

WSH071 

14/JUM 

22:30 

01:09 

15/JUN 

13:10 

01:09 

16/JUN 

13:10 

22:41 

WSH072 

19/JUN 

15:20 

01:09 

20-25/JUN 

13:10 

01:09 

WSH073 

26/JUN 

15:20 

01:09 

27/JUN 

13:10 

21:18 

WSH074 

28/JUN 

16:40 

01:09 

29-30/JUN 

13:10 

01:09 

1-3/JUL  13:10  19:09 

4/JUL  13:10  20:09 

WSH075  5/JUL  15:30  01:09 
6-7/JUL  13:10  01:09 

8/JUL  13:10  21:09 
9/JUL  13:10  22:59 

WSH076  11/JUL  15:40  01:09 
12/JUL  13:10  16:54 

WSH077  12/JUL  18:30  01:09 
13-14/JUL  13:10  01:09 
15/JUL  13:10  00:09 
16-18/JUL  13:10  01:09 

WSH078  19/JUL  15:10  01:09 
20-25/JUl,  13:10  01:09 


WSH079  26/JUL  15:10  01:09 
27-31/JUL  13:10  01:09 
1/AUG  13:10  01:09 

WSH080  2/AUG  16:20  01:09 
3-8/AUG  13:10  01:09 

WSH081  9/AUG  14:30  01:09 
10/AUG  13:10  01:09 
11/AUG  13:10  00:59 

WSH082  15/AUG  16:20  01:09 
16/AUG  13:10  01:09 
17/aUG  13:48  01:09 
18-21/AUG  13:10  01:09 

WSH083  22/AUG  19:30  01:09 
23/AUG  13:10  16:16 

VISH084  24/AUG  17:40  19:11 

WSH085  24/AUG  19:14  01:09 


25-30/AUG 

13:10 

01:09 

HSM086 

31/AUG 

15:50 

01:09 

1-2/SEP 

13:10 

01:09 

3/SEP 

13:10 

22:27 

WSH087 

5/SEP 

15:50 

00:59 

6-8/SEP 

13:00 

00:59 

9/SEP 

13:00 

00:19 

WSH088 

12/SEP 

15:20 

00:59 

13-18/SEP 

13:00 

00:59 

WSH089 

19/SEr 

15:10 

00:59 

20-25/SEP 

13:00 

00:59 

WSH090 

26/SEP 

15:20 

00:59 

27-28/SEP 

13:00 

00:59 

29/SEP 

13:00 

13:05 

WSH091 

3/OCT 

13:50 

18:31 

WSH092 

3/OCT 

18:40 

00:49 

4-5/OCT 

12:50 

00:49 

6/OCT 

12:50 

16:08 

W3H093 

6 /OCT 

16:20 

00:49 

7-12/OCT 

12:50 

00:49 

APrKNlHX  A 


Will  to  Snnds  Missilr  RnnRe/liri.STl',  Nrw  Moxiro 


lAEE.i  PATE _ 

_2IARI__5T0P 

XAPi:,l 

DAIE 

START 

9Toe 

W3H094 

13/OCT 

15:30 

00:49 

WSH108 

3/JAN 

16:30 

01:08 

H-18/OCT 

12:50 

00:49 

4-9/JAN 

13:10 

01:08 

19/OCT 

12:50 

16:27 

WSH109 

17/JAN 

17:10 

01:18 

WSH095 

20/OCT 

19:50 

00:49 

18-23/JAN 

13:20 

01:18 

21-26/OCT 

12:50 

00:49 

WSHllO 

10/JAN 

15:40 

01:08 

lfSH096 

27/OCT 

15:40 

00:49 

11-15/JAN 

13:10 

01:08 

28-31/OCT 

12:50 

00:49 

16/JAN 

13:20 

01:18 

1-2/NOV 

12:50 

00:49 

WSHlll 

24/JAN 

16:10 

01:18 

WSH097 

3/NOV 

16:00 

00:49 

25-30/JAN 

13:20 

01:18 

4-9/NOV 

12:50 

00:49 

WSH112 

31/JAN 

15:20 

01:18 

WSH098 

lO/NOV 

17:40 

00:49 

1-6/FEB 

13:20 

01:18 

11-12/NOV 

12:50 

00:49 

13/NOV 

12:50 

21:29 

WSH113 

7/FEB 

16:10 

01:18 

8-9/FEB 

13:20 

01:18 

WSH099 

15-21/NOV 

12:50 

00:49 

10/FEB 

13:20 

21:22 

WSHIOO 

22/NOV 

16:50 

00:49 

WSH114 

12/FEB 

16:10 

01:18 

23-28/NOV 

12:50 

00:49 

13-18/FBB 

13:20 

01:18 

WSHlOl 

29/NOV 

15:10 

00:49 

WSH115 

20/FEB 

16:20 

01:18 

30/NOV 

13:00 

00:59 

21-26/FEB 

13:20 

01:18 

1-3/DEC 

13:00 

00:59 

4/DEC 

13:00 

23:28 

WSH116 

28/FEB 

15:50 

01:18 

1-6/MAR 

13:20 

01:18 

WSH102 

5/DEC 

13:00 

17:40 

WSH117 

7/MAR 

16:10 

01:18 

VSH103 

6/DEC 

22:30 

00:58 

8-13/MAR 

13:20 

01:18 

7-12/DEC 

13:00 

00:58 

WSH118 

14/MAR 

17:20 

01:18 

VSH104 

13/DEC 

15:20 

00:58 

15-20/MAR 

13:10 

01:08 

14-19/DEC 

13:00 

00:58 

WSH119 

21/MAR 

16:10 

01:08 

VSH105 

20/DEC 

16:00 

00:58 

22-26/MAR 

13:10 

01:08 

21-22/DEC 

13:00 

00:58 

27/MAR 

13:10 

17:01 

23-26/OEC 

13:10 

01:08 

WSH120 

27/MAR 

20:20 

01:08 

VSH106 

27/DEC 

15:39 

18:15 

28-31/MAR 

13:10 

01:08 

1-2/APR 

13:10 

01:08 

irsHio? 

27/DEC 

19:50 

01:08 

28-31/DEC 

13:10 

01:08 

WSH121 

3/APR 

16:30 

01:08 

4-9/APR 

13:10 

01:08 

tEAR..= 

■19?  Q. 

WSH122 

10/APR 

14:50 

01:08 

11-16/APR 

13:10 

01:08 

irsHio? 

1-2/JAN 

13:10 

01:08 

VSH123 

17/APR 

14:30 

01:08 

18-23/APR 

13:00 

00:58 

APPENDIX  A 


White  Sands  Missile  Rnnge/HELSTF,  New  Mexico 


lAfE  i  DATE _ START  STOP 


WSH124 

24/APR 

15:50 

00:58 

25-26/APR 

13:00 

00:58 

27/APR 

13:00 

22:06 

WSH125 

30/APR 

14:40 

22:17 

ffSH126 

1/MAY 

16:30 

00:58 

2/HAY 

13:00 

00:58 

3/MAY 

13:00 

17:39 

WSH127 

10/HAY 

16:10 

00:58 

11-15/MAY 

13:00 

00:58 

16/HAY 

13:00 

16:59 

»rSH128 

17/MAY 

16:20 

00:58 

18-20/MAY 

13:00 

00:58 

WSH129 

24/HAY 

19:40 

00:58 

25/MAY 

13:00 

15:24 

WSHUO 

25/MAY 

20:00 

00:59 

26-31/MAY 

13:00 

00:59 

«SH13l  S/JUH 

6-8/JUN 

9/JUN 

10/JUN 


16:50  00:5$ 
13:00  00:69 
13:10  01:09 
13:10  14:09 

16:00  01:09 
13:10  01:09 


WSH132  11/JUN 

12-17/JUN 


WSH133  18/JUN  14:50  01:09 
19-24/JUN  13:10  01:09 


WSHllS  9/JUL 
10/JUL 
11/JUL 

IISH136  2/JUL 
3-6/JUL 
7/JUL 
8/JUL 


17:20  01:09 
13:10  01:09 
13:10  15:36 

16:10  01:09 
13:10  01:09 
14:07  01:09 
13:10  01:09 

20:10  01:09 
13:10  01:09 


14SH137  11/JUL 

12-17/JUL 


TAPE  I 

DATE 

START 

STOP 

WSH140 

27/JUL 

16:10 

01 : 09 

28-31/JUL  13:10 

01:09 

1-2/AUG 

13:10 

01:09 

WS!il41 

3/AUG 

14:40 

17:23 

WSH142 

3/AUG 

19:10 

01:09 

4-7/AUG 

13:10 

01:09 

8/AUG 

13:10 

14:10 

WSH143 

8/AUG 

14:20 

01:09 

9-14/AUG 

13:10 

01:09 

WSM144 

15/AUG 

15:40 

01:09 

16/AUG 

13:10 

01:09 

17/ AUG 

13:10 

15:53 

WSH145 

20/AUG 

14:30 

01:09 

21-26/AUG 

13:10 

01:09 

VSH146 

27/AUG 

16:10 

01:09 

28-31/AUG 

13:10 

01:09 

1-2/SEP 

13:10 

01:09 

WSH147 

5/SEP 

15:50 

00:59 

6-8/SEP 

13:00 

00:59 

WSH148 

10/SEP 

17:50 

00:59 

11-16/SEP 

13:00 

00:59 

WSH149 

18/SEP 

22:10 

00:59 

19-24/SBP 

13:00 

00:59 

WSH150  25/SEP  16:30  00:59 
26-27/SEP  13:00  00:59 
28/SEP  13:00  18:46 


WSH151  28/SEP  19:10  00:59 
29-30/SEP  13:00  00:59 
1-4/OCT  12:50  00:49 


WSH152  5/OCT  16:40  00:49 

6-8/OCT  12:50  00:49 

9/OCT  12:50  15:21 


WSH138  23/JUL  14:50  01:09 
24-26/JUL  13:10  01:09 


WSH153  9/OCT  15:50  00:49 
10-15/OCT  12:50  00:49 


lfSH139  18/JUL  15:30 

19-22/JUL  13:10 
23/JUL  13:10 


WSH154 

01:09 

01:09 

14:39 


16/OCT  16:30  00:49 
17-22/OCT  12:50  00:49 


AI'PENIM':  A 


White  Sands  Missile  Range/IIELSTK,  New  Mexico 


TAPE  i  DATE 


START 


W3H155  23/OCT 

14:50 

00:49 

24-29/OCT 

’  12:50 

00:49 

WSH156  30/OCT 

18:40 

00:49 

31/OCT 

12:50 

00:49 

1/NOV 

12:50 

00:49 

2/NOV 

12:50 

00:29 

W3H157  5/NOV 

16:30 

00:49 

6-8/NOV 

12:50 

00:49 

W3H158  13/NOV 

17:40 

00:49 

14/NOV 

12:50 

19:29 

WSH159  14/NOV 

21:20 

00:49 

15/NOV 

12:50 

00:49 

16/NOV 

12:50 

17:08 

W3H160  16/NOV 

17:40 

00:49 

17-19/NOV 

12:50 

00:49 

20/NOV 

12:50 

21:42 

WSH161  20/NOV 

22:00 

00:49 

21-26/NOV 

12:50 

00:49 

WSH162  27/NOV 

16:40 

00:49 

28/NOV 

12:50 

00:49 

29-30/NOV 

13:00 

00:59 

1-3/DEC 

13:00 

00:59 

WSH163  4/DEC 

16:40 

00:59 

5-10/DEC 

13:00 

00:59 

NSH164  11/DEC 

16:00 

00:59 

12-17/DEC 

13:00 

00:59 

N3H165  18/DEC 

16:10 

00:59 

19-21/DEC 

13:00 

00:59 

22-24/DEC 

13:10 

01:09 

W3H166  26/DEC  16:00  01:09 
27-31/DEC  13:10  01:09 


'45 


APPENDIX  A 


2.3  Kirlland  Air  Force  Base,  Albuquerque,  New  Mexico 


TAPE  t 

DATE 

START 

STOP 

TAPE  1 

DATE 

START 

STOP 

YEAR  = 

JIM 

KAAOie 

17/AUG 

13:50 

01:09 

18-23/AUG 

13:10 

01:09 

KAAOOl 

17/MAY 

19:30 

01:08 

18-23/MAY 

13:10 

01:08 

KAA019 

24/AUG 

13:10 

20:30 

KAA002 

24-30/MAY 

13:10 

01:08 

KAA020 

25/AUG 

14:10 

01:09 

26-31/AUG 

13:10 

01:09 

KAA003 

31/MAY 

13:10 

01:08 

1-2/JUN 

13:10 

01:08 

KAA021 

1/SEP 

13:10 

20:50 

3/JUN 

13:10 

KAA022 

1/SEP 

21:20 

01:09 

KAA004 

3/JUN 

16:40 

01:08 

2-7/SEP 

13:10 

01:09 

4-9/JUN 

13:10 

01:08 

KAA023 

8/SEP 

13:50 

01:09 

KAA005 

10-13/JUN 

13:10 

01:08 

9/SEP 

13:10 

16:30 

14/JUN 

13:10 

KAA024 

12/SEP 

14:00 

01:09 

KAA006 

14-20/JUN 

13:10 

01:08 

13/SEP 

13:10 

16:30 

KAA007 

21-27/JUN 

13:10 

01:08 

KAA025 

14/SEP 

16:10 

01:09 

15/SEP 

13:10 

13:34 

KAA008 

28/JUN 

13:10 

18:02 

KAA026 

15/SEP 

13:50 

01:09 

KAA009 

28/JUN 

20:40 

22:40 

16/SEP 

13:10 

13:34 

29-30/JUN 

13:10 

00:10 

17-21/SEP 

13:00 

00:59 

1-4/JUL 

13:10 

00:10 

KAA027 

22/SEP 

13:00 

20:40 

KAAOlO 

5/JUL 

13:10 

01:09 

6-8/JUL 

13:20 

01:18 

KAA028 

23/SEP 

16:30 

19:00 

9/JUL 

13:20 

00:58 

KAA030 

1/OCT 

20:30 

00:59 

KAAOll 

19/JUL 

20:30 

01:19 

2-3/OCT 

13:00 

00:59 

20/JUL 

13:10 

01:19 

4/OCT 

13:00 

16:20 

21/JUL 

13:10 

19:42 

KAA031 

6 /OCT 

16:40 

00:59 

KAA012 

21/JUL 

20:30 

01:19 

7/OCT 

13:00 

00:59 

22-23/JUL 

13:20 

01:19 

8/OCT 

13:00 

18:50 

24/JUL 

13:20 

00:00 

KAA032 

13/OCT 

13:20 

00:59 

KAA013 

25/JUL 

16:40 

01:19 

14-19/OCT 

13:00 

00:59 

26/JUL 

13:20 

17:14 

KAA033 

20/OCT 

13:00 

15:34 

KAAOH 

1/AUG 

20:30 

23:30 

KAA034 

2/NOV 

23:40 

00:49 

KAA013 

2/AUG 

17:30 

01:19 

3-7/NOV 

12:50 

00:49 

3-7/AUG 

13:20 

01:19 

8 /AUG 

14:20 

21:04 

KAA035 

9-13/NOV 

12:60 

00:49 

14-15/NOV 

13:00 

00:59 

KAA016 

9/AUG 

16:20 

01:19 

10-15/AUG 

13:20 

01:19 

KAA036 

16-20/NOV 

13:00 

00:59 

21/NOV 

13:00 

17:35 

KAA017 

16/AUG 

13:20 

23:00 

AI’I’I'NDIX  A 


Kinlnml  Ait  r'n!  i  r  ll.irjc,  A I  l)ii')iii'rqu(’ ,  Nt'w  Mi’xitd 


TAPE  1 

DATE 

START 

STOP 

lAEEJL 

DATE _ 

START 

STOP 

KAA037 

2i/N0V 

22:50 

00:59 

KAA056 

29-30/HAR 

13:20 

01:19 

22-27/NOV 

13:00 

00:59 

31/MAR 

13:10 

01:09 

1-4/APR 

13:10 

01:09 

KAA036 

28/NOV 

21:20 

00:59 

29/NOV 

13:00 

00:59 

KAA057 

5/APR 

13:30 

18:50 

30/NOV 

13:00 

23:59 

KAA05e 

5/APR 

18:40 

01:09 

KAA039 

5-10/DEC 

13:00 

00:59 

6-11/APR 

13:10 

01:09 

11/DEC 

13:00 

00:06 

KAA059 

12/APR 

13:40 

01:09 

KAA040 

12-14/DEC 

13:00 

00:59 

13-15/AF.’: 

13:10 

01:09 

15-17/DEC 

13:10 

01:09 

16/APR 

13:20 

01:09 

17-18/APR 

13:10 

01:09 

KAA041 

19/DEC 

16:00 

01:09 

20-25/DEC 

13:10 

01:09 

KAA061 

20/APR 

14:40 

01:09 

21-26/APR 

13:10 

01:09 

KAA042 

26/DEC 

14:20 

01:09 

27-31/DEC 

13:10 

01:09 

KAA062 

27-30/APR 

13:10 

01:09 

1-3/MAY 

13:10 

01:09 

ms^ 

1.3  9  i 

KAA063 

5-10/MAY 

13:10 

01:09 

KAA064 

11-17/MAY 

13:10 

01:09 

KAA042 

1/JAN 

13:10 

01:09 

KAA065 

18-24/MAY 

13:10 

01:09 

KAA047 

30/JAN 

20:10 

01:19 

31/JAN 

13:20 

01:19 

KAA066 

25/MAY 

13:20 

01:09 

1-5/FEB 

13:20 

01:19 

26/MAY 

13:10 

01:09 

27/MAY 

13:10 

13:30 

KAA048 

6/FEB 

13:30 

01:19 

7/FEB 

13:20 

17:22 

KAA067 

27/MAY 

23:20 

01:09 

28-31/MAY 

13:10 

01:09 

KAA049 

9/FEB 

18:00 

01:19 

1-2/JDN 

13:10 

01:09 

10-15/FEB 

13:20 

01:19 

KAA068 

3-9/JUN 

13:10 

01:09 

KAA050 

16-20/FEB 

13:20 

01:19 

21/FEB 

13:20 

16:20 

KAA069 

10-14/JUN 

13:10 

01:09 

15/JUN 

13:10 

20:17 

KAA051 

22/FEB 

16:40 

01:19 

23-28/FEB 

13:20 

01:19 

KAA070 

15/JUN 

22:40 

01:09 

16-20/JUN 

13:10 

01:09 

KAA052 

1-7/MAR 

13:20 

01:19 

21/JUN 

13:10 

20:39 

KAA053 

8-14/MAR 

13:20 

01:19 

KAA072 

29/JUN 

15:20 

01:09 

30/JUN 

13:10 

01:09 

KAA054 

15-19/MAR 

13:20 

01:19 

1-5/JUL 

13:10 

01:09 

20/MAR 

13:20 

22:48 

KAA073 

6/JUL 

13:10 

20:30 

KAA0S5 

22-28/MAR 

13:20 

01:19 

KAA075 

14/JUL 

22:20 

01:09 

15-18/JUL 

13:10 

01:09 

19/JUL 

13:10 

18:46 

APPENDIX  A 


K  i  r  1 1 ,11 

TAPE  1 

ul  Air  I'orce  Base,  Albtifjiicrqiic, 

DATE _ START  STOP 

New  ilex  Leo 

TAPE  1 

DATE 

START 

STOP 

KAA076 

7/JUL 

23:00 

01:09 

KAA094 

2/NOV 

12:50 

16:41 

8-13/JUL 

13:10 

01:09 

KAA095 

2/NOV 

18:10 

00:49 

KAA077 

19/JUL 

22:00 

01:09 

3-8/NOV 

12:50 

00:49 

20-25/JUL 

13:10 

01:09 

KAA096 

10-16/NOV 

12:50 

00:49 

KAA078 

26-31/JUL 

13:10 

01:09 

1/AUG 

13:10 

01:09 

KAA097 

17-23/NOV 

12:50 

00:49 

KAA079 

2-8/AUG 

13:10 

01:09 

KAA098 

24-26/NOV 

12:50 

00:49 

27-29/NOV 

13:00 

00:59 

KAAO0O 

9-15/AUG 

13:10 

01:09 

30/NOV 

13:00 

00:05 

KAA081 

16-22/AUG 

13:10 

01:09 

KAA099 

1-6/DEC 

13:00 

00:59 

KAA082 

23/AUG 

13:20 

01:09 

KAAIOO 

7/DEC 

23:20 

00:59 

24-27/AUG 

13:10 

01:09 

8-13/DEC 

13:00 

00:59 

28/AUG 

13:10 

16:40 

KAAlOl 

14/DEC 

23:50 

00:59 

KAA083 

29/AUG 

17:00 

01:09 

15-20/DEC 

13:00 

00:59 

30-31/AUG 

13:10 

01:09 

1-4/SEP 

13:10 

01:09 

KAA102 

21-26/DEC 

13:10 

01:09 

27/DEC 

13:10 

22:21 

KAA084 

5-6/SEP 

13:10 

01:09 

7-11/SEP 

13:00 

00:59 

KAA103 

27/DEC 

22:50 

01:09 

28-31/DEC 

13:10 

01:09 

KAA085 

12/SEP 

13:00 

19:28 

KAA086 

13/SEP 

21:30 

00:59 

YEAR  = 

1990 

14-19/SEP 

13:00 

00:59 

KAA087 

20-26/SEP 

13:00 

00:59 

KAAlOl 

1/JAN 

13:10 

01:09 

2/JAN 

13:10 

00:06 

KAA088 

27-30/SEP 

13:00 

00:59 

1-3/OCT 

13:00 

00:59 

KAA104 

3-8/JAN 

13:10 

01:09 

KAA089 

4/OCT 

13:20 

00:59 

KAA105 

9/JAN 

14:20 

01:09 

5/OCT 

12:50 

13:49 

10-12/JAN 

13:10 

01:09 

13-15/JAN 

13:20 

01:19 

KAA090 

5/OCT 

19:50 

00:49 

6-11/OCT 

12:50 

00:49 

KAA106 

16/JAN 

14:10 

01:19 

17-22/JAN 

13:20 

01:19 

KAA091 

12/OCT 

21:20 

00:49 

13-18/OCT 

12:50 

00:49 

KAA107 

6/FEB 

14:00 

01:19 

7-8/FEB 

13:20 

01:19 

KAA092 

19/OCT 

13:20 

00:49 

9/FEB 

13:20 

00:19 

20-25/OCT 

12:50 

00:49 

KAA108 

30/JAN 

14:00 

01:19 

KAA093 

26/OCT 

13:10 

00:49 

31/JAN 

13:20 

01:19 

27-30/OCT 

12:50 

00:49 

1-5/FEB 

13:20 

01:19 

31/OCT 

12:50 

17:49 

<'.8 


API’F.NDIX  a 


KirLiaiiil  Air  I'di'f'  A I  hii<(iip  rriiic ,  New  Mcy.  iin 


TAPE  i  DATE 

-SIARI. 

STOP 

TAPE  I  DATE 

START 

STOP 

KAA109  12/FEB 

23:20 

01:19 

KAA122  20/JUN 

20:10 

01:09 

13-15/FEB 

13:20 

01:19 

21/JUN 

13:10 

21:42 

16/FEB 

13:20 

22:59 

KAA123  2^/JUN 

17:16 

17:31 

KAAllO  16/FEB 

23:10 

01:19 

17/FEH 

13:20 

01:19 

KAA124  18/JUf. 

22:40 

01:09 

18/FEB 

13:20 

19:13 

19-24/JUL 

13:10 

01:09 

KAAlll  20/FEB 

17:50 

01:19 

KAA125  3/AUG 

15:40 

01:09 

21-23/FEB 

13:20 

01:19 

4-8/AUG 

13:10 

01:09 

24/FEB 

13:20 

23:36 

9/AUG 

13:10 

17:00 

KAA112  26/FEB 

15:00 

01:19 

KAA126  26/AUG 

20:40 

01:09 

27-28/FEB 

13:20 

01:19 

27-31/JUL 

13:10 

01:09 

1-4/MAR 

13:20 

01:19 

1/AUG 

13:10 

01:09 

KAA113  5/MAR 

14:20 

01:19 

KAA127  9/AUG 

17:20 

01:09 

6-7/MAR 

13:20 

01:19 

10-15/AUG 

13:10 

01:09 

8/MAR 

13:20 

14:20 

KAA128  20/AUG 

14:20 

01:09 

KAA114  9/MAR 

21:20 

01:19 

21-26/AUG 

13:10 

01:09 

10-15/MAR 

13:20 

01:19 

KAA129  27/AUG 

14:50 

01:09 

KAA115  16/MAR 

14:10 

01:19 

28-31/AUG 

13:10 

01:09 

17-18/MAR 

13:20 

01:19 

1-2/SEP 

13:10 

01:09 

19/MAR 

13:10 

01:09 

20/MAR 

13:10 

22:33 

KAA130  4/SEP 

18:50 

01:09 

21/MAR 

13:10 

17:28 

5/SEP 

13:10 

01:09 

6-7/SEP 

13:00 

00:59 

KAA116  6/APR 

21:20 

01:09 

8/SEP 

13:00 

15:00 

7-9/APR 

13:10 

01:09 

10/APR 

13:10 

19:06 

KAA131  12/SEP 

17:30 

00:59 

11/APR 

13:10 

20:12 

13-18/SEP 

13:00 

00:59 

12/APR 

13:10 

00:34 

KAA132  20/SEP 

14:30 

00:59 

KAA117  17/MAY 

22:20 

00:59 

21-26/SEP 

13:00 

00:59 

18/MAY 

13:00 

16:07 

KAA133  1/OCT 

15:40 

00:59 

KAA118  18/MAY 

17:20 

00:59 

2-3/OCT 

13:00 

00:59 

19-24/MAY 

13:00 

00:59 

4 /OCT 

12:50 

00:49 

5/OCT 

12:50 

17:07 

KAA119  29/MAY 

16:00 

00:59 

30-31/MAY 

13:00 

00:59 

KAA134  11/OCT 

14:50 

00:49 

1-2/JUN 

13:00 

00:59 

12-16/OCT 

12:50 

00:49 

3-4/JUN 

13:10 

01:09 

17 /OCT 

12:50 

14:16 

KAA120  6/JUN 

15:20 

01:09 

KAA135  17/OCT 

14:30 

00:49 

7-12/JUN 

13:10 

01:09 

18/OCT 

12:50 

16:32 

KAA121  13/JUN 

21:20 

01:09 

KAA136  24/OCT 

13:50 

00:49 

14-19/JL'N 

13:10 

01:09 

25/OCT 

12:5C 

00:49 

26/OCT 

12:50 

20:59 

'.9 


APPENDIX  A 


Kirtland  Air  Force  Base,  Albuquerque,  New  Mexico 


KAA141  11/DEC  21:30  00:59 
12-17/DEC  13:00  00:59 


KAA142  19/DEC  15:10  00:59 
20-25/DEC  13:10  01:09 


AI’PF.NDIX  A 


2.4  Navnl  Weapons  ('enter,  China  hake,  California 


TAPE  I 

DATE 

START 

STOP 

TAPE  1 

DATE 

START 

STOP 

YEAR  = 

1988 

CL4019 

5-8/SEP 

13:50 

01:48 

9/SEP 

13:50 

18:30 

CL4001 

23/JUN 

19:30 

01:08 

CL4023 

4/NOV 

20:40 

01:39 

CL4002 

23/JUN 

23:00 

01:52 

6-7/NOV 

13:40 

01:39 

8/NOV 

13:40 

17:14 

CL4004 

1-4/JUL 

14:00 

01:58 

5/JUL 

14:00 

19:26 

CL4024 

8 /NOV 

21:00 

01:39 

9-14/NOV 

13:40 

01:39 

CL4005 

5/JUL 

19:50 

01:59 

6-11/JUL 

14:00 

01:12 

CL4025 

15-19/NOV 

13:40 

01:39 

20/NOV 

14:50 

01:39 

CL4006 

12-18/JUL 

14:00 

01:12 

22/NOV 

13:40 

01:39 

CL4007 

19/JUL 

19:10 

01:59 

CL4026 

23-24/NOV 

13:40 

01:39 

20/JUL 

14:00 

00:10 

25/NOV 

13:40 

14:59 

26/NOV 

15:00 

01:39 

CL4008 

21-22/JUL 

14:00 

01:12 

CL4027 

28/NOV 

20:10 

01:39 

CL4009 

22/JUL 

23:30 

01:30 

23-28/JUL 

14:00 

01:12 

CL4028 

29/NOV 

17:50 

01:39 

CL4010 

29-31/JUL 

14:00 

01:58 

CL4029 

30/NOV 

19:00 

01:39 

1/AUG 

23:59 

01:58 

2-3/AUG 

14:00 

01:58 

CL4030 

1/DEC 

16:00 

01:39 

4/AUG 

14:00 

19:59 

CL4031 

2 /DEC 

17:30 

01:39 

CL4011 

5-11/AUG 

14:00 

01:58 

CL4032 

3/DEC 

15:20 

01:49 

CL4012 

12-13/AUG 

14:00 

01:59 

4/DEC 

13:50 

01:49 

CL4013 

14-15/AUG 

14:00 

01:58 

CL4033 

5/DEC 

15:40 

01:49 

16/AUG 

14:00 

00:58 

17/AUG 

14:00 

01:58 

CL4034 

6/DEC 

18:20 

01:49 

18/AUG 

14:00 

00:58 

CL4035 

7 /DEC 

16:00 

01:49 

CL4014 

19-20/AUG 

14:00 

01:58 

8-9/DEC 

13:50 

01:49 

21/AUG 

14:00 

23:58 

10/DEC 

21:51 

21:00 

22/AUG 

14:00 

00:58 

2 3/ AUG 

14:00 

01:58 

CL4036 

11/DEC 

23:10 

15:58 

CL4015 

24/AUG 

14:00 

01:59 

CL4037 

13/DEC 

16:20 

17:19 

CL4016 

25/AUG 

17:20 

01:59 

CL4038 

13/DEC 

17:30 

01:49 

26/AUG 

14:00 

20:50 

14/DEC 

23:56 

00:36 

16/DEC 

13:50 

01:49 

CL4017 

26/AUG 

21:10 

01:58 

17/DBC 

15:40 

23:27 

27/AUG 

14:00 

01:58 

28/AUG 

14:00 

01:04 

CL4039 

21/DEC 

15:10 

01:49 

22-24/DEC 

13:50 

01:49 

CL4018 

29/AUG 

16:50 

01:59 

30-31/AUG 

14:00 

01:58 

51 


APPENDIX  A 


Naval  Weapons  Center,  China  Lake,  California 


TAPE  •  DATE 

STA8T 

STOP 

TAPE  •  DATE 

START 

STOP 

CL4040  27/DEC 

18:50 

01:59 

CL4047  20/FEB 

18:40 

02:09 

28/DEC 

15:00 

00:51 

21/FEB 

15:10 

02:09 

29/DEC 

14:00 

15:20 

22/FEB 

14:50 

02:09 

23/FEB 

14:15 

02:09 

24-25/FEB 

14:08 

02:09 

YEAR  =  1989 

26/FEB 

02:09 

02:09 

CL4049  6/MAR 

14:45 

02:09 

CL4041  4/JAN 

14:58 

01:59 

7/MAR 

22:50 

02:09 

5/JAN 

14:19 

01:59 

8/MAR 

14:08 

02:09 

6/JAN 

13:58 

01:59 

9-10/MAR 

14:40 

02:09 

7/JAN 

15:43 

00:48 

11 /MAR 

16:50 

02:09 

8/JAN 

16:24 

00:48 

9/JAN 

15:11 

10:50 

CL4050  12/MAR 

14:08 

14:10 

10/JAN 

14:00 

01:59 

13/HAR 

13:00 

01:59 

CL4042  11/JAN 

15:20 

01:59 

CL4051  16/MAR 

17:00 

00:15 

12/JAN 

15:46 

06:14 

13/ JAN 

23:00 

07:18 

CL4053  5/APR 

23:20 

01:59 

14/JAN 

15:07 

01:08 

6-9/APR 

13:58 

01:59 

15/ JAN 

15:07 

08:01 

10/APR 

21:29 

21:29 

CL4043  17/JAN 

13:58 

01:59 

CL4054  11/APR 

18:40 

01:59 

18-19/JAN 

15:18 

01:59 

12-14/APR 

14:00 

01:59 

15-17/APR 

13:50 

01:49 

CL4044  20/JAN 

16:04 

01:59 

21/JAN 

15:20 

02:09 

CL4055  17-23/APR 

13:50 

01:49 

22/JAN 

15:50 

02:09 

23/JAN 

15:34 

02:09 

CL4056  24/APR 

13:50 

01:49 

24/Jan 

15:20 

02:09 

25/JAN 

11:40 

11:40 

CL4057  25/APR 

13:50 

01:49 

26/JAN 

14:10 

02:09 

26/APR 

13:50 

21:02 

CL4045  27/JAN 

14:08 

06:01 

CL4058  29/APR 

20:10 

01:49 

28/JAN 

15:40 

02:09 

30/APR 

13:50 

01:49 

29/JAN 

15:56 

02:09 

1-5/MAY 

13:50 

01:49 

30/JAN 

15:46 

02:09 

31/JAN 

02:09 

02:09 

CL4059  6/MAY 

13:50 

01:49 

1/FEB 

15:10 

02:09 

7/MAY 

13:50 

21:23 

2/FEB 

14:20 

02:09 

CL4061  15-20/HAY 

13:50 

01:49 

CL4046  3/FEB 

14:08 

02:09 

4/FEB 

15:36 

01:05 

CL4062  21-27/MAY 

13:50 

01:49 

6/FEB 

14:56 

02:09 

7/FEB 

15:57 

00:28 

CL4063  31/MAy 

17:00 

01:49 

8/FEB 

15:14 

04:17 

1-5/JUN 

13:50 

01:49 

9/FEB 

14:10 

02:09 

11/FEB 

15:04 

01:28 

CL4064  7/JUN 

17:00 

01:49 

8-13/JUN 

13:50 

01:49 

Al’PKNDIX  A 


Navnl  Weapons  Oiilcr,  Chinn  I.ako,  California 


TAPE  1 

DATE 

START 

STOP 

TAPE  1  DATE 

START 

STOP 

CL4065 

14-16/JUN 

13:50 

23:59 

CL40ei  6/OCT 

22:30 

01:39 

17/JUN 

13:50 

01:49 

7-8/OCT 

13:40 

01:39 

18-20/JUN 

13:50 

23:59 

9 /OCT 

13:40 

23:39 

lO-ll/OCT 

13:40 

01:38 

CL4066 

23-29/JUN 

13:50 

01:49 

12/OCT 

13:40 

21:05 

CL4067 

30/JUN 

20:00 

01:49 

CL4082  12/OCT 

23:10 

01:39 

1-4/JUL 

14:00 

01:59 

13/OCT 

13:40 

14:50 

5/JUL 

14:00 

21:45 

CL4083  13/OCT 

16:10 

01:39 

CL4068 

9/JUL 

17:10 

01:59 

14 /OCT 

13:40 

01:39 

10-15/JUL 

14:10 

01:58 

15/OCT 

15:11 

22:56 

CL4069 

18-24/JUL 

14:00 

01:58 

CL4084  17/OCT 

20:19 

21:19 

CL4070 

25-31/JUL 

14:00 

01:58 

CL4085  18/OCT 

21:30 

01:39 

19 /OCT 

13:40 

15:49 

CL4071 

1/AUG 

15:10 

01:58 

2-7/AUG 

14:00 

01:58 

CL4086  20/OCT 

16:50 

01:39 

21/OCT 

13:40 

01:39 

CL4072 

8/AUG 

17:30 

01:58 

22/OCT 

13:40 

15:02 

9-14/AUG 

14:00 

01:58 

CL4087  24/OCT 

23:50 

01:39 

CL4073 

15/AUG 

15:00 

01:59 

25/OCT 

13:40 

19:46 

16-19/AUG 

14:00 

01:59 

20-21/AUG 

13:50 

01:49 

CL4088  25/OCT 

23:00 

01:39 

26/OCT 

13:40 

01:39 

CL4074 

22-28/AUG 

13:50 

01:48 

27-31/OCT 

13:30 

01:29 

CL4075 

29/AUG 

15:10 

01:49 

CL4089  1/NOV 

13:30 

15:07 

30-31/AUG 

13:50 

01:49 

1-4/SEP 

13:50 

01:48 

CL4090  11/NOV 

23:40 

01:39 

12-15/NOV 

13:40 

01:39 

CL4076 

5/SEP 

14:10 

01:49 

16/NOV 

13:40 

17:21 

6-8/3EP 

13:50 

01:48 

9/SEP 

14:30 

01:48 

CL4091  20/NOV 

18:00 

01:39 

10-11/SEP 

13:50 

01:48 

21-26/NOV 

13:40 

01:39 

CL4077 

12-18/SEP 

13:50 

01:49 

CL4092  27/NOV 

17:10 

01:39 

28-30/NOV 

13:40 

01:39 

CL4078 

19/SEP 

14:20 

01:39 

1-3/DEC 

13:40 

01:39 

20-22/3EP 

13:40 

01:39 

23/SEP 

13:53 

01:39 

CL4093  4/DEC 

18:30 

01:39 

24/SEP 

14:00 

01:39 

5-7/DEC 

13:40 

01:39 

25/SEP 

13:40 

01:39 

8/DEC 

13:40 

21:35 

CL4079 

26/SEP 

17:30 

01:39 

CL4094  10/DEC 

20:50 

01:39 

27-28/SEP 

13:40 

01:39 

11/DEC 

13:40 

18:00 

29/SEP 

13:40 

20:14 

CL4095  12/DEC 

23:50 

01:49 

CL4080 

1-4/OCT 

13:40 

01:39 

13/DEC 

13:50 

01:49 

14/DEC 

13:50 

16:38 

APPENDIX  A 


Naval  Weapons  Center,  China  Lake,  California 


iTlwWl 

CL4096 

15-18/DEC 

13:50 

01:49 

CL4112 

29/MAY 

19:10 

01:49 

19/DEC 

13:50 

17:01 

30-31/MAY 

13:50 

01:49 

1-3/JUN 

13:50 

01:49 

CL4097 

21/DEC 

20:20 

01:49 

4/JUN 

13:50 

18:02 

22-27/DEC 

13:50 

01:49 

CL4113 

5-8/JUN 

13:50 

01:49 

CL4098 

28-31/DEC 

13:50 

01:49 

9/JUN 

13:50 

19:57 

CL4114 

12/JUN 

16:00 

01:49 

xm  = 

I??? 

13/JUN 

13:50 

01:49 

14/JUN 

13:50 

16:03 

CL4098 

1/JAN 

13:50 

01:49 

CL4115 

3/JUL 

20:30 

01:59 

2-3/JAN 

14:00 

01:59 

4-9/JUL 

14:00 

01:59 

CL4099 

4/JAN 

21:30 

01:59 

CL4116 

19/JUL 

21:20 

01:59 

5-10/JAN 

14:00 

01:59 

20/JUL 

14:00 

21:59 

CL4100 

11/JAN 

21:20 

01:59 

CL4117 

2/AUG 

20:00 

01:59 

12-17/JAN 

14:00 

01:59 

3 /AUG 

14:00 

16:15 

CL4101 

21/JAN 

21:30 

01:59 

CL4118 

3/AUG 

18:00 

01:59 

22-27/JAN 

14:00 

01:59 

4-9/AUG 

14:00 

01:59 

CL4102 

29/JAN 

21:50 

01:59 

CL4119 

10/AUG 

16:40 

01:59 

30-31/JAN 

14:00 

01:59 

11-16/AUG 

14:00 

01:59 

1-4/FEB 

14:00 

01:59 

CL4120 

17/AUG 

16:30 

01:59 

CL4103 

5-11/FEB 

14:10 

02:09 

18/AUG 

14:00 

01:59 

19-23/AUG 

13:50 

01:49 

CL4104 

7/MAR 

19:40 

01:59 

8-13/MAR 

14:00 

01:59 

CL4121 

26/AUG 

21:20 

01:49 

27-31/AUG 

13:50 

01:49 

CL4105 

14-17/MAR 

14:00 

01:59 

1/SEP 

13:50 

01:49 

18/MAR 

14:00 

21:38 

CL4122 

4-10/SEP 

13:50 

01:48 

CL4106 

18/APR 

17:40 

01:49 

CL4123 

11-17/SEP 

13:50 

01:48 

CL4107 

23-29/APR 

13:50 

01:49 

CL4124 

18-24/SEP 

13:40 

01:39 

CL4108 

1-7/MAY 

13:50 

01:49 

CL4125 

25-30/SEP 

13:40 

01:39 

CL4109 

8/MAY 

22:30 

01:49 

1/OCT 

13:40 

01:39 

9-14/MAY 

13:50 

01:49 

CL4126 

2-7/OCT 

13:40 

01:39 

CL4110 

16/MAY 

16:00 

01:49 

8 /OCT 

23:00 

01:39 

17-18/MAY 

13:50 

01:49 

19/HAY 

13:50 

16:51 

CL4127 

9-15/OCT 

13:40 

01:39 

CL4111 

21-27/MAY 

13:50 

01:49 

CL4128 

16-22/OCT 

13:40 

01:39 

APPENDIX  A 


Naval 

U ('a  pons  Con  tor. 

Cli  i  na 

TAPE  1 

DATE 

START 

STOP 

CLn29 

23-25/OCT 

13:40 

01:39 

26-29/OCT 

13:30 

01:29 

CL4130 

30-31/OCT 

13:30 

01:29 

1-4/NOV 

13:30 

01:29 

5/NOV 

13:40 

01:39 

CL4131 

6-12/NOV 

13:40 

01:39 

CL4132 

13/NOV 

15:40 

01:39 

14-19/NOV 

13:40 

01:39 

CL4133 

20-24/NOV 

13:40 

01:39 

CL4134 

28-30/NOV 

13:40 

01:39 

1-3/DEC 

13:40 

01:39 

CL4135 

5/DEC 

22:00 

01:39 

6/DEC 

13:40 

21:28 

C;i  1  i  font  i  .1 


2.5  Malmstrom  Air  Force  Base,  Montana 


APPENDIX  A 


TAPE  i  DATE 

START  STOP 

TAPE  1  DATE 

START  .STOP 

YEAR  =  1988 

MAG018  6/DEC 

18:40  01:19 

7-12/DEC 

13:20  01:19 

MAGOOl  30/AUG 

20:40  01:29 

31/AUG 

13:30  14:58 

MAG019  13/DEC 

15:10  01:19 

14-17/DEC 

:  13:20  01:19 

MAG002  31/AUG 

16:20  21:10 

18-19/DEC 

;  13:30  01:29 

MAG003  31/AUG 

22:00  01:29 

MAG020  20/DEC 

16:40  01:29 

1-6/AUG 

13:30  01:29 

21-26/DEC 

13:30  01:29 

HAG004  8/SBP 

16:50  01:29 

MAG021  27/DEC 

15:10  01:29 

9-12/SEP 

13:30  01:29 

28-31/DEC 

13:30  01:29 

13-H/SEP 

13:20  01:19 

NAG005  16/SEP 

15:40  01:19 

YEAR  =  1989 

17-22/SEP 

13:20  01:19 

HAG006  23/SBP 

13:50  01:19 

HAG021  1-2/JAN 

13:30  01:29 

24-29/SEP 

13:20  01:19 

MAG022  3/JAN 

15:10  01:29 

MAG007  30/SEP 

14:50  01:19 

4-9/JAN 

13:30  01:29 

1-6 /OCT 

13:20  01:19 

MAG023  IC/JAN 

14:40  01:39 

MAG008  7/OCT 

14:20  01 : 19 

11-16/JAN 

13:40  01:39 

8-12/OCT 

13:20  01:19 

13/OCT 

13:10  01:09 

MAG024  17/JAN 

15:00  01:39 

18-23/JAN 

13:40  01:39 

MAG009  14/OCT 

14:50  01:19 

15-20/OCT 

13:10  01:09 

MAG025  24/JAN 

15:00  01:39 

25-26/JAN 

13:40  01:39 

MAGOlO  21/OCT 

14:20  01:09 

27/JAN 

13:40  19:04 

22-27/OCT 

13:10  01:09 

MAG026  31/JAN 

15:30  01:39 

MAGOll  28/OCT 

15:30  01:09 

1/FEB 

13:40  15:31 

29-31/OCT 

13:10  01:09 

MAG027  1/FEB 

21:00  01:39 

MAG012  1/NOV 

17:50  01:09 

2-6/FEB 

13:40  01:39 

MAG013  9/NOV 

23:30  01:09 

MAG028  7/FEB 

14:30  01:39 

10-15/NOV 

13:10  01:09 

8-13/FEB 

13:40  01:39 

MAG014  16/NOV 

18:10  01:09 

MAG029  14/FEB 

14:40  01:39 

17/NOV 

13:10  01:09 

15-20/FEB 

13:40  01:39 

18/NOV 

13:10  13:46 

MAG030  22/FBB 

14:40  01:39 

MAG016  21/NOV 

15:30  01:19 

23-28/FEB 

13:40  01:39 

22-23/NOV 

13:20  01:19 

24/NOV 

13:20  16:51 

MAG031  1/MAR 

15:30  01:39 

2-7/MAR 

13:40  01:39 

MAG017  29/NOV 

22:00  01:19 

30/NOV 

13:20  01:19 

MAG032  8/MAR 

15:30  01:39 

1-5/DEC 

13:20  01:19 

9-14/MAR 

13:40  01:39 

APPENDIX  A 


M.ilms 

iAeiL_t 

In 'in  Air  1' 

DATE 

on  ('  11; 

■  START 

ISP,  Mi'll t.ina 

__STOP 

TAPE  1 

DATE 

START 

STOP 

HAG033 

15/MAR 

14:40 

01:39 

MAG048 

27-30/JUN  13:30 

01:29 

16-21/MAR  13:40 

01:39 

1-3/JUL 

13:30 

01:29 

HAG034 

22/MAR 

15:10 

01:39 

NAG049 

5/JUL 

16:40 

01:29 

23-25/MAR  13:40 

01:39 

6-11/JUL 

13:30 

01:29 

26-27/MAR  13:30 

01:29 

26/HAR 

13:30 

15:29 

MAG050 

12/JUL 

16:00 

01:29 

13-17/JUL 

,  13:30 

01:29 

HAG035 

29/MAR 

15:50 

01:29 

18/JUL 

13:40 

01:39 

30-31/MAR 

t  13:30 

01:29 

1/APR 

13:30 

01:19 

M AGO 51 

19/JUL 

15:00 

01:39 

20-21/JUL 

13:40 

01:39 

HAG036 

4/APR 

20:50 

01:29 

5-10/APR 

13:30 

01:29 

MAG052 

16/AUG 

17:10 

01:29 

17-22/AUG 

13:30 

01:29 

MAG037 

11/APR 

14:30 

01:29 

12-17/APR 

13:30 

01:29 

MAG053 

23/AUG 

14:40 

01:29 

24-28/AUG 

13:30 

01:29 

HAG038 

18/APR 

14:20 

01:29 

29/AUG 

13:30 

23:41 

19-24/APR 

13:30 

01:29 

HAG0S4 

8/SEP 

13:20 

15:37 

MAG039 

25/APR 

20:30 

01:29 

26-30/APR 

13:30 

01:29 

MAG055 

8 /SEP 

16:50 

01:19 

l/MAY 

13:30 

01:29 

9-14/SEP 

13:20 

01:19 

HAG040 

2-8/HAY 

13:30 

01:29 

NAG0S6 

15/SEP 

13:50 

01:19 

16-21/SEP 

13:20 

01:19 

HAG041 

9/MAy 

14:00 

01:29 

10/MAY 

13:30 

01:29 

MAG057 

22/SEP 

14:30 

01:19 

11-14/NAY 

13:20 

01:19 

23-28/SEP 

13:20 

01:19 

15/MAY 

13:20 

13:46 

HAG058 

29/SEP 

16:50 

01:19 

NAG042 

16/MAY 

16:40 

01:19 

30/SBP 

13:10 

01:09 

17/MAY 

13:20 

01:19 

1/OCT 

13:10 

01:09 

18-22/MAY 

13:30 

01:29 

MAG059 

2/OCT 

14:00 

01:09 

HAG043 

23/NAY 

14:00 

01:29 

3/OCT 

13:10 

01:09 

24-29/MAY 

13:30 

01:29 

MAG061 

11/OCT 

19:20 

01:09 

HAG044 

30/MAY 

14:30 

01:29 

12-17/OCT 

13:10 

01:09 

31/MAY 

13:30 

01:29 

1-5/JUN 

13:30 

01:29 

MAG062 

18 /OCT 

19:30 

01:09 

19-24/OCT 

13:10 

01:09 

HAG045 

6/JUN 

14:10 

01:29 

7-12/JUN 

13:30 

01:29 

NAG063 

25/OCT 

15:20 

01:09 

26-31/OCT 

13:10 

01:09 

HAG046 

13/JUN 

14:50 

01:29 

14-19/JUN 

13:30 

01:29 

NAG064 

1/NOV 

16:30 

01:09 

2-7/NOV 

13:10 

01:09 

MAG047 

20/JUN 

16:00 

01:29 

21-26/JUN 

13:30 

01:29 

Malmstrom  Air  Force  Base,  Montana 


APPENDIX  A 


lAPE  I 

DATE 

START 

STOP 

MAG067 

18/NOV 

18:50 

01:09 

19/NOV 

13:10 

01:09 

20/NOV 

13:10 

17:52 

NAG068 

25/NOV 

19:10 

00:02 

MAG069 

27/NOV 

22:20 

01:09 

28-30/NOV 

13:10 

01:09 

1/DEC 

13:10 

01:09 

2-3/DEC 

13:20 

01:19 

HAG070 

6/DEC 

23:10 

01:19 

7-12/DEC 

13:20 

01:19 

HAG071 

14/DEC 

15:10 

01:19 

15-20/DEC 

13:20 

01:19 

ArrRNDix  A 


2.  A  Malah.ir  TransmUtor  Amipx,  I’alm  Uay,  I'loridi 


TAPE  i  DATE  START  STOP 


BAR008  8/DEC  17:00  23:19 
9-14/DEC  11:20  23:19 

BAR009  15/DEC  13:40  23:19 
16-21/DEC  11:20  23:19 

BAROlO  22/DEC  14:10  23:19 
23/DEC  11:20  23:19 
24-28/DEC  11:30  23:29 

BAROll  29/DEC  14:30  23:29 
30-31/DEC  11:30  23:29 


YEAR  =  1989 


BAROll  1-4/JAN  11:30  23:29 

BAR012  5/JAN  13:30  23:29 

6/JAN  11:30  18:38 


BAR013  6/JAN 

7-11/JAN 

12/JAN 


20:30  23:29 
11:30  23:29 
11:30  17:21 


BAR014  12/JAN  18:00  23:29 
13-15/JAN  11:30  23:29 
16-18/JAN  11:40  23:39 


BAR015  19-25/JAN  11:40  23:39 


BAR023  15/MAR  13:50  23:39 
16-20/MAR  11:40  23:39 
21/MAR  11:30  23:39 


BAR024  23/MAR  15:50  23:29 
24-29/MAR  11:30  23:29 


BAR025  30/MAR 
31/MAR 
1-5/APR 


13:30  23:29 
11:30  23:29 
11:30  23:29 


BAR026  6/APR 
7/APR 

BAR027  7/APR 

8-10/APR 

11/APR 


12:00  23:29 
11:30  19:01 

19:10  23:29 
11:30  23:29 
11:30  14-29 


TAPJ  iDaiE 

-SlARl 

_ STOP 

BAR028 

25/APR 

15:40 

23:19 

26/APR 

11:20 

23:19 

27/APR 

17:50 

23:19 

28-30/APR 

11:20 

23:19 

1/MAY 

11:20 

23:19 

BAR029 

3/MAY 

11:50 

23:19 

4/MAY 

11:20 

23:19 

BAR030 

5/MAY 

13:10 

23:19 

6-11/HAY 

11:20 

23:19 

BAR031 

12/MAY 

12:30 

23:19 

13-18/HAY 

11:20 

23:19 

BAR032 

19/MAY 

12:30 

23:19 

BAR033 

22/MAY 

11:50 

23:19 

23-28/MAY 

11:20 

23:19 

BAR034 

30/MAY 

21:20 

23:19 

1-5/JUN 

11:20 

23:19 

BAR035 

6/JUN 

12:20 

23:19 

7-12/JUN 

11:30 

23:29 

BAR036 

13/JUN 

13:30 

23:29 

14-19/JUN 

11:30 

23:29 

BAR037 

20/JUN 

12:40 

23:29 

21-26/JUN 

11:30 

23:29 

BAR038 

27/JUN 

18:50 

23:29 

28-30/JUN 

11:30 

23:29 

1-3/JUL 

11:30 

23:29 

BAR039 

5/JUL 

13:10 

23:29 

6-11/JUL 

11:30 

23:29 

BAR040 

12/JUL 

15:40 

23:29 

13-18/JUL 

11:30 

23:29 

BAR041 

19/JUL 

13:30 

23:29 

20-25/JUL 

11:30 

23:29 

BAR042 

26/JUL 

11:30 

23:29 

27/JUL 

11:30 

12:11 

BAR441 

7/AUG 

13:10 

23:29 

8/AUG 

11:30 

11:53 

59 


APPENDIX  A 


Malabar  Transmitter  Annex,  Palm  Bay,  Florida 


TAPE  i  DATE  START  STOP 


TAPE  I  DATE  START  STOP 


BAR443  9/AUG 
10/AUG 


19:30  23:39 
11:30  12:30 


BAR063  15-21/DEC  11:20  23:19 


BAR444  10/AUG  14:50  18:09 

BAR045  10/AUG  18:50  23:29 
11-16/AUG  11:30  23:29 

BAR046  17/AUG  11:50  23:29 
18-23/AUG  11:30  23:29 


BAR064  22/DEC 
23/DEC 
24/DEC 


15:50  23:19 
11:20  23:19 
11:20  12:51 


BAR065  27/DEC  12:30  23:19 
28/DEC  11:20  23:19 
29-31/DEC  11:30  23:29 


BAR047  24-26/AUG  11:30  23:29 
28-30/AUG  11:20  23:19 


YEAR  -  1990 


BAR048  1/SEP  12:00  23:19 
2-7/SEP  11:20  23:19 

BAR049  8/SEP  11:50  23:19 
9-14/SEP  11:20  23:19 

BAR050  15-21/SEP  11:20  23:19 

BAR051  22/SEP  12:00  23:19 
23/SEP  11:20  23:19 
24-28/SEP  11:10  23:09 

BAR052  29-30/SEP  11:10  23:09 
1-5/OCT  11:10  23:09 

BAR053  6/OCT  11:40  23:09 
7-12/OCT  11:10  23:09 

BAR054  20-26/OCT  11:10  23:09 

BAR055  13/OCT  12:10  23:09 
14-19/OCT  11:10  23:09 

BAR056  27-31/OCT  11:10  23:09 
1  2/NOV  11:10  23:09 

BAR057  3-9/NOV  11:10  23:09 

BAR058  10-16/NOV  11:10  23:09 

BAR059  17-23/NOV  11:10  23:09 

BAR060  24-30/NOV  11:10  23:09 

BAR061  1-6/DEC  11:10  23:09 

7/DEC  11:20  23:19 

BAR062  8  14/DEC  11:20  23:19 


BAR065 

1-2/JAN 

11:30 

23:29 

BAR066 

3/JAN 

14:00 

23:29 

4-9/JAN 

11:30 

23:29 

BAR067 

10/JAN 

13:30 

23:29 

11-16/JAN 

11:30 

23:29 

BAR068 

17/JAN 

12:20 

23:29 

18-23/JAN 

11:30 

23:29 

BAR069 

24/JAN 

14:30 

23:29 

25-30/JAN 

11:40 

23:39 

BAR070 

31/JAN 

15:20 

23:39 

1-6/FEB 

11:40 

23:39 

BAR071 

14/FEB 

14:10 

23:39 

BAR072 

7/FEB 

14:00 

23:39 

0-13/FEB 

11:40 

23:39 

BAR073 

15/FEB 

14:50 

23:39 

16-21/FEB 

11:40 

23:39 

BAR074 

22/FEB 

12:40 

23:39 

23-20/FEB 

11:40 

23:39 

BAR075 

1/HAR 

12:10 

23:39 

2-4/HAR 

11:40 

23:39 

5-6/HAR 

11:30 

23:29 

7/MAR 

11:40 

23:29 

BAR076 

8/MAR 

14:20 

23:29 

9-14/HAR 

11:30 

23:29 

BAR077 

1  5/MAR 

14:00 

23:29 

16  21/MAR 

11:30 

23:29 

AIM'I'NDIX  A 


M.il.il);n  I  t  .iiisni  i  I  I  i>r  Annex,  !’n  I  rii  ll.iv,  I’ I  o  r  i  <!.i 


TAPE  1 

date _ 

START 

STOP 

TAPE  1 

DATE 

START 

STOP 

BAR078 

22/MAR 

14:50 

23:29 

BAR095 

17/SEP 

13:50 

23:19 

23-25/MAR 

11:30 

23:29 

18-22/SEP 

11:20 

23:19 

26/MAR 

11:30 

14:00 

23/SEP 

11:10 

23:09 

BAR079 

28/MAR 

15:40 

23:29 

BAR096 

24/SEP 

13:00 

23:09 

29-31/MAR 

11:30 

23:29 

25-30/SEP 

11:10 

23:09 

1-3/APR 

11:30 

23:29 

BAR097 

1/OCT 

11:10 

16:32 

BAR080 

4/APR 

12:10 

23:29 

5-7/APR 

11:30 

23:29 

BAR098 

1/OCT 

18:10 

23:09 

8-10/APR 

11:20 

23:19 

2-7/OCT 

11:10 

23:09 

BAR081 

11/APR 

12:10 

23:19 

BAR099 

8/OCT 

11:20 

23:09 

12-17/APR 

11:20 

23:19 

9/OCT 

11:10 

23:09 

10/OCT 

11:10 

20:32 

BAR082 

18/APR 

11:30 

23:19 

19-24/APR 

11:20 

23:19 

BARIOO 

11/OCT 

12:30 

23:09 

12-17/OCT 

11:10 

23:09 

BAR083 

25/APR 

11:40 

23:19 

26-30/APR 

11:20 

23:19 

BARlOl 

18/OCT 

11:40 

23:09 

1/MAY 

11:20 

23:19 

19-24/OCT 

11:10 

23:09 

BAR084 

6/JUL 

15:00 

23:29 

BAR102 

25-31/OCT 

11:10 

23:09 

7/JUL 

11:30 

13:52 

BAR103 

1/NOV 

11:10 

16:25 

BAR08S 

7/JUL 

15:20 

23:29 

8-13/JUI. 

11:30 

23:29 

BAR104 

16/NOV 

16:10 

23:09 

17-18/NOV 

11:10 

23:09 

BAR086 

15-21/JUL 

11:30 

23:29 

19/NOV 

11:10 

19:56 

BAROe? 

23/JUL 

15:40 

23:29 

24-29/JUL 

11:30 

23:29 

BAR088 

30-31/oJL 

11:30 

23:29 

1-5/AUG 

11:30 

23:29 

BAR089 

6/AUG 

15:00 

23:29 

7-12/AUG 

11:30 

23:29 

BAR090 

13/AUG 

13:10 

23:29 

14-19/AUG 

11:30 

23:29 

BAR091 

20/AUG 

13:10 

23:29 

21-26/AUG 

11:20 

23:19 

BAR092 

27-31/AUG 

11:20 

23:19 

1-2/SEP 

11:20 

23:19 

BAR093 

3-9/SEP 

11:20 

23:19 

BAR094 

10/SEP 

12:40 

23:19 

11-16/SEP 

11:20 

23:19 

A  I 


APPENDIX  A 


2,7  Natloaal  Weather  Facility,  Columbia,  Missouri 


EmTi 

■Q/TyMlHIl 

mm 

YEAR  = 

1989 

COL018 

29-30/JUN 

12:10 

00:09 

1-5/JUL 

12:10 

00:09 

COLOOl 

9/FEB 

16:20 

00:29 

10-14/FEB 

12:30 

00:29 

COL019 

6-11/JUL 

12:10 

00:09 

15/FEB 

13:10 

00:29 

12/JUL 

12:20 

00:19 

COL002 

16-21/FEB 

12:30 

00:29 

COL020 

13-19/JUL 

12:20 

00:19 

22/FEB 

13:10 

00:29 

COLO 21 

20-26/JUL 

12:20 

00:19 

COL003 

23/FEB 

14:00 

00:29 

24-28/FEB 

12:30 

00:29 

COL022 

27/JUL 

20:30 

00:19 

1/HAR 

13:10 

00:29 

28-30/JUL 

12:20 

00:19 

1-2/AUG 

12:20 

00:19 

COLOOL 

2-3/MAR 

12:30 

00:29 

4-7/MAR 

12:20 

00:19 

COL023 

3-9/AUG 

12:20 

00:19 

COL005 

8-14/MAR 

12:20 

00:19 

COLO  2 4 

10-16/AUG 

12:10 

00:09 

COL006 

15-21/MAR 

12:20 

00:19 

COLO 2 5 

17-23/AUG 

12:10 

00:09 

COL007 

22-28/MAR 

12:20 

00:19 

COLO  2 7 

24-30/AUG 

12:10 

00:09 

COLOO0 

29-31/MAR 

12:20 

00:19 

COL028 

31/AUG 

14:10 

00:09 

1-3/APR 

12:20 

00:19 

1-6/SEP 

12:10 

00:09 

COL009 

5/APR 

12:20 

00:19 

COL029 

7/SEP 

15:10 

00:09 

6/APR 

12:40 

00:19 

8-13/SEP 

12:10 

00:09 

7/APR 

12:20 

00:19 

8-11/APR 

12:10 

00:09 

COLO 30 

14/SEP 

14:00 

23:59 

15-20/SEP 

12:00 

23:59 

COLOlO 

12-17/APR 

12:10 

00:09 

18/APR 

12:10 

22:09 

COL031 

21-27/SEP 

12:00 

23:59 

COLOll 

19-25/APR 

12:10 

00:09 

COLO  3 2 

28-30/SEP 

12:00 

23:59 

1-4/OCT 

12:00 

23:59 

COL012 

26/APR 

12:10 

00:09 

27/APR 

12:10 

22:44 

COL033 

5/OCT 

17:30 

23:59 

6-11/OCT 

12:00 

23:59 

COLO  13 

30/APR 

19:40 

00:09 

1  4/MAY 

12:10 

00:09 

COLO  3 4 

12-13/OCT 

12:00 

23:59 

5/MAY 

12:20 

00:09 

14-18/OCT 

11:50 

23:49 

6/MAY 

12:10 

00:09 

COL035 

19-25/OCT 

11:50 

23:49 

COLOH 

16/MAY 

12:10 

20:54 

COLO 36 

26-31/OCT 

11:50 

23:49 

COL015 

16/MAY 

21:20 

00:09 

1/NOV 

11:50 

23:49 

17-22/MAy 

12:10 

00:09 

COL037 

2/NOV 

22:40 

23:49 

COL016 

23-29/MAY 

12:10 

00:09 

3/NOV 

11:50 

16:21 

COLO 17 

22/JUN 

21:20 

00:09 

COL038 

3/NOV 

16:24 

23:49 

2J-28/JUN 

12:10 

00:09 

4-8/NOV 

11:50 

23:49 

M’l’I'.NDtX  A 


Natioiinl  V/onLlior  lATcility,  (lolumbin,  Missouri 


TAPE  i 

_J2AI£ _ 

JilAEIL 

STOP 

IAE£  1  DATfL 

START 

STOP 

COLO  39 

9-15/NOV 

11:50 

23: 49 

COL058 

15-21/MAR 

12:20 

00:19 

COL040 

16-17/NOV 

11:50 

23:49 

COL059 

22-25/MAR 

12:20 

00:19 

18-22/NOV 

12:00 

23:59 

26-28/MAR 

12:10 

00:09 

COL041 

23-29/NOV 

12:00 

23:59 

COL060 

2/APR 

12:50 

00:09 

3-8/APR 

12:10 

00:09 

COLO 4 2 

30/NOV 

12:00 

23:59 

1-6/DEC 

12:00 

23:59 

COL061 

9-11/APR 

12:10 

00:09 

COLO  4 3 

14/DEC 

12:23 

23:59 

COL062 

12-18/APR 

12:10 

00:09 

COL044 

15/DEC 

12:00 

23:59 

COL063 

19-25/APR 

12:10 

00:09 

16-21/DEC 

12:10 

00:09 

COL064 

26-30/APR 

12:10 

00:09 

COL045 

22/DEC 

14:30 

00:09 

1-2/MAY 

12:10 

00:09 

23-27/DEC 

12:10 

00:09 

COL065 

3-9/MAY 

12:10 

00:09 

COLO  4 6 

28-31/DEC 

12:10 

00:09 

COL066 

10-16/MAY 

12:10 

00:09 

YEAR  = 

1990 

COL067 

17-23/MAY 

12:10 

00:09 

COLO 6 8 

24-30/MAY 

12:10 

00:09 

COL046 

1-3/JAN 

12:10 

00:09 

COL069 

31/MAY 

12:10 

00:09 

COLO  4 7 

4/JAN 

13:30 

00:09 

1-4/JUN 

12:10 

00:09 

5-6/JAN 

12:10 

00:09 

7-10/JAN 

12:20 

00:19 

COLO 70 

5/JUN 

13:50 

00:09 

6/JUN 

12:10 

00:09 

COL048 

11-17/JAN 

12:20 

00:19 

COL071 

7- 13/JUN 

12:10 

00:09 

COLO  4 9 

18-24/JAN 

12:20 

00:19 

COLO 7 5 

1/AUG 

17:50 

00:19 

COL050 

25-31/JAN 

12:20 

00:19 

2-7/AUG 

12:20 

00:19 

COL051 

1-7/FEB 

12:20 

00:19 

COL076 

8/AUG 

12:20 

00:19 

9-14/AUG 

12:10 

00:09 

COL052 

8-14/FEB 

12:20 

00:19 

COL077 

15-21/AUG 

12:10 

00:09 

COL053 

15/FEB 

15:16 

00:19 

16-21/FEB 

12:20 

00:19 

COLO 7 8 

22-28/AUG 

12:10 

00:09 

COL054 

22-28/FEB 

12:20 

00:19 

COLO 7 9 

29-31/AUG 

12:10 

00:09 

1-4/SEP 

12:10 

00:09 

COL055 

1-7/MAR 

12:20 

00:19 

COL080 

5-11/SEP 

12:10 

00:09 

COL056 

8-10/MAR 

12:20 

00:19 

11/MAR 

12:20 

00:51 

COLO 81 

19-25/SEP 

12:00 

23:59 

COL057 

12/MAR 

20:24 

20:24 

COL082 

26-30/SEP 

12:00 

23:59 

13-14/MAR 

12:20 

00:18 

1-2/OCT 

12:00 

23:59 

A3 
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National  Weather  Facility,  Columbia,  Missouri 


TAPE  1 

DATE 

START 

STOP 

COLO 8 3 

3-9/OCT 

12:00 

23:59 

COLO 8 4 

10-12/OCT 

12:00 

23:59 

13-16/OCT 

11:50 

23:49 

COLO 05 

17-23/OCT 

11:50 

23:49 

COLO06 

24-30/OCT 

11:50 

23:49 

COLO 8 7 

31/OCT 

11:50 

23:49 

1-6/NOV 

11:50 

23:49 

COL088 

7-13/NOV 

11:50 

23:49 

COL089 

14-16/NOV 

11:50 

23:49 

17-20/NOV 

12:00 

23:59 

COL090 

21-27/NOV 

12:00 

23:59 

COL091 

28/NOV 

19:40 

23:59 

29-30/NOV 

12:00 

23:59 

1-4/DEC 

12:00 

23:59 

COL092 

10/DEC 

J7;10 

23:59 

11-14/DEC 

12:00 

23:59 

15-16/DEC 

12:10 

00:09 

COLO 9 3 

17-23/DEC 

12:10 

00:09 

COL094 

24-30/DEC 

12:10 

00:09 

COL095 

31/DEC 

12:10 

00:09 

OPTICAL  SYSTEMS  GROUP 
TECHNICAL  NOTE  NO.  232 
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Bl.O  INTRODUCTION 

This  appendix  provides  a  brief  overview  of  the  soft¬ 
ware  used  to  produce  cloud  decision  imagery  from  the 
raw  radiance  imagery  collected  in  the  field.  An  overview 
of  the  procedure  used  to  produce  the  1 -minute  and  10- 
minute  cloud  decisions  for  the  contract  is  shown  in  Fig. 
Bl-1.  The  following  three  sections  describe  the  three 
main  processing  programs,  TAPRATPL,  CLDDECM 
and  CMPDECTP.  Section  B5  describes  some  of  the 
other  programs  used  in  intermediate  steps.  Section  B6 
lists  brief  descriptions  of  utility  routines  written  at  MPL 
that  are  accessed  by  the  three  major  processing  routines. 

B2.0  TAPRATPL  -  The  Composite  Ratio  Process¬ 
ing  Program 

TAPRATPL  is  the  program  that  takes  a  set  of  four 
radiance  images  from  a  field  tape,  and  produces  a  com¬ 
posite  ratio  image.  The  ratio  image,  along  with  its 
corresponding  dark  red  image,  is  then  copied  onto  an 
output  tape  that  is  usually  referred  to  as  a  ratio  tape. 
TAPRATPL  is  run  either  for  1  or  10  minute  images. 
Thus,  a  given  field  tape  must  be  processed  twice  if  both 
1  and  10  minute  ratio  images  are  needed. 

To  run  TAPRATPL,  the  AT  type  computer  must  be 
equipped  with  an  IT!  FG-KX),  a  PL-1250  accelerator 
board,  and  2  ExaByte  tape  drives  (one  for  input,  the  other 
for  output).  The  PL  board  allowed  us  to  cut  the  compu¬ 
tation  time  by  a  factor  greater  than  3,  enabling  a 
TAPRATPL  run  to  usually  take  less  than  1 2  hours.  It  also 
increased  the  complexity  ofthc  programming  effort.  The 
following  description  provides  a  general  overview  of  the 
TAPRATPL  main  program. 

The  input  files  needed  by  TAPRATPL  are  described 
in  Tech  Notes  230  and  23 1 .  The  program  also  requires 
user  input  regarding  the  field  tape  number,  the  type  of 
images  being  proces.scd,  thin  and  opaque  values  for  a 
color  ratio  display  option,  the  site  number  for  the  field 
tape,  and  the  expected  start  date  on  the  input  tape. 

The  program  makes  several  consistency  checks  be¬ 
tween  input  parameters  cither  from  the  user  or  the  input 
files,  and  the  input  tape,  and  notifies  the  user  of  any 
di.scrcpancy.  Other  parts  of  the  code  will  check  for 
consistency  witliin  the  four  image  sets,  and  the  main 
program  also  notifies  the  user  of  any  unexpected  time 
shifts.  All  peculiarities  are  written  to  an  output  diagnos¬ 
tic  file.  A  log  file  containing  date,  time,  quality  and  ratio 
histogram  information  is  also  produced  by  TAPRATPL. 

figure  B2-1  .shows  a  conceptual  flowchart  for  the 
I’Ai’RATPL  program.  Basically,  after  some 
initializations  and  input  definition  steps,  the  program 
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consists  of  an  outer  day  loop  and  an  itmer  image  loop. 
Due  to  the  occasional  date  shifts  in  our  field  systems,  date 
changes  are  detected  by  backward  time  skips  of  greater 
than  four  hours. 

B2.I  TAPRATPL  Subroutine  Descriptions 

TAPRATPL  is  a  complex  program,  linking  together 
more  than  90  subprograms.  The  names  of  the  routines 
developed  here  at  MPL  for  TAPRATPL  are  listed  in 
Figure  B2-2.  The  remaining  MPL  utility  programs  are 
described  in  a  later  section.  Additional  subprograms 
perform  specific  functions  on  the  FG-100  and  PL-1250 
boards,  and  are  in  documentation  provided  by  the  ven¬ 
dors  (ITl  and  Eighteen  Eight  respectively.) 

The  subroutines  written  specifically  for  TAPRATPL 
are  described  below  in  their  order  of  appearance.  The  * 
symbol  indicates  those  routines  that  require  PL- 1250 
access. 

SETPL*  -  initializes  the  PL- 1250  processor,  allocates 
PL  memory  locations  for  subsequent  programs,  and 
transfers  some  of  the  constants  used  in  later  calcula¬ 
tions  to  the  PL  board. 

CLRCLD  -  creates  the  color  lookup  tables  for  the  color 
ratio  display  hotkey  option  from  input  thin  and 
opaque  thresholds. 

CALSEL  -  requests  the  site  number,  opens  the  diagnos¬ 
tic  file,  requests  the  start  date  of  the  field  tape. 
Confirms  the  date  and  field  tape  .sequence  number 
with  data  from  D  ATE_.LST.  Retrieves  the  hard  ware 
version  numbers  and  data  quality  information  from 
VERSION.LOG  for  each  date  on  the  input  tape. 

FRSTHD  -  reads  initial  header  from  field  tape,  end 
determines  the  software  version  that  created  the  tape 
and  whether  the  first  image  is  a  full  resolution,  10 
minute  image, ora. subset  1  minute  image.  Confinns 
the  initial  date  from  CALSEL  input  and  the  site. 
Reads  the  iris,  occultor  and  neutral  density  settings 
for  the  image. 

RDCAL*  -  reads  the  calibration  file  (CAL_V _ .RAT) 

specified  by  CALSEL  (or  a  given  date.  Performs 
some  preliminary  compulations  related  to  the  image 
offset  and  magnifications,  and  the  ratio  multiplier 
factors.  Transfers  this  information  to  the  PL- 1250. 

EQTIME  -  using  code  similarlo  that  u.sed  in  the  FIELD 
program,  compute  the  start  time  (6  hours  before 
LAN)  forthcspccified  Julian  day  and  longitude.  The 
start  time  is  used  to  determine  the  expected  occultor 
position. 
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TIMCOR  -  reads  the  time  correction  from  the 
TIME_.COR  file  for  die  date  specified. 

ARMSET  -  reads  the  occultor  arm  number  and  angle 
offset  from  the  OCC_.DAT  file  for  the  date  speci¬ 
fied. 

GETIMG  -  transfers  the  remainderofthe  field  image  set 
from  ExaByte  to  the  EG- 100  board. 

CHCKHD  -  performs  consistency  checks  on  the  neutral 
density  and  spectral  filtcrdesignations,  and  the  quad¬ 
rant  number  specified  in  all  four  headers  in  an  image 
set. 

KEYCK  -  checks  for  new  hotkey  entries  during  process¬ 
ing.  Activates  flags  that  are  acted  upon  in  the  main 
program. 

LINDIG  -  applies  the  linearity  correction  to  the  radiance 
byte  values  via  an  input  LUT  on  the  FG-100. 

RATCOM*  -  transfers  the  radiance  information  from 
the  EG- 1 00  to  the  PL- 1 250  using  a  double  buffering 
technique  described  in  Eighteen  Eight  (1988).  Calls 
RATPL  which  performs  the  ratio  operations  on  the 
set  of  4  lines  transferred.  Returns  the  resulting 
composite  ratio  line  to  the  FG-100  board. 

RATPL*  -  performs  the  ratio  operations  on  the  4 
image  lines  provided.  Includes  setting  offscale 
dark  points  on  the  bright  red-blue  pair  to  0, 
offscale  bright  points  on  the  dark  red-blue  pair  to 
240,  and  the  remaining  ratio  byte  values  between 
1  and  239. 

OCCMSK  -  determines  the  arm  and  edge  locations  in  3- 
dimensions.  Calls  POSIMA  to  locate  points  on  the 
edge  in  image  coordinates,  and  DRAWLN  to  iden¬ 
tify  the  pixel  locations  connecting  consecutive  edge 
points.  Overlays  the  resulting  occultormask  onto  the 
ratio  image  with  byte  value  0.  Currently,  the  arm 
mask  is  drawn  twice  its  actual  diameter,  and  the 
frame  is  expanded  by  15%  to  compensate  for  me¬ 
chanical  play  in  the  occultor  assembly. 

POSIMA  -  converts  the  3-dimension  location  of  a 
point  on  the  occultor  assembly  to  ratio  image 
coordinates. 

DRAWLN  -  determines  the  image  pixels  along  a 
line  segment  defined  by  any  two  points  in  the 
image  coordinates. 

OBDRAW  -  overlays  the  obstacle  masks  onto  the  ratio 
image  with  byte  value  0  using  edges  defined  in 
RDCAL 


SMOHST*  -  performs  two  passes  of  a  smoothing  filter 
over  the  raw  ratio  values.  Passes  a  double  buffered 
set  of  lines  from  the  FG-100  to  the  PL- 1250.  Calls 
SMTHPL  to  perform  the  smoothing  operations,  then 
returns  the  smoothed  results  to  the  FG-100.  Com¬ 
putes  the  resulting  ratio  histogram  distribution. 

SMTHPL*  -  passes  a  2-dimensional  kernel  over  a 
full  resolution  10-minute  image,  or  a  1 -dimen¬ 
sion  kernel  over  the  subset  lines  of  a  1 -minute 
image.  If  more  than  half  the  weight  of  the  filter 
is  at  a  pi  xel  value  of  0  or  240  (offscale  bright),  no 
smoothing  is  performed  at  that  point. 

WRITHD  -  imbeds  the  ratio  header  information  onto  the 
ratio  image  on  the  FG-100  board. 

PUTIMG  -  writes  a  DOS  header  block  for  this  image 
onto  the  output  ExaByte  drive,  then  copies  the  ratio- 
dark  red  image  pair  from  the  FG-100  to  ExaByte. 

LOOP  -  loops  from  the  colored  ratio  image  to  the  red 
image  when  a  certain  hotkey  is  activated. 

NEXTHD  -  transfer  the  next  header  on  the  input 
ExaByte  to  the  FG-100,  determines  the  image  type, 
and  extracts  the  iris,  occultor  and  neutral  deasity 
settings. 


B3.0  CLDDECM  -  The  Fixed  Threshold  Cloud 
Decision  Program 

In  early  1990,  we  were  tasked  to  provide  preliminary 
cloud-no  cloud  decision  data  for4  sites  forupto  a  yeario 
The  Analytic  Sciences  Corporation  (TASC).  The  cloud 
decisions  for  the  1  -minute  subset  imagery  processed  for 
this  requirement  were  run  using  tJte  fixed  threshold 
algorithm  named  CLDDECM.  This  program  is  much 
simpler  than  TAPRATPL  and  subsequently  much  faster. 
The  cloud  decision  is  made  by  passing  the  ratio  image 
through  an  input  look-up  table  and  no  PL- 1250  compu¬ 
tations  are  required.  The  thin  and  opaque  decision 
thresholds  arc  constant  for  the  entire  tape,  and  arc  speci¬ 
fied  as  user  inputs.  In  addition  to  the  output  tape,  one 
output  file  is  generated,  and  it  contains  sky  cover  results 
only  from  those  minutes  evenly  divisible  by  10. 

The  conceptual  flow  chart  for  CLDDECM  is  shown  in 
Fig.  B3-L  While  not  shown  in  the  diagram,  cnd-of-file 
marks  are  written  to  the  output  ExaByte  drive  at  the  end 
of  each  day.  No  dummy  DOS  headers  arc  placed  before 
the  images,  as  done  in  TAPRATPL. 
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B3.1  CLDDECM  Subroutines 

CLDDECM  links  together  more  than  30  subpro¬ 
grams.  The  7  subroutines  written  specincally  for 
CLDDECM  arc  dc-scribcd  below.  Twelve  MPL  utility- 
programs  are  also  used,  along  with  ,10  Image  Tool 
Library  routines,  (f-'ig.  3-2) 

CLRDEC  -  requests  Uie  thin  and  opaque  threshold 
inputs  and  computes  the  cloud  decision  LUT.  if  the 
ratio  input  value  is  0,  its  corres[XMiding  cloud  deci¬ 
sion  value  will  also  be  /ero.  Ratio  values  between  1 
and  the  thin  threshold  inclusive  will  lx;  assigned  the 
cleardeci.sion  value  of  1  (K).  The  remaining  values  up 
to  and  including  the  opaque  threshold  will  be  as¬ 
signed  the  thin  cloud  decision  value  of  150.  Those 
points  with  ratios  greater  than  the  opaque  threshold 
up  to  die  value  239  will  be  given  the  opaque  cloud 
value  of  200.  Finally,  the  off-scale  bright  points, 
with  a  ratio  value  of  240,  w  ill  be  assigned  a  cloud 
decision  value  of  202.  The  color  LUT  for  hotkey 
cloud  decision  displays  are  also  determined. 

CUKTAP  -  tc.sls  of  the  source  and  target  ExaByte  drives 
arc  ready,  and  can  write  an  EOF  to  the  target  ExaByte. 

1-RATHD  -  finds  the  DOS  header  and  imbedded  ratio 
licadcrs  in  die  ratio  ta[X'  images. 

ONfiRAT  -  reads  the  remaining  65  hues  ol  the  1 -minute 
ratio  images. 

Cl  TIKE  Y  -  checks  for  hotkey  ciiiries  during  processing, 
and  .sets  appropriate  (lags 

WRCDUD  -  writes  the  cloud  decision  tieader  informa¬ 
tion  to  the  f-Ci  KX). 

ONESAV  -  copies  the  cloud  decision  image  from  the 
FG-KK)  to  the  target  ExaByic 

B4.0  CMPDECI  P  -  The  ( 'oniposite  <  ieud 
Decision  Program 

<  lurcxpericncc  with  the  cloud  dei  is.nrw  produ'.'-d  b\ 
I'lxed  threshold  algoritliins  such  ,i  -  Cl.l  M )l  C'\1,  w:!-  ih  it 
while  opa(]uccl{Mjdwcn’  reasonabls  well  id'  otificd.lhln 
cloud  in  the  dt-iwnsun  region  of  the  sky  weie  urideier- 
mincd,  while  regions  near  the  tion/on,  pardcularly  at 
near  hori/.on  sun  .ingles  were  ineorrei’ll v  iden'ified  as 
cloudy.  To  overcome  these  problems,  .i  new  algondiin 
was  developed.  fixed  threshold  is  still  uw'd  to  identily 
opaque  clouds,  but  the  thin  cloml  delemmiaiion  is  now- 
made  against  an  estimate  of  the  clear  sky  ratio  that  v.ines 
in  both  space  and  lime.  A  more  detailed  desenirtion  of 
the  algonthm  is  provided  in  .Section  5  in  the  mam  Ixxly 
of  the  report 


C.MPDHCPL  also  requires  the  PL- 1250  board  to  get 
the  nin  ume  lor  a  7  day  tape  down  to  about  3  1/2  hours. 
(Without  the  FI  ,  the  run  time  would  exceed  15  hours). 
For  .1  paiti.-ii-ar  fio!  ;  ta;x-.  three  input  files  must  be 
pR'ser.i  CM!';  -t T'  pni'  -,p,-eifies  the  opaque  thawnolds 
for  each  o!  the  a  ,;al;..jl  density  .sellings,  the  a/imulh 
offset,  acsi  ifte  niiu  cf  ud  ,r.  \  i,  piance  level  for  a  particular 
hardware  version  Auoilicr  file,  SKY.D.'kT 

contains  li.e  nonua'i  -ed  Reii.-'Hlue  clear  sky  radiance 
e.'timair  s  .j;  5  '  meii-m.  nis  lO  solar  /enilh  angle  that  were 
extracted  Irom  clear  day  samples  for  the  site.  The  first  3 
letters  in  this  file  name  an:  the  1  letter  site  identifier,  (For 
example,  COLSKY  DA  I  would  contain  the  Columbia, 
MO  tables  )  The  final  file  needed  is  _  _  _  ###.REF, 
containing  the  refeamice  factors  that  vary  as  a  function  of 
boundary-  layer  ha/.e  loading.  File  KAA050.REF  would 
contain  the  releamcc  values  for  tape  number  50  from 
Kinland  AFB,  NM.  These  files  arc  generated  in  the 
Sin'RFF  program  described  in  the  next  section. 

A  conceptual  flo-Achart  for  CMPDECTP  is  shown  in 
Fig.  B4- 1  CMPDECTP  only  processes  the  full  resolu¬ 
tion  imagery  taken  „l  lO-minuie  intervals  in  the  field,  in 
contrast  to  the  1 -minute  imagery  processed  by 
CLDDTiC.M.  A  1  -minute  v,;rsion  of  the  new  algorithm 
has  yet  to  Iv  dcvcio|x  J  As  m  T.APRATPL,  dummy 
DOS  hetidots  ;;[u  v  riiten  to  tlie  output  ExaByte  before 
each  cloud  decision  image 

B4.1  -  C  MPDF.f  !  P ‘‘TibroLitine  Description 

I  he  .subrt'-itm;".  c  ,'ded  to  tun  ('MPDECTP  are  de¬ 
scribed  in  Eig  H4  1  IX'scnpiions  for  those  programs 
w'ntten  s|x-citicall\  lor  the  ('MPDECT’P  routine  arc 
pros  Ided  belo'vx  Ih.ose  Mjbn'ulines  interacting  with  the 
PL- 1250  fx-aid  are  again  marked  witli  an  *. 

INfECI)  -  mitiali/e  the  I  (!- liX!  f^xiard 

SI  I  i’l  ( '■'  '.■Klip;  T'!  .'-.Mn '.vcs  and  pas.ses  some  of 
Ihe  const. tiu  .md  lived  '.  alue  arrays  to  the  PL- 
:  ^so 

Ci  IK'.ME  rm-'  r  .  ol.i’  :  2,  "f  tor  elouti  decision 
im;  d: -;-i  e. 

f'l'R'l  Hi  )  rea;*'  td.:  nci  mHo  header  from  the  input 
f,-'  (Hyn 

P.-MvCMf'  on  idic  hi-.;  time  ...nieil.  sets  the  ratio  image 
si/e  [viiamciiTK  .C'mt'un's  l.mtude  and  longitude 
Irom  111'.-  r.ttio  lieidc  mlnrmalion  and  reads  tlic 
1  MPDi'C  INI’  ti!;.  and  *  ki  E  file  forthis  tape.  For 
e.ieh  new  d.iie  ■,  ncoiinieivd.  ihe  solar  declination 
angle,  needed  lo  .  on'iniie  sol.ii  /enith  and  a/imuth, 
is  (leiermTied 
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RDTRAT  -  reads  the  remaining  479  lines  of  the  10-min 
ratio  image. 

CMPCMT  -  computes  the  solar  zenith  and  azimuth 
angles  forthe  image,  and  interpolates  the  normalized 
clear  sky  reference  tables  to  the  image  zenith.  Inter¬ 
polates  the  haze-related  reference  factor  in  time  to 
image  time.  Passes  upper  and  lower  ratio  image  row 
pairs  to  the  PL  for  processing,  then  places  the  result¬ 
ing  cloud  decision  rows  into  the  FG- 1 00.  At  comple¬ 
tion,  computes  sky  cover  parameters. 

CMPPLT  -  interpolates  the  sky  ratio  table  to  the 
image  row  pixel  locations,  and  multiplies  by  the 
reference  value  producing  a  clear  sky  pixel  esti¬ 
mate.  A  cloud  decision  value  is  then  assigned  at 
each  pixel  in  the  following  manner. 

No  data:  RATIO  =  0 

DECIS  =  0 

Off-Scale  bright:  RATIO  =  240 

DECIS  =  201 

Indeterminate:  CLEAR  >  OPAQUE 

DECIS  =  202 

Opaque  Cloud:  RATIO  >  OPAQUE 

DECIS  =  100  (RATIO/OPAQUE) +  40 
Thin  Cloud:  RATIO/CLEAR  >  ACCEPT 

DECIS  =  40  [RAT10/(ACCEPT  •  CLEAR)|+60 
Clear:  RATIO/CLEAR  <  ACCEPT 

DECIS  =  100  (RATIO/OPAQUE) 

where  RATIO  -  ratio  byte  value 

CLEAR  -  clear  sky  ratio  estimate 


CLDDECM.  Following  the  cloud  decision  in 
subroutine  CMPPLT,  the  decision  histogram  is 
then  updated. 

SKYINT  -  first  time  through  loads  the  normalized 
clear  sky  rcd/blue  ratio  estimates  for  each  5° 
.solar  zenith  angle.  Then  interpolates  the  tables 
to  the  observed  solar  zenith  angle  (all  calls). 

SUNPOS  -  determines  the  sun  position  in  WSI 
zenith  and  azimuth  space  for  the  .specified  time 
and  azimuth  offset. 

REFCOM  -  interpolates  the  reference  values  forthe 
date  specified  to  the  time  specified. 

CLDKEY  -  checks  for  hot  key  entries  during  processing 
and  sets  appropriate  flags. 

WRCTHD  -  writes  the  cloud  decision  header  informa¬ 
tion  to  the  FG-1(X)  board. 

CMPSVT  -  writes  a  DOS  header  block  and  decision 
image  on  the  output  ExaByte  drive. 

B5.0  OTHER  RELATED  ROUTINES 

Several  computer  routines  have  been  written  by  our 
group  for  WSI  data  display  and  manipulations.  The 
routines  needed  to  perform  many  of  the  intermediate 
tasks  in  the  ratio  and  cloud  decision  processing  sequence 
are  described  in  this  section. 

CHECKAZM  -  determines  the  azimuth  offset  of  the 
sensor  with  respect  to  true  North.  Done  fora  particu¬ 
lar  hard  wane  version  from  a  raw  radiance  tape.  Houiiy 
sun  positioas  for  an  entire  clear  day  are  extracted 
from  full-resolution  image  with  the  aid  of  a  mouse. 
These  positions  are  then  plotted  on  a  nominal  WSI 
image.  The  solar  path  fora  user-input  azimuth  offset 
is  then  plotted.  The  offset  can  then  be  changed  until 
the  best  visual  fit  to  the  extracted  sun  positions  is 
found. 


OPAQUE  -  opaque  ratio  threshold 

ACCEPT  -  thin  cloud  acceptance  level 

DECIS  -  cloud  decision  byte  value 

Thus,  clear  pixels  fall  between  1  and  99,  thin 
cloud  pixels  between  100  and  139,  and  opaque 
cloud  pixels  between  140  to  200  in  the  compos¬ 
ite  decision  imagery.  Some  of  the  original  vari¬ 
ability  in  the  ratio  image  remains  in  the  decision 
image,  in  contrast  to  the  CLDDECM  output  that 
had  single  byte  values  for  each  category.  The 
indeterminate  category  also  was  not  available  in 


CHECKLAN  -  allows  the  u.serto  determine  the  occultor 
arm  and  offset  values,  and  the  time  offset  for  a  LAN 
image  on  disk.  Other  quality-related  comments 
regarding  the  image,  and  estimated  cloud  cover 
values  are  also  manually  recorded  during  these  runs. 
This  information  forms  the  basis  for  several  of  the 
TAPRATPL  input  files. 

NDRCALC  -  computes  the  rcdA)lue  ratios  from  raw 
radiance  quadrants  1  and  2,  and  then  from  quadrants 
3  and  4.  Then  forms  the  ratio  between  these  values, 
and  writes  the  histogram  of  the  final  ratio  pixel 
values.  The  peak  in  this  histogram,  along  with  the 
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absolute  calibration  results,  provide  a  good  first 
cstimatcofthe  NDR  correction  factorinTAPRATPL. 

OBSBETA  -  using  measured  nitio  images,  computes 
normalized  clear  sky  ratio  distribution  from  fuil 
resolution  ratio  tape  input,  usually  from  completely 
clear  days.  Using  a  mouse,  nonreprescniativc  points 
in  regions  of  stray  light  or  covered  h\  cloud  elemer.ts 
can  be  flagged  for  removal  from  the  sample  (5^ 
zenith  and  15°  relative  azimuth  mtcrvals)  extracted 
Irom  the  full  image. 

COMBNORM  -  combines  the  nonnahzed  clear  sky 
ratio  flics  from  OBSBETA  for  a  given  site.  Tables 
from  the  same  solar  zenith  angle  arc  averaged  to¬ 
gether  for  all  the  days  with  OBSBETA  output.  The 
resulting  averaged  tables  arc  examined  for  consis¬ 
tency,  and  after  some  manual  adjustments,  provide 

the _ SK  /.DAT  files  needed  by  CiviPDECTP  and 

GETREF. 

GETREF  -  determines  a  representative  dear  sky  refer¬ 
ence  value  estimate  for  each  ratio  image  on  a  10- 
minute  rat  0  tape.  TTic  observed  ratio  is  divided  by 
the  norma'  iz.ed  clear  sky  ratio  for  all  points  within  die 
70°  zenith  angle  circle,  and  a  histogram  of  the  result¬ 
ing  value  •'hosen.  A  seven  count  windov,  is  passed 
across  the  low  value  end  of  the  histogram  to  find  the 
window  \.ith  die  largest  image  area  contribution. 
The  cent!  "ol  the  maximum  wnuU'w  and  the  area  of 
the  wind('-‘''s contribution  are  written  tothcGF  fREF 
output  fi! '. 

SETREF  -  pi  Its  the  GETREF  reference  values  for  up  to 
four  days  in  different  colors,  and  alUnss  the  user  to 
select  rt.  rcsentativc  reterence  values  for  the  day 
Value  selection  is  performed  widi  a  mouse.  Results 
arc  writfi  n  to  *.REF  files  for  use  in  '  ’MFIiHGl'P. 

RA'fRFiF  -  ■  e.ids  a  ratio  iin.iU  ' Irom  ,1  lO  tnar.ii i  itio 
tape  and  t)nslrucl^  a  raiio  hisioL.a.m  iV.  use  oi  a 
mouse,  tiv  User  can  then '  d'llc'e...  tl"  'sh-  ds 
The  region  of  the  image  b.'liw,  n.,u  ihre-.b.  ■!,.!  is 
cok'red  blue  on  the  ratio  image  d  h  ■  ■'evuliing  ’•ano 
image  e.in  lIk'ii  be  compared  tr*  d  .■  ci.i-t,  i.-d  unagi 
that  accompanies  Ib.e  ratio  imj.ee  oii  i.ipe  1  be  is.i.d 
compansons  enable  the  selection  I'fop.uju  ■  levels  ,n 
different  neutral  dciisit)  setimgs  needed  by 
CMF’DEfTB 

F.XCOBY  -  hacks  up  the  \V, SI  image  !.i|x->  eol'.ci  mv., 
ratio  or  cloud  decision 

EXQC  -  provides  a  summary  of  die  images  av.iilable  on 
a  I -minute  cloud  decision  tape 


ISOTES'f  -  allows  examination  of  different  thin  and 
opaque  threshold  combinations  on  the  images  stored 
on  disk.  Used  primarily  for  CLDDECM  threshold 
dcicrmin.gion 

Se\  er,t!  niu-aiis  ■  f.as  c  been  written  for  viewing  WSl 
imagery  i.'c.m  disi;  and  evaluating  them  with  a  mouse. 
These  include  SHOiMG.  VIEWFMG  and  SAMPIMG. 
Ar.oihcr  set  of  routines  manipulate  the  WSl  images  on 
la[v  include  !,M;  IN'D.  RDIMGS  and  WRIMGS. 

B6.0  .MPL  Utility  Routines  Used  in  Ratio  and 
Cloud  Decision  Processing 

Several  utility  routines  have  been  written  for  use  with 
our  processing  programs.  Brief  descriptions  of  these 
routines  arc  given  below, 

BEEP  -  causes  an  audible  sound  to  be  made  by  the 
computer. 

CALEND  -  computes  the  month,  day  and  year  from  the 
relative  Julian  date,  and  constructs  the  appropriate 
date  siring. 

Cl.RSt.'R  clears  .all  4  quadrants  of  the  EG- 100  board. 

DEF.W  -  p.iu'.’s  tlie  nmsircam  for  the  specified  number 
of  secends 

EOFCT  \x  n'c-  .;  1  Ft  ilM  l  AuByie  unit  1,  (Usually  the 
wriie  dii  .ci. 

!  .\SNS0  -  decodes  the  Ex  iBy  tc  unit  0  sense  data. 

EXSNSl  -  decodes  ihc  ExaByie  unit  1  sense  data. 

FA'’i-G  1  -  p  !vs  s  a  I  K  bkKk  liom  ExaByie  to  the  FG- 
10(1  beard. 

F'GKI-GS  ‘--Is  ihf'  i-fi  ioo  Tvgistcrs  before  data  is 
p,''.-  '.  1  voceii  en  inpu:  i  I'T. 

i '  •SI,  J  s(  : ,  ;he  f  '  '  b  U  ,'cgisiers  for  normal  ExaByie 

Hii'm  i  Aiu-.’s  M.'.  l,  ii'MU  Ihc  FGKKl board  to 
the  i  '  F I  uw  i 

l-Gi-.:"  "i  s  ’-j  ■ ' '  'K  b;  •.  k  tn'-m  the  EG  -  1(X)  board  to 

;  ASfi  !■  Jill  i 

I  1  11  i.  .  •eai"' dun  iiiy  DOS  ficadcr  block  for  a  file 
oi  q Vv  liieu  si  ■  ■ 

I  II  I  ..  reaies  j  J,  a  iio  F)(  ).'k  header  for  a  ratio  image 
loruit’.cr  a  1  or  i;i  mm  miagc. 

INI  I  FG  -  miiiali/cs  the  1  O  fiX)  board 
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JULIAN  -  computes  ihc  relative  Julian  date  from  cither 
a  date  string  or  integer  month,  day  and  year  values. 
The  relative  Julian  date  for  l/JAN/88  is  1. 

OUTEOF  -  writes  an  EOF  to  the  ExaByte  write  drive. 

READY  -  determines  whether  the  ExaByte  unit  speci¬ 
fied  is  ready  for  a  command. 

RSENSE  -  decodes  the  ExaByte  sense  data  from  the  read 
drive. 

SKBLK  -  skips  the  specified  numberof  IK  blocks  on  the 
read  Exabyte. 

SUN  -  compute  the  solar  declination  angle,  radius 
vector,  and  equation  of  time  for  a  particular  longi¬ 
tude  and  date.  Based  on  the  cphcmcris  algorithm  of 
Wilson  (1980). 

WHBLK  -  writes  the  1 K  DOS  header  to  ExaByte. 

WRFMK  -  write  an  EOF  to  ExaByte. 

WSENSE  -  decodes  the  ExaByte  .sense  data  from  the 
write  Exabyte. 
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Figure  B1-1 
Overview  of  Ratio 
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Processing  Steps 


Predeployment 
linearity,  absolute 
and  geometric  calibrations 


Raw  radiance 
field  tapes 


Composite  ratio 

processing 

(TAPRATPL) 

z 

_ : 

1  _ 

Ratio  input  file 
preparation  steps 


Ratio  processing 
input  files 


1 -minute,  subset 
resolution  ratio  tapes 


> 

r 

Fixed  thin  and 
opaque  threshold 
determination 

1 

f 

Fixed  threshold 
cloud  decision 
processing 
(CLDDECM) 

> 

10-minute,  full 
resolution  ratio  tapes 


Fixed  opaque 
threshold 
determination 


Clear  sky  ratio 
parameter 
definitions 


Composite  cloud 
decision  input 
files 


1  minute  cloud 
decision  tapes 


Composite  cloud 
decision  processing 
(CMPDECTP) 


10-min.  cloud 
decision  tapes 


72 


Fig.  B2-1  TAPRATPL  Conceptual  Flowchart 
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TAPRATPL- 


Fig.  B2-2  TAPRAT  PL  Subroutines 
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Fig.  B3'1  CLDDECM  Conceptual  Flowchart 
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Fig.  B3-2  CLDDECM  Subroutines 
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Fig.  B4-1  CMPDECTP  Conceptual  Flowchart 
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CMPDECTP  - 


Fig.  B4-2  CMPDECTP  Subroutines 
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